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SECTION  1 


INTRODUCTION 

I 


1.1  BACKGROUND 

Several  recent  studies  have  demonstrated  that  the  structural 
efficiency  of  future  USAF  aircraft  can  be  improved  by  taking  advantage  of  the 
postbuckled  strength  of  stiffened  panels.  A  majority  of  these  studies,  per¬ 
formed  under  the  auspices  of  the  USAF,  the  U.S.  Navy  and  NASA,  have  provided 
a  sizeable  data  base  on  the  static  and  fatigue  behavior  of  postbuckled  stiff¬ 
ened  panels.  The  test  data  from  the  experimental  studies  have  been  used  to 
verify  the  predominantly  empirical  design  methodology  for  postbuckled  panels 
and  to  establish  static  and  fatigue  failure  characteristics  of  metal  and 
composite  panels.  The  results  of  the  correlation  between  the  semiempirical 
design  methodology  and  the  experimental  data  indicate  several  shortcomings 
in  the  analysis  capabilities  as  well  as  lack  of  essential  test  data. 

The  application  of  existing  postbuckling  methodology  to  the  design 
of  advanced  composite  panels  has  resulted  in  unconservative  designs  due  to 
the  presence  of  additional  failure  modes,  such  as  delamination  in  the  skin, 
stiffener/web  separation  and  compression  failure  of  the  skin.  Such  failure 
modes  are  not  accounted  for  in  the  existing  methodology.  These  deficiencies 
in  the  analysis  and  design  methodology  have  to  be  corrected  to  realize  the 
full  weight  savings  potential  of  postbuckled  designs. 

Aircraft  panels  operating  in  the  postbuckling  range  are  usually 
curved.  These  panels  are  subjected  to  shear,  compression  (or  tension),  and 
a  combination  of  compression  (or  tension)  and  shear  loads.  Extensive  test 
data  are  available  for  flat  metal  and  composite  panels  subjected  to  shear  or 
compression  loading.  However,  the  test  data  for  curved  shear  or  compression 
panels  are  minimal  and  insufficient  to  develop  and  verify  improved  analysis 
methods.  These  gaps  in  the  available  test  data  need  to  be  filled. 
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Damage  tolerance  and  durability  requirements  for  the  design  of 
U.S.  military  aircraft  call  for  the  economic  life  of  an  aircraft  to  be  greater 
than  the  design  service  life  when  subjected  to  the  design  service  loads.  A 
fatigue  design  methodology  for  metal  or  composite  postbuckled  panels  is  not 
available  and  test  data  are  sketchy  or  incomplete.  Additional  experimental 
data  must  be  generated  to  identify  and  meet  the  fatigue  requirements  of  post- 
buckled  designs. 

1.2  PROGRAM  OBJECTIVES 

In  view  of  the  above  postbuckling  technology  needs,  a  combined 
analytical  and  experimental  program  was  undertaker,  to  develop  a  unified  de¬ 
sign  methodology,  design  validation  data  and  fatigue  life  data  for  composite 
and  metal  panels  operating  in  the  postbuckled  regime.  The  specific  objec¬ 
tives  of  the  program  were  to  develop  an  experimentally  validated  analysis 
capability  and  simple  to  use,  yet  accurate,  design  procedures  for  curved 
metal  and  composite  postbuckled  panels  loaded  in  compression  or  shear.  In¬ 
herent  in  this  program  objective  was  the  need  to  develop  techniques  to  pre¬ 
dict  the  initial  buckling  load,  ultimate  failure  load  and  failure  mode,  and 
fatigue  life  of  postbuckled  panels. 

The  results  of  the  design  methodology  development  program  are 
documented  in  this  report. 

1 .3  PROGRAM  SUMMARY 

The  program  plan  was  to  first  review  the  available  analysis  and 
design  techniques  for  metal  and  composite  panels  and  establish  a  design  meth¬ 
odology  for  curved  postbuckled  panels  loaded  in  compression  or  shear.  This 
methodology  was  then  used  in  designing  curved  panels  for  a  test  program  to 
generate  design  validation  and  fatigue  life  data. 

All  panels  were  cylindrically  curved  and  had  a  radius  of  45  inches. 
The  composite  panels  were  stiffened  with  hat  section  stringers  and  in  the  case 
of  the  composite  shear  panels  J-section  frames  were  used  as  circumferential 
stiffeners.  The  composite  compression  panels  were  not  circumferentially 
stiffened  in  the  test  section  since  they  simulated  the  region  between  two 
adjacent  bulkheads.  The  metal  panels  were  stiffened  with  Z-section  stringers 
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and  in  the  case  of  shear  panels  circumferential  Z-section  frames.  Fabrication 
of  the  composite  panels  was  accomplished  using  specially  designed  tooling  with 
careful  attention  paid  to  details  at  the  stringer/frame  intersection. 

The  test  plan  called  for  static  and  fatigue  tests  on  the  curved 
metal  and  composite  panels  loaded  either  in  compression  or  in  shear.  The 
static  test  data  were  used  to  verify  the  semiempirical  design  methodology, 
whereas  the  fatigue  test  data  were  utilized  to  determine  the  fatigue  failure 
modes  and  obtain  load  versus  life  data  to  formulate  fatigue  analyses  approaches. 

Comparison  of  the  static  test  data  with  the  semiempirical  design 
methodology  demonstrated  a  need  to  develop  a  more  rigorous  analysis  procedure 
for  postbuckled  panels.  In  addition,  the  fatigue  analysis  method  proposed 
requires  that  the  local  skin  displacements  be  known.  Therefore,  rigorous 
analysis  methods  using  the  principle  of  minimum  potential  energy  were  developed 
for  the  compression  and  the  shear  panels.  The  predictions  from  the  rigorous 
analysis  are  compared  with  the  test  data. 

Finally,  the  semiempirical  design  methodology  for  curved  composite 
and  metal  panels  subjected  to  compression  or  shear  loads  in  the  postbuckled 
regime  was  documented  in  a  Preliminary  Design  Guide. 

1.4  REPORT  OUTLINE 

The  program  was  performed  in  four  tasks.  Task  I  consisted  of 
selecting  analysis  methods  and  design  procedures  for  postbuckled  metal  and 
composite  panels.  These  methods  were  selected  from  the  available  technology 
on  postbuckled  structures  design.  The  selected  design  methodology  and  an 
assessment  of  the  technology  available  prior  to  the  start  of  this  program  are 
documented  in  Section  2.  Design  of  postbuckled  composite  and  metal  panels, 
panel  fabrication  and  testing  were  accomplished  in  Task  II  entitled  "Experi¬ 
mental  Test  Program."  A  description  of  Task  II  and  the  test  data  obtained  are 
documented  in  Section  3.  Task  III  consisted  of  comparing  the  test  results 
with  the  predictions  based  on  the  design  methodology  selected  in  Task  I,  and 
of  developing  the  more  rigorous  strain  energy  based  analyses  for  the  compres¬ 
sion  and  shear  panels.  Under  Task  III,  a  fatigue  analysis  approach  for  post¬ 
buckled  panels  was  also  developed  on  the  basis  of  the  test  data  generated  in 
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Task  II.  Development  of  the  compression  and  shear  panel  analyses  based  on  the 
principle  of  minimum  potential  energy  is  detailed  in  Section  4  along  with  the 
proposed  fatigue  analysis  approaches .  Correlation  of  test  data  with  results  of 
the  analysis  is  discussed  in  Section  5.  The  Preliminary  Design  Guide  developed 
in  Task  IV  is  published  separately  as  Volume  III  of  the  Final  Report.  The  pro¬ 
gram  conclusions  and  recommendations  for  future  work  are  summarized  in  Section 


SECTION  2 


SELECTION  OF  ANALYTICAL  METHODS 

2.1  BACKGROUND 

The  large  deflections  associated  with  postbuckled  structures 
make  the  elementary  theories  of  structural  analysis  inapplicable  to  determining 
the  detailed  stress  field  in  such  structures,  and  in  general,  closed-form  or 
analytical  solutions  cannot  be  obtained.  In  current  design  practice,  semi- 
empirical  analysis  methods  are  most  widely  used  in  sizing  stiffened  panels 
operating  in  the  postbuckled  range.  The  semiempirical  analysis  methods  are 
attractive  for  use  in  actual  design  situations  due  to  their  simplicity,  ease 
of  application  and  the  built-in  conservatism  in  the  static  analysis  results. 

The  main  drawback  of  these  methods,  besides  the  weight  penalty  associated  with 
the  conservatism  in  the  analysis,  is  that  they  do  not  provide  a  detailed 
stress  or  displacement  field  in  the  postbuckled  stiffened  panels.  Determina¬ 
tion  of  the  stress  field  is  essential  in  formulating  a  viable  fatigue  analysis 
approach.  A  few  numerical  solutions  have  been  attempted,  and  although  these 
techniques  provide  the  local  displacements  and  stresses,  they  are  cumbersome 
to  use,  and  too  expensive  in  terms  of  computer  costs  to  be  considered  viable 
design  tools. 

The  design  and  analysis  approach  adopted  in  this  program  was  to 
review  the  available  technology  base  on  postbuckled  metal  and  composite  panels 
and  then  select  a  semiempirical  methodology  for  modification  and  subsequent 
use  in  the  test  program.  The  purpose  of  the  modifications  was  to  extend  the 
applicability  of  the  analysis  techniques,  developed  for  metal  panels,  to  com¬ 
posite  panels.  The  semiempirical  methodology  was  used  to  design  the  program 
test  panels  which  in  turn  provided  verification  test  data.  Recognizing  the 
previously  mentioned  drawbacks  in  the  semiempirical  techniques,  development 
of  a  new  non-empirical  analysis  methodology  with  the  objective  of  predicting 
the  total  postbuckling  behavior  of  compression  or  shear  loaded  panels  was  also 
undertaken.  It  was  envisioned  that  the  results  of  this  rigorous  analysis 
would  find  immediate  application  in  formulating  approaches  to  fatigue  analysis 
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of  postbuckled  stiffened  panels. 


In  this  section,  an  assessment  of  the  current  technology  related 
to  design,  analysis  and  fatigue  life  prediction  of  postbuckled  structures  is 
presented  and  the  analytical  metl  ods  selected  for  use  in  the  program  are  de¬ 
tailed.  The  objectives  of  this  t«.  hnology  assessment  were:  (a)  to  enable 
selection  of  promising  analytical  methods  for  further  verification,  and  for 
test  panel  design;  and  (b)  to  compile  test  data  that  could  be  used  in  charac¬ 
terizing  the  fatigue  behavior  of  postbuckled  panels.  The  technology  assess¬ 
ment  was  also  used  to  select  an  approach  to  developing  a  non-empirical  analy¬ 
sis  methodology  for  postbuckled  compression  and  shear  panels. 

2.2  TECHNOLOGY  ASSESSMENT 

An  exhaustive  review  of  several  preliminary  design  studies  and 
test  programs  was  conducted  towards  selecting  design  methods  for  postbuckled 
metal  and  composite  panels  loaded  in  compression  or  shear.  The  currently 
available  design  and  analysis  methodologies  were  assessed  for  simplicity, 
accuracy  and  generality.  The  following  paragraphs  describe  in  a  summary 
form  the  state-of-the-art  methodology  for  design,  analysis,  and  durability 
validation  of  postbuckled  metal  and  composite  panels. 

2.2.1  Semi-Empirical  Static  Analysis  Methods 

The  semi-empirical  analysis  methods  evolved  from  an  extensive 
data  base  for  metal  panels  and  their  final  form  is  a  result  of  several  modifi¬ 
cations.  The  numerous  sources  of  postbuckled  panel  test  data  surveyed  and 
their  contributions  to  the  development  of  the  semi-empirical  analysis  tech¬ 
niques  are  summarized  in  Table  2.1.  A  detailed  discussion  of  the  significant 
contributions  is  presented  in  the  following  paragraphs. 

Shear  Panels  -  The  semi-empirical  analysis  and  design  method 
currently  used  for  stringer  stiffened  panels  loaded  in  shear  or  a  combination 
of  shear  and  compression  loading  has  evolved  over  the  years  from  the  "tension 
field  theory."  This  theory  was  originally  conceived  by  Wagner  (Reference  1) 
for  thin  flat  metallic  shear  webs.  Based  on  the  results  of  several  hundred 
tests,  some  empirical  constants  were  introduced  in  this  theory  to  broaden  its 
applicability  by  Kuhn  (Reference  2) .  The  theory  was  extended  to  the  analysis 
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of  curved  panels  as  well.  The  empirical  constants  were  chosen  such 
that  the  theory  yields  conservative  results  over  the  entire  range  of  its 
applicability. 


The  tests  conducted  by  Kuhn  (Reference  2)  at  NACA  generated 
elaborate  experimental  data  which  were  used  to  formulate  the  theory,  while 
tests  conducted  throughout  the  industry  verified  the  strength  predictions. 

The  flat  beams  tested  by  NACA  may  be  divided  into  three  groups:  small  but 
heavily  loaded  beams  (12  inches  deep) ;  medium-sized  beam  (25  or  40  inches 
deep),  which  formed  the  largest  group;  and  large  beams  (75  inches  deep). 

Small  but  heavily  loaded  beams  were  tested  using  a  simply 
supported  beam  loaded  at  the  beam  center.  The  lateral  movement  of  the  com¬ 
pression  chords  was  restricted  to  prevent  lateral  buckling.  The  medium  sized 
and  large  beams  were  tested  as  eccentrically  loaded  cantilever  beams.  Up¬ 
rights  or  stiffeners  on  these  panels  consisted  of  single  or  double  angles 
or  Z-sections.  A  number  of  variations  in  upright  and  web  geometry  were 
tested.  The  experimental  data  generated  from  these  tests  are  well  docu¬ 
mented  in  Reference  2  and  are  quite  useful  for  a  design  methodology  develop¬ 
ment  over  a  large  spectrum  of  panel  configurations. 

In  the  late  1960's,  the  results  of  14  full-scale  shear  beams 
were  reported  in  Reference  3.  The  beams  had  thin,  chem-milled  aluminum  webs 
with  lands  to  which  lipped  and  unlipped  Z-stiffeners  and  T-flanges  of  the 
same  material  were  riveted.  The  web  and  stiffener  thicknesses  ware  varied 
to  verify  the  tension  field  theory.  The  panels  were  tested  at  Room  Temperature 
Dry  (RTD)  conditions  through  the  use  of  an  eccentrically  loaded  cantilever 
beam.  The  load  was  incremented  in  small  steps.  Between  increments  the  load 
was  dropped  to  zero  to  observe  permanent  buckling  in  the  web.  Based  upon  the 
results  of  static  tests  some  modifications  to  the  tension  field  theory  were 
recommended . 

As  part  of  an  SST  technology  program  a  number  of  aluminum  and 
titanium  shear  panel  tests  were  conducted  as  discussed  in  Reference  4. 

Results  of  18  aluminum  and  titanium  shear  panel  tests  were  reported.  Panels 
were  tested  through  the  use  of  eccentrically  loaded  cantilever  beams.  The 
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test  panels  consisted  of  angle,  J-  and  Z-section  stiffeners  which  were  riveted 
to  the  panel  web.  The  aluminum  and  titanium  panels  had  identical  configura¬ 
tions  to  evaluate  the  usefulness  of  tension  field  theory  to  accommodate  dif¬ 
ferent  material  systems.  The  tests  conducted  on  aluminum  and  titanium  panels 
under  this  program  pointed  out  deficiencies  of  tension  field  theory  and  pro¬ 
vided  very  useful  data  for  future  analysis  development. 

Almost  all  of  the  information  available  on  curved  web  systems  in 
diagonal  tension  has  been  obtained  on  circular  cylinders  tested  in  torsion. 
Because  cylinders  are  more  expensive  to  manufacture  and  test  than  plane  web 
beams,  the  ^otal  number  of  cylinder  tests  is  rather  small.  The  tension  field 
theory  for  curved  panels  has  thus  seen  a  limited  verification.  A  total  of  12 
tests  were  reported  in  Reference  2.  In  these  tests  an  unconventional  arrange¬ 
ment  of  double  stringers  was  used  to  eliminate  bending  stresses  in  the  stringer. 
The  rings  were  made  relatively  large  to  preclude  ring  failure.  All  the  cyl¬ 
inders  tested  were  15  inches  in  radius,  a  dimension  which  is  not  quite  repre¬ 
sentative  of  real  aircraft  panel  configuration.  However,  variety  of  ring 
spacings  and  web  thicknesses  were  tested.  All  cylinders  were  tested  static¬ 
ally  under  RTD  conditions.  The  results  of  these  cylinder  tests  are  well 
documented  in  Reference  2  and  have  been  used  to  modify  the  tension  field 
theory  for  application  to  curved  web  systems. 

The  only  results  available  for  metal  panels  subjected  to  combined 
loading  are  reported  in  Reference  5.  Five  cylinders  were  tested  under  com¬ 
bined  compression  and  shear  loading.  All  cylinders  were  identical  in  con¬ 
struction  and  geometry.  The  cylinders  were  15  inches  in  radius  and  contained 
Z-section  stringers  and  rings  which  were  bolted  to  the  cylinder  web.  The 
ratio  of  torsional  and  compression  load  was  varied  for  each  cylinder.  All 
the  cylinders  were  tested  statically  under  RTD  conditions.  The  results  of 
these  tests  are  uniquely  suited  to  verification  of  analyses  for  metal  curved 
panels  under  combined  loading. 

Data  from  the  tests  described  in  the  preceding  paragraphs 
(References  2  through  5)  were  used  to  examine  the  validity  of  the  tension 
field  theory.  The  results  showed  that  the  internal  loads  predicted  by  the 
tension  field  theory  were  conservative  by  as  much  as  50  percent  for  alumi¬ 
num  beams  and  by  as  much  as  90  percent  for  the  titanium  beams.  In  some 
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cases,  failure  load  predictions  based  on  local  crippling  of  stiffeners  were 
found  to  be  unconservative.  The  primary  reason  for  this  unconservatism  in 
predicting  stiffener  failure  was  the  inability  of  tension  field  theory  to 
accurately  compute  stresses  in  the  eccentric  stiffeners,  especially  for  panel 
loads  several  times  the  initial  skin  buckling  load. 

The  introduction  of  advanced  composite  materials  as  viable  :an- 
didates  in  airframe  usage  has  stimulated  a  large  number  of  studies  in  recent 
years  to  determine  their  postbuckling  behavior.  One  of  the  earliest  demon¬ 
strations  of  the  postbuckling  strength  of  composite  materials  was  presented 
in  Reference  6.  Several  boron/epoxy  unstiffened  shear  panels  were  tested 
statically  and  under  fatigue  loading  at  load  magnitudes  several  times  their 
initial  buckling  load.  In  similar  studies  in  References  7,  8  and  10  the 
postbuckling  strength  of  graphite/epoxy  and  Kevlar/epoxy  materials  was  demon¬ 
strated.  The  test  data  for  composite  flat  shear  panels  were  obtained  in  these 
studies  through  the  use  of  picture  frame  test  setups.  This  test  setup  results 
in  the  introduction  of  severe  stress  concentrations  in  the  diagonal  corners, 
which  influence  the  panel  failure  load  as  well  as  the  mode  of  failure.  Thus, 
the  test  data  generated  in  the  above  studies  are  of  limited  use  for  the  pur¬ 
pose  of  failure  analysis  development  and  strength  verification.  However,  the 
test  data  are  useful  in  determining  the  panel  response  before  failure.  Some 
of  the  unstiffened  composite  shear  panel  data  were  obtained  from  tests  in  an 
eccentrically  loaded  cantilever  beam  setup.  Although  the  strain  data  gene¬ 
rated  in  these  tests  are  not  extensive,  the  data  are  more  reliable  and  suited 
to  failure  analysis  verification. 

Several  composite  stiffened  flat  panels  have  been  tested  under 
shear  loading  over  the  past  decade  through  the  use  of  various  test  methods. 

In  Reference  7,  several  panels  containing  bonded  hat  and  angle  section  stiff¬ 
eners  were  examined.  The  panel  web  thickness  was  varied  to  study  its  effect 
on  postbuckling  behavior.  The  panels  were  tested  in  a  picture  frame.  Al¬ 
though  severe  stress  concentrations  affected  most  of  the  panel  failure  modes, 
some  panels  failed  away  from  the  region  of  stress  concentrations.  An  exami¬ 
nation  of  these  results  indicates  that  variations  in  panel  thickness  are  most 
influential  in  determining  the  magnitude  of  postbuckling  deformations  which 
a  panel  can  sustain.  Panels  containing  thin  webs  can  sustain  loads  much 
higher  than  their  initial  buckling  load  without  failure,  as  compared  to 


panels  containing  thicker  webs.  For  example,  a  panel  containing  a  4-ply  web 
failed  at  a  load  of  10  times  its  initial  buckling  load,  whereas  a  similar 
panel  containing  an  8-ply  web  failed  at  a  load  5  times  its  initial  buck¬ 
ling  load.  The  failures  in  most  cases  were  due  to  stif fener/web  separation. 
The  stiffener  shape,  as  would  be  expected  in  practical  structural  designs, 
does  not  seem  to  affect  the  failure  load  or  mode  of  failure  in  any  signifi¬ 
cant  way. 

In  Reference  13,  the  influence  of  various  design  parameters  on 
the  behavior  of  composite  tension  field  panels  was  investigated.  The  stiff¬ 
ener  spacing  was  varied  to  change  the  onset  of  initial  buckling  load.  The 
stiffener  shape  was  varied  to  study  the  differences  between  the  behavior  of 
a  panel  containing  a  closed  section  ("hat")  and  an  open  section  ("I")  stiff¬ 
ener.  The  graphite-epoxy  stiffened  panels  in  this  study  were  designed  using 
the  metal  panel  tension  field  theory  with  modifications  that  account  for  the 
directional  dependence  of  the  composite  web  moduli.  The  panels  were  tested 
under  static  and  constant -amplitude  fatigue  loading.  The  test  results  were 
used  to  verify  the  design  methodology.  Two  failure  modes  somewhat  different 
from  metal  panel  failure  modes  were  discovered.  One  mode  of  failure,  stiff- 
ener/web  separation,  was  due  to  separation  of  stiffeners  from  the  panel  web. 
The  second  mode,  compression  failure,  was  due  to  deep  buckles  resulting  in 
large  compressive  stresses  in  the  web  corners.  The  test  fixture  used  in 
this  study  was  an  eccentrically  loaded  cantilever  beam  which  in  Reference  2 
has  beeu  demonstrated  to  apply  a  uniform  shear  to  the  panel.  The  stiffener 
shape  (Hat  or  I-section)  did  not  seem  to  affect  pane]  posfbuckling  behavior. 
Additional  static  tests  on  similar  panels  were  performed  in  Reference  14. 
Heavily  instrumented  and  carefully  replicated  tests  established  the  static 
behavior  of  the  panels.  In  a  similar  manner,  the  fatigue  behavior  of  iden¬ 
tical  panels  was  studied  for  panels  subjected  to  fully  reversed  fatigue 
loading.  The  panels  consisted  of  cocured  hat  stiffeners.  The  panel  failure 
was  due  to  stiffener/web  separation.  These  test  data  proved  valuable  for 
analysis  verification  and  modifications. 

In  Reference  15,  J-stiffened  composite  shear  webs  were  examined. 
The  test  specimen  contained  three  integrally  cocured  J-stiffeners.  The  panels 
were  tested  statically  as  well  as  under  constant-amplitude  fatigue  loading 
(with  no  load  reversal).  The  test  setup  used  for  testing  these  panels  was 
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a  "picture  frame."  The  statically  tested  panels  failed  prematurely  due  to 
separation  of  the  stiffeners  from  the  panel  skin  in  the  diagonal  corners  of 
the  test  fixture.  The  results  of  these  tests  may  be  useful  in  predicting 
the  lower  bound  of  panel  strength. 

In  Reference  16,  the  test  specimens  contained  two  hat  stiffeners 
(simulating  a  fuselage  panel),  and  two  blade  stiffeners  (simulating  the  fuse¬ 
lage  frames  or  bulkheads) .  The  panels  were  tested  through  the  use  of  a  mod¬ 
ified  picture  frame  in  which  the  load  is  applied  along  one  edge,  as  opposed 
to  along  the  panel  diagonal  in  a  conventional  picture  frame  test  setup.  The 
panels  tested  statically  under  this  program  failed  due  to  separation  of  the 
blade  stiffeners  from  the  skin.  The  hat  stiffeners  were  unaffected.  From 
phenomenological  considerations  it  would  seem  logical  for  stiffener/web  sep¬ 
aration  to  occur  between  hat  stiffeners  and  the  panel  skin  where  maximum 
peel  and  interlaminar  shear  stresses  are  introduced  due  to  the  buckles.  The 
relative  magnitude  of  peel  and  interlaminar  stresses  should  be  small  near 
the  blade  stiffeners.  Thus,  the  blade  stiffener/skin  separation  mode  of 
failure  was  not  anticipated. 

In  a  recently  completed  Advanced  Composite  Center  Fuselage  Pro¬ 
gram  (References  17  and  18),  a  few  graphite /epoxy  stiffened  and  unstiffened 
shear  panels  were  tested.  The  tests  were  conducted  on  panels  subjected  to 
different  environmental  conditions  as  well  as  on  panels  containing  impact 
damage.  None  of  the  test  conditions  was  replicated.  Because  the  panels  were 
tested  as  part  of  a  design  verification  program,  no  attempt  was  made  to  ob¬ 
tain  extensive  strain  distributions  or  to  replicate  test  conditions.  Al¬ 
though  the  test  results  have  limited  value  in  analysis  verification,  they 
do  provide  valuable  design  data.  Evaluation  of  a  variety  of  design  concepts 
and  stiffener/web  interface  improvements  to  increase  the  strength  of  post- 
buckled  composite  panels  was  conducted  in  References  19  through  25. 

In  Reference  19,  a  Kevlar  panel  with  embedded  graphite  plies 
in  the  stiffener  flanges  was  examined  for  postbuckling  strength.  However, 
uncertainty  in  the  panel  quality  casts  some  doubt  on  the  reliability  of  the 
test  results. 

The  shear  panels  tested  in  Reference  20  incorporated  a  design 
improvement  at  the  stiffener/web  interface.  The  improvement  consisted  of 
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applying  a  single  layer  of  FM-300  film  adhesive  between  the  interface  be¬ 
fore  the  panel  was  cured.  Another  panel  was  fabricated  with  an  extra  layer 
of  3501-6  resin  at  the  interface.  Both  panels  when  tested  showed  signifi¬ 
cant  improvement  in  the  stiffener /web  interface  strengths. 

The  postbuckling  behavior  of  flat  stiffened  multibay  composite 
panels  was  examined  in  Reference  21.  The  panels  were  subjected  to  a  combi¬ 
nation  of  inplane  shear,  compression,  and  bending  loads.  Because  test  re¬ 
sults  for  panels  under  combined  load  are  rare  in  the  literature,  the  results 
of  this  study  are  valuable  for  preliminary  analysis  verification.  The 
panels  tested  under  this  program  consisted  of  a  10-ply  graphite/epoxy  web 
with  cocured  hat  section  stringers  and  bead  section  frames.  The  panels  were 
designed  to  buckle  at  30  percent  of  the  design  limit  load.  Extensive  in¬ 
strumentation  was  provided  to  measure  the  stress  distribution  as  well  as  the 
postbuckling  behavior.  A  significant  loss  in  the  panel  initial  buckling  load 
was  reported  due  to  fatigue  loading. 

The  postbuckling  behavior  of  graphite/epoxy  polyimide  panels 
was  examined  under  a  Navy  contract  (Reference  22)  .  The  panels  in  this  pro¬ 
gram  were  identical  to  the  panels  tested  earlier  in  Reference  14  except  for 
the  materials.  The  results  of  this  program  showed  that  the  strength  of  poly¬ 
imide  materials  under  out-of-plane  loads  was  rather  poor.  This  poor  strength 
was  manifested  as  skin/stiffener  separation  which  was  a  severe  problem  in 
these  panels. 

As  part  of  the  Wing/Fuselage  program  (Reference  25) ,  several 
composite  stiffened  flat  panels  were  designed,  fabricated  and  tested.  These 
panels  consisted  of  two  cocured  hat  stiffeners.  The  overall  panel  dimensions 
were  22  by  26  inches.  A  total  of  12  panels  were  tested  to  determine  the 
postbuckling  behavior  under  combined  compression  and  shear  loading.  Repli¬ 
cated  tests  were  used  to  examine  the  static  as  well  as  fatigue  behavior  under 
RTD  conditions  as  well  as  under  Elevated  Temperature  Wet  (ETW)  conditions. 

A  flight  by  flight  spectrum  loading  typical  of  a  Mach  2  fighter  aircraft  was 
used  to  study  durability  of  these  panels  for  two  lifetimes  of  fatigue  load¬ 
ing.  The  results  of  this  study  are  potentially  useful  in  developing  a  design 
methodology  i or  postbuckled  panels  subjected  to  combined  loading. 
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Another  design  improvement  for  the  stiffener/web  interface  was 
examined  in  Reference  24.  The  results  of  this  study  show  that  tailoring  the 
skin/stiffener  interface  results  in  improved  fatigue  life.  The  tailoring 
consisted  of  tapering  the  flanges  of  the  stiffeners  at  the  interface  with 
the  skin.  This  tapering  was  accomplished  by  selectively  dropping  off  the 
flange  plies.  The  influence  of  environment  and  combined  loads  on  flat  com¬ 
posite  stiffened  panels  is  examined  in  Reference  27.  This  study  is  currently 
in  progress  and  will  provide  much  needed  data  for  analysis  development  and 
verification. 

The  analysis  methodology  for  curved  composite  panels  under  shear 
loading  is  still  in  its  infancy.  This  is  primarily  due  to  a  lack  of  reliable, 
detailed  test  data.  The  postbuckling  strength  of  curved  composite  shear 
panels  was  demonstrated  in  Reference  7  where  a  graphite/epoxy  and  a  boron/ 
epoxy  panel  were  tested  through  the  use  of  a  cylinder  torque  test.  These 
data  in  conjunction  with  data  from  Reference  27  and  those  being  generated 
in  Reference  32  will  prove  valuable  in  verifying  the  applicability  of  the 
modified  tension  field  theory  to  curved  composite  shear  panels. 

Compression  Panels  -  The  semiempirical  postbuckling  analysis  of 
flat  and  curved  stiffened  panels  loaded  in  compression  is  generally  done  in 
steps  because  it  involves  several  complexities  which  are  difficult  to  account 
for  simultaneously.  The  method  normally  used  to  analyze  metal  panels  is  in 
four  parts: 

1.  Determine  the  panel  initial  buckling  load. 

2.  Determine  the  compressive  strength  of  the 
stiffener  alone. 

3.  Determine  the  effective  width  of  the  skin  for 
a  load  equal  to  the  compressive  strength  of 
the  stiffener  alone. 

4.  Determine  the  total  load  carried  by  the  panel 
by  taking  into  account  the  load  on  the  stiff¬ 
ener  plus  the  effective  width  of  the  skin, 
plus  the  critical  buckling  load  of  the  skin. 
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The  above  method  requires  calculation  of  the  panel  buckling  load, 
the  behavior  of  the  skin  after  buckling  and  the  failure  strength  of  a  column 
of  any  arbitrary  shape.  For  design  purposes,  the  above  process  is  generally 
repeated  several  times  to  obtain  positive  margins  on  all  structural  elements. 

Analytical  and  semiempirical  methods  for  predicting  the  initial 
buckling  load  of  the  skin  are  well  developed  for  metal  and  composite  panels 
and  are  extensively  documented  in  References  34  through  39,  for  example. 

These  methods  vary  widely  in  rigor  and  accuracy.  The  analysis  methods  de¬ 
veloped  for  metals  up  to  the  early  1950's  and  documented  in  References  34 
and  35,  are  semiempirical  in  nature  and  are  based  on  the  results  of  an  ex¬ 
tensive  test  data  base. 

Several  more  rigorous  analytical  methods  were  made  possible  by 
the  evolution  of  high-speed  digital  computers.  Some  of  the  most  recent  and 
advanced  analysis  methods  for  linear  bifurcation  buckling  analysis  are  de¬ 
scribed  in  References  40,  41  and  42  (Computer  codes  SX8,  BUCLASP2,  and 
VIPASA) .  Computer  code  SX8  is  based  upon  the  Rayleigh-Ritz  energy  principle 
for  analysis  of  flat  composite  and  metal  stiffened  panels.  The  stiffeners 
are  assumed  to  be  axial  members  and  their  effect  on  panel  behavior  is  in¬ 
cluded  by  taking  into  consideration  the  bending  and  torsional  stiffnesses 
of  the  discrete  stiffeners  in  the  energy  expression  for  the  panel.  In  doing 
so,  the  effects  of  stiffener  shape  are  neglected.  The  main  advantage  of 
this  method  is  that  arbitrary  boundary  conditions  along  the  panel  edges  can 
be  analyzed. 

The  computer  codes  BUCLASP-2  and  VIPASA  (References  41  and  42) 
are  quite  similar  to  each  other  except  in  their  ability  to  solve  cases  with 
different  loading  and  boundary  conditions.  These  methods  are  based  on  solv¬ 
ing  exact  force-displacement  relations  for  a  plate-strip  element  with  the 
assumption  of  simple  support  boundary  conditions  along  the  edges  normal  to 
the  longitudinal  direction  and  arbitrary  boundary  conditions  along  the  long¬ 
itudinal  edges.  An  assembly  method,  similar  to  the  one  used  in  finite  ele¬ 
ment  analysis,  is  used  to  generate  any  desired  panel  configuration.  The  ad¬ 
vantage  of  this  method  of  analysis  is  that  both  general  and  local  instability 
modes  can  be  simultaneously  predicted.  Additionally,  any  combination  of 
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local  and  general  Instability  modes,  which  results  in  lower  buckling  loads, 
can  also  be  predicted.  This  method  of  analysis  was  used  very  successfully 
to  correlate  experimental  data  for  hat-stiffened  and  J-stiffened  graphite- 
epoxy  compression  panels  in  References  43  and  44. 

The  analytical  methods  for  determining  the  onset  of  buckling  for 
curved-stiffened  panels  are  limited  in  number  and  have  developed  along  lines 
similar  to  the  corresponding  methods  for  flat  panels.  The  methods  described 
in  References  41,  42  and  45  (Computer  codes  BUCLASP-2,  VIPASA,  and  SS8)  have 
been  used  more  recently  throughout  the  aircraft  industry  for  metal  and  com¬ 
posite  panels.  These  methods  can  also  be  used  for  flat-stiffened  panels. 

In  order  to  predict  stiffened  panel  behavior  after  initial  buck¬ 
ling  and  the  ultimate  compressive  strength  of  stiffeners,  semiempirical  tech¬ 
niques  have  to  be  utilized.  The  semiempirical  methodology  for  predicting 
local  buckling  and  crippling  strength  of  metal  stiffeners  has  been  derived 
from  a  large  data  base  and  is  documented  in  References  34  and  39  for  various 
stiffener  configurations.  Test  data  on  the  crippling  strength  of  composite 
stiffeners  are  sparse  and  as  a  result  no  definitive  analysis  techniques  exist. 
One  approach  suggested  in  Reference  46  is  based  on  tests  performed  on  several 
composite  stiffener  elements.  These  tests  included  plates  with  both  sides 
simply  supported  as  well  as  plates  with  one  edge  free  and  the  other  edge 
simply  supported.  In  Reference  46,  the  following  empirical  equation  was 
suggested  for  use  in  predicting  the  crippling  strength  of  composite  plates: 
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In  the  above  equation,  Fcc  is  the  crippling  stress,  F^r  is  the  theoretical 
initial  buckling  stress,  and  F£u  is  the  plate  compressive  strength.  The 
values  of  a  and  3  are  mate.rlali-fldependent  constants.  The  values  of  these 
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constants  for  AS/3501-6  and  T300/5208  graphite/epoxy  systems  were  obtained 
in  Reference  46  by  fitting  a  curve  to  the  test  data.  Validity  of  the  em¬ 
pirical  equation  in  predicting  the  crippling  strength  of  built-up  sections 
was  demonstrated  in  Reference  46  by  tests  on  square  composite  tubes. 
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Generally,  metal  compression  panels  sized  using  the  above  ap¬ 
proach  are  conservatively  designed.  In  recent  studies  (References  29  and 
30),  composite  compression  panels  were  designed  using  the  above  approach. 

The  panels  consisted  of  cocured  hat-stiffeners  and  were  designed  to  buckle 
at  loads  significantly  below  the  panel  ultimate  strength.  The  methods  used 
for  determining  initial  buckling  load  and  stiffener  crippling  strength  were 
those  mentioned  above.  The  replicated  tests  indicated  that  the  analysis 
methods  were  sufficiently  accurate,  although  on  the  conservative  side.  How¬ 
ever,  additional  verification  of  the  semiempirical  methodology  for  composite 
panels  was  necessary  and,  therefore,  carried  out  in  this  program.  Test  data 
obtained  in  References  26,  29  through  31,  47  and  48  were  used  for  verifica¬ 
tion.  This  process  of  correlating  the  data  base  with  the  semiempirical  pre¬ 
dictions,  resulted  in  some  modifications  to  the  crippling  strength  predic¬ 
tion  equation  given  above.  These  results  and  modifications  are  discussed 
in  Section  2.3. 

2.2.2  Non-empirical  Static  Analysis  Methods 

Analysis  of  stringer  stiffened  panels  loaded  beyond  skin  buck¬ 
ling  requires  solution  of  nonlinear  equations  and  no  closed-form  expressions 
describing  the  response  can  be  obtained.  Thus,  numerical  methods  must  be 
resorted  to.  The  numerical  methods  of  solution  are  generally  iterative  in 
nature,  and  their  efficiency  and  utility  are  limited  by  the  number  of  unknown 
variables  in  the  solution  process.  Two  of  the  most  commonly  used  nonempiri- 
c.al  solution  methods  for  postbuckling  analysis  are  Finite-Difference/Finite 
Element  Methods  (FD/FEM)  and  Rayleigh-Ritz  type  methods.  Several  applica¬ 
tions  of  these  techniques  to  postbuckled  structures  are  discussed  in  the 
following  paragraphs. 

Finite  Difference/Finite  Element  Methods  -  The  main  advantage 
of  using  FD/FEM  is  that  several  large  computer  codes  exist  (e.g.,  NASTRAN, 
STAGS,  ANSYS)  and  are  easy  to  access  in  the  aerospace  industry.  However, 
a  few  basic  problems  in  using  these  techniques  make  them  undesirable  for 
design  purposes.  The  problems  are  related  to  questions  of  accuracy,  effi¬ 
ciency,  and  convergence. 
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The  accuracy  of  the  finite  element  solution  method  depends  upon 
the  type  and  number  of  finite  elements  used  to  model  structural  behavior. 
Since  the  structural  behavior  is  not  known  in  advance,  a  general  practice 
is  to  start  with  a  reasonably  fine  mesh  size  and  refine  it  subsequently  to 
establish  the  accuracy  of  the  solution.  Thus,  the  solution  process  can  be 
expensive  if  several  such  iterations  are  needed.  The  analytical  solutions 
of  problems  involving  large  deflections  are  iterative  in  nature  by  themselves 
and  this  added  iteration  can  make  the  solution  economically  unfeasible. 

The  efficiency  of  the  finite  element  solution  depends  upon  the 
nature  of  the  problem  and  also  upon  the  mesh  size  used  tc  model  structural 
response.  For  example,  in  order  to  model  the  stiffener /web  separation, 
three-dimensional  finite  elements  must  be  used  in  the  interface  area.  Since 
the  thickness  of  the  interface  is  quite  small,  a  large  number  of  elements 
will  be  required  in  this  region  to  avoid  numerical  difficulties  as  well  as 
to  model  the  behavior  accurately.  Recent  experience  at  Northrop  and  results 
in  *"he  literature  indicate  that  the  interface  stresses  are  sensitive  to 
element  size  in  structures  involving  small  displacements.  It  can  be  antic¬ 
ipated  that  similar  or  even  worse  difficulties  will  be  encountered  xn  panels 
subjected  to  postbuckling  deformations.  Again  due  to  the  iterative  nature 
of  the  solution  process,  a  large  number  of  elements  can  make  the  solution 
economically  unfeasible. 

The  governing  nonlinear  equations  for  postbuckled  structures 
are  solved  incrementally.  A  general  practice  is  to  increase  the  applied 
load  or  displacement  in  small  increments  with  the  size  of  these  increments 
determining  the  progress  of  the  solution  towards  the  maximum  applied  load. 
Since  the  increment  size  is  not  known  in  advance,  several  trials  are 
generally  needed  to  obtain  a  convergent  solution.  Convergence  difficulties 
are  also  encountered  due  to  the  size  of  the  mesh  used  in  making  a  geometric 
representation  of  the  problem.  Several  such  convergence  difficulties  have 
been  reported  in  References  49  through  57  and  are  discussed  below.  Another 
simplifying  assumption  made  in  the  application  of  FD/FEM  to  postbuckled 
structures  is  to  model  the  skin  separately  without  regard  to  the  interaction 
effects  of  the  s>.in  and  the  stringers.  Similar  drawbacks  also  apply  to 
finite  difference  solution  methods. 
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Due  to  the  difficulties  encountered  in  the  use  of  FD/FEM  methods, 
most  of  the  studies  reported  in  the  literature  on  postbuckled  structures  have 
I  been  conducted  making  several  simplifying  assumptions.  One  of  the  most  com- 

»  mon  assumptions  made  is  to  model  the  skin  of  the  skin-stringer  panel  separate 

ly,  without  regard  to  the  interacting  effect  of  skin  and  stringers. 

Sharifi  (Reference  49)  modeled  the  behavior  of  an  unstiffened 
shear  panel  subjected  to  postbuckling  loads,  using  finite  elements,  and 
showed  that  the  postbuckling  deflections  could  be  predicted  quite  accurate¬ 
ly.  In  a  similar  study  by  Bhatia  (Reference  9)  the  STAGS  program  was  used 
to  predict  the  postbuckling  behavior  of  composite  shear  webs.  In  these 
studies  no  attempt  to  predict  failure  was  made.  Turney  and  Wittrick  (Refer¬ 
ence  50)  and  Rushton  (Reference  51)  used  a  finite  difference  Iterative 
method  known  as  "dynamic  relaxation"  for  the  postbuckling  analysis  of  square 
plates  subjected  to  uniaxial  compression  and  shear.  Rectangular  finite 
elements  with  bi-cubic  Hermitian  interpolation  functions  vere  used  in 
Reference  52  to  study  the  postbuckling  behavior  of  uniaxially  compressed 
sandwich  panels. 

Postbuckling  behavior  of  stiffened  shear  webs  was  also  studied 
by  Stein  and  Starnes  (Reference  53)  through  the  use  of  the  STAGS  computer 
code.  They  conducted  parametric  comparisons  on  the  efficiencies  of  metal 
and  composite  shear  webs  loaded  up  to  about  twice  the  buckling  load.  Sev¬ 
eral  convergence  difficulties  were  pointed  out.  This  study  demonstrated 
the  usefulness  of  STAGS  in  performing  postbuckling  stress  analysis,  but 
failed  to  establish  the  accuracy  of  the  solution  process. 

Vestergren  and  Knutsson  (Reference  54)  also  used  STAGS  to  study 
the  postbuckling  behavior  of  unstiffened  compression  and  shear  panels.  The 
initial  buckling  loads  were  predicted  quite  accurately  for  compression  as 
well  as  shear  panels  using  the  bifurcation  analysis.  However,  the  failure 
load  predicted  for  compression  panels  was  twice  the  experimentally  obtained 
value.  No  data  were  presented  for  the  failure  load  predictions  of  shear 
panels.  Again,  several  difficulties  in  the  use  of  computer  code  STAGS  for 
analysis  purposes  were  acknowledged. 
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The  change  in  stiffness  of  a  plate  loaded  under  combined  load¬ 
ing  was  discussed  by  Rothwell  and  Allahyari  (Reference  55).  The  finite 
element  analysis  used  in  this  reference  provided  guidelines  for  minimizing 
the  loss  of  the  beam  flexural  stiffness  as  a  result  of  web  buckling. 

The  postbuckling  behavior  of  a  composite  panel  shear  web  (ex¬ 
cluding  the  stiffeners)  was  modeled  in  Reference  13  using  the  computer  code 
NASTRAN.  In  this  study,  the  web  was  assumed  to  be  simply  supported  on  all 
four  edges.  The  analytically  obtained  displacements  and  strains  in  the 
postbuckling  regime  were  shown  to  compare  favorably  with  experimental 
values.  Web  rupture  due  to  compressive  stresses  resulting  from  deep  buckles 
was  shown  to  be  predictable.  The  analytical  results  in  this  study  indi¬ 
cated  a  concentration  of  out-of -plane  constraint  forces  in  the  diagonal 
corners  where  failure  due  to  stiffener/web  separation  was  observed.  Since 
the  stiffener/web  interface  was  not  modeled  in  the  NASTRAN  analysis,  no 
accurate  prediction  of  failure  due  to  stiffener/web  separation  was  made. 
Further  attempts  to  enlarge  the  model  to  include  the  total  panel  behavior 
had  to  be  aborted  due  to  convergence  difficulties.  A  NASTRAN  finite  element 
analysis  of  postbuckled  shear  panels  was  also  attempted  in  Reference  16. 

The  initial  buckling  load  was  shown  to  agree  quite  well  with  the  experi¬ 
mental  data.  However,  the  analysis  attempt  was  again  aborted  above  150 
percent  of  the  initial  buckling  load  due  to  convergence  difficulties. 

In  Reference  29,  the  postbuckling  behavior  of  a  composite  stiff¬ 
ened  compression  panel  was  modeled  through  the  use  of  the  large  deflection 
theory  of  NASTRAN.  A  convergence  difficulty  resulting  from  using  a  rela¬ 
tively  coarse  mesh  size  was  encountered  after  the  load  exceeded  twice  the 
initial  buckling  load,  and  the  solution  attempt  was  aborced.  However,  a 
fairly  good  correlation  with  experimental  data  was  observed  for  the  results 
obtained. 

Several  recent  attempts  have  been  made  at  improving  the  effi¬ 
ciency  and  reducing  the  convergence  difficulties  of  finite  difference/finite 
element  methods  (References  56  and  57),  but  their  implementation  as  a  de¬ 
sign  tool  in  the  near  future  is  unlikely. 

Rayleigh-Ritz  Type  Methods  -  These  methods  are  widely  used  to 
model  the  behavior  of  complex  structures,  since  they  are  conceptually  simple 
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to  use.  The  displacement  field  is  approximated  by  assumed  functions  with 
unknown  coefficients.  The  number  of  equations  to  be  solved  is  reduced  sig¬ 
nificantly  when  compared  to  the  finite  element  methods,  thus,  reducing  the 
computation  time  and  cost.  Furthermore,  once  the  problem  is  formulated, 
parametric  studies  can  be  conducted  with  little  additional  cost.  However, 
there  are  some  difficulties  in  using  this  approach  also.  The  main  diffi¬ 
culty  arises  from  selecting  the  deformed  structural  shape.  This  is  usually 
resolved  by  selecting  a  shape  which  is  a  combination  of  several  possible 
shapes.  The  experimentally  observed  behavior  greatly  enhances  this  selec¬ 
tion  process.  The  difficulty,  which  is  common  to  all  numerical  methods 
(finite  element  methods  as  well),  is  that  the  computed  deflections  may  be 
quite  accurate,  but  the  computed  strains  tend  to  be  in  error.  This  diffi¬ 
culty  can  be  rectified  by  increasing  the  number  of  terms  in  the  assumed 
displacement  function.  This  increases  the  computation  time  required,  but 
the  relative  magnitude  of  the  increase  is  quite  small  for  the  Rayleigh- 
Ritz  type  of  analysis  as  compared  to  the  finite  element  analyses.  However, 
since  the  advantages  overshadow  the  disadvantages  in  a  Rayleigh-Ritz 
solution,  several  such  solution  methods  have  been  attempted  over  the  years. 
As  discussed  below,  these  studies  have  addressed  different  aspects  of  the 
postbuckling  problem. 

The  Rayleigh-Ritz  technique  was  used  to  analyze  an  incomplc.  e 
tension  field  stiffened  beam  by  Denke  (Reference  58).  The  wave  form  of  the 
buckled  surface  was  approximated  by  a  function  that  contained  the  wave 
length,  wave  angle,  and  wave  depth  as  parameters.  Four  additional  para¬ 
meters  -  namely,  stiffener  compressive  strain,  chord  compressive  strain, 
the  chord  bending  deflection  and  the  panel  shearing  strain  -  were  introduced 
to  account  for  the  effect  of  inplane  membrane  forces.  The  resultant 
governing  equations  were  solved  to  predict  the  principal  midplane  stresses, 
maximum  web  bending  stresses,  stiffener  and  chord  compressive  stress.  A 
comparison  with  limited  experimental  data  showed  reasonably  good  correla¬ 
tion.  This  analysis  was  limited  in  scope,  as  it  failed  to  include  rotation 
and  out-of-plane  bending  in  the  chords  and  stiffeners.  In  addition,  the 
limited  number  of  unknown  terms  (seven)  restricted  the  accuracy  of  the 
solution  process. 


24 


I 


Levy,  et  al  (References  59  and  60)  used  the  Von  Karman  large 
deflection  plate  equations  to  study  the  postbuckling  behavior  of  unstiff¬ 
ened  metal  plates  leaded  in  shear.  The  Von  Kansan  equations  were  solved 
approximately  by  assuming  a  truncated  Fourier  double  sine  series.  The 
maximum  number  of  terns  in  the  solution  was  relatively  small,  which  limited 
the  accuracy  of  the  results  obtained. 

The  analyses  in  the  above  studies  (References  58  through  60) 
were  limited  to  a  few  terms  in  the  assumed  functions  primarily  due  to  the 
absence  of  high  speed  computers.  Modem  computers  have  increased  the  feas¬ 
ibility  of  introducing  a  significantly  larger  number  of  terms  in  the  analy¬ 
sis  at  a  reasonable  cost. 

Several  later  studies  used  considerably  larger  numbers  of  terms 
to  obtain  more  accurate  solutions.  Mayers  and  Budiansky  (Reference  61)  used 
a  combination  of  algebraic  and  trigonometric  functions  to  represent  the  in¬ 
plane  and  out-of-plane  deflections  of  a  flat  plate  loaded  in  compression 
beyond  initial  buckling.  The  analytical  results  were  shown  to  agree  reasonably 
well  with  the  test  results.  Since  material  plasticity  was  not  included  in 
these  analyses,  a  small  difference  in  the  experimental  and  analytical  results 
was  to  be  expected  which  did,  in  fact,  occur. 

Chia  and  Prabhakara  (Reference  62)  presented  an  analysis  based 
on  the  Von  Karman  type  of  large  deflection  equations.  These  equations  were 
solved  by  expressing  the  force  function  and  transverse  deflection  as  a  double 
Fourier  series  in  terms  of  approximate  beam  Eigen  functions  for  unsymmetric- 
ally  laminated  rectangular  plates.  These  plates  were  subjected  to  uniaxial 
and  biaxial  compression.  Both  simply  supported  and  clamped-boundary  condi¬ 
tions  were  considered.  Harris  (References  63  and  64)  presented  approxi¬ 
mate  analytical  expressions  for  the  inplane  stiffness  immediately  after 
buckling  for  rectangular  composite  plates  subjected  to  biaxial  compression. 

Chan  (Reference  65)  presented  a  slightly  different  form  of 
Rayleigh-Ritz  analysis  to  obtain  the  postbuckling  behavior  of  compression 
loaded  composite  flat  plates.-.  The  solution  was  carried  out  with  only  the 
transverse  displacement  mode  assumed.  The  inplane  displacements  were  ob¬ 
tained  exactly  from  the  two  membrane  displacement  equations.  Although  this 
method  of  solution  reduces  the  number  of  simultaneous  nonlinear  equations 
to  be  solved,  it  introduces  a  few  additional  restrictions.  The  inplane 
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displacements  become  dependent  on  the  assumed  shape  of  transverse  deflec¬ 
tion  and  the  resulting  displacements  are  not  in  general  capable  of  accommo¬ 
dating  the  imposed  inplane  boundary  conditions.  In  addition,  this  analysis 
method  cannot  be  extended  to  the  analysis  of  stringer  stiffened  panels 
under  general  loading  conditions.  Further,  the  advantage  offered  by  this 
method  for  postbuckling  analysis  is  not  as  significant  because  in  conven¬ 
tional  applications  of  the  Rayleigh-Ritz  method,  most  of  the  resulting 
equations  involving  inplane  displacements  are  linear  in  terms  of  transverse 
displacement  coefficients  and,  thus,  can  easily  be  solved  in  terms  of  these 
transverse  deflection  coefficients. 

The  effect  of  inplane  flexural  and  axial  rigidity  of  stiffeners 
on  the  postbuckling  behavior  of  square  metal  panels  was  examined  in  Refer¬ 
ence  66.  The  panels  were  loaded  in  shear  and  compression.  The  Von-Karman 
equations  were  solved  by  assuming  the  shape  of  the  normal  displacement  and 
the  stress  function.  The  analysis  failed  to  include  the  torsion  and  out- 
of-plane  flexibility  of  the  stiffeners. 

In  an  attempt  to  improve  the  empirical  analysis  of  Kuhn,  a 
rather  rigorous  analysis  of  flat  tension  field  beams  was  formulated  during 
the  development  of  SST  technology  (Reference  67) .  The  problem  was  formu¬ 
lated  for  a  cantilever  beam  consisting  of  internal  stiffeners.  The  beam 
displacement  functions  were  selected  to  satisfy  the  boundary  conditions 
imposed  during  the  tests.  This  resulted  in  the  inclusion  of  the  shear 
and  bending  deformations  of  the  beam;  inplane  bending,  out-of-plane  bending, 
axial  deformation,  rotation,  and  warping  of  internal  stiffeners;  lnplane 
bending,  axial  deformation,  rotation,  and  warping  of  the  chords;  and  axial 
deformation,  rotation,  and  warping  of  the  two  edge  vertical  stiffeners. 

The  analysis  was  the  first  attempt  to  duplicate  the  exact  mechanism  of 
load  introduction  for  the  tension  field  shear  beam  specimen.  The  analysis 
was  formulated  and  some  results  were  presented  for  a  single-bay  panel.  The 
cancellation  of  the  SST  program  and  some  convergence  difficulties  in  the 
solution  process  halted  further  development  and  verification  of  the  analy¬ 
sis  for  multibay  panels. 

Khot  (References  68  and  69)  demonstrated  the  usefulness  of  this 
approach  by  studying  the  postbuckling  behavior  and  imperfection  sensitivity 
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of  composite  cylindrical  shells  loaded  under  axial  compression.  Dickson,  et 
al,  References  70  through  72,  formulated  the  problem  of  composite  stiffened 
panels  in  the  postbuckling  range  using  the  Rayleigh-Ritz  approach.  In  Refer¬ 
ence  72,  the  Rayleigh-Ritz  solution  procedure  haa  been  used  in  conjunction 
with  an  optimization  routine  to  design  a  curved  composite  stiffened  panel. 

The  analysis  has  also  been  extended  to  predicting  the  local  stress  state  at 
the  skin/stif fener  interface  in  stiffened  composite  panels  (Reference  73) . 
Experimental  evaluation  of  this  predictive  methodology,  however,  has  not 
been  carried  out. 

In  References  74  and  75  the  analysis  suggested  in  Reference  58 
for  metal  panels  was  modified  for  use  with  composite  panels  loaded  in  shear. 
The  results  of  this  analysis  were  compared  with  several  existing  composite 
panel  tests.  A  fairly  good  correlation  between  measured  and  predicted  maxi¬ 
mum  out-of-plane  deflection-  and  ir,;dane  strains  in  the  panel  center  was 
demonstrated.  This  analysis,  once  gain,  demonstrated  the  usefulness  of 
the  Rayleigh-Ri  tz  type  of  solutions  to  predict  the  postbuckling  behavior 
of  composite  panels.  Because  of  the  limited  number  of  terns  used  in  the 
assumed  functions,  the  predicted  and  measured  values  were  not  in  as  good 
agreement  near  the  panel  edges  as  at  the  center.  This  shortcoming  can  be 
improved  by  taking  additional  terms  in  the  solution.  Rapid  convergence 
coupled  with  nominal  computer  run  times  makes  the  approach  attractive  for 
design  purposes.  The  increasing  popularity  of  the  Rayleigh-Ritz  approach 
is  manifested  in  several  recent  studies  (References  76  through  78)  where 
attempts  have  been  made  to  develop  the  technique  into  a  design  tool. 

2.2.3  Fatigue  Analysis  Methods 

The  phenomenon  of  fatigue  crack  initiation  and  propagation  in 
metallic  structures  has  been  studied  by  many  investigators.  As  a  result, 
several  useful  and  practical  damage  propagation  models  have  been  established 
on  the  basis  of  classical  linear  elastic  fracture  mechanics.  One  such  ex¬ 
ample  is  the  Forman  equation  (Reference  79)  which  is  useful  in  predicting 
the  fatigue  crack  growth  life  of  metallic  structures.  Before  applying 
these  techniques  to  postbuckled  metal  panels,  however,  test  data  are  re¬ 
quired  to  determine  the  fatigue  failure  modes  of  these  panels.  In  addi¬ 
tion,  analytical  techniques  are  needed  to  predict  the  local  stress  intens¬ 
ity  factors. 
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The  only  fatigue  test  data  available  in  the  literature  for  post- 
buckled  metal  shear  panels  are  shown  in  Figure  2.1.  These  data  were  obtained 
in  Reference  3  from  constant  amplitude  fatigue  tests  on  multibay  shear  panels. 
The  four  panels  tested  in  Reference  3  exhibited  very  short  fatigue  lives  (500 
to  4000  cycles)  due  to  cracks  initiating  at  the  corners  of  the  chem-railled 
stiffener  attachment  pads  on  the  skin.  These  data  show  that  fatigue  is  a 
serious  concern  in  the  design  of  postbuckled  metal  panels.  Additional  test 
data  are  required,  however,  to  establish  the  fatigue  failure  modes  and  S-N 
curves  for  curved  metal  shear  panels. 

Compression  fatigue  test  data  for  flat  stiffened  panels  loaded 
in  the  postbuckling  range  have  been  obtained  in  References  80,  81  and  82. 

In  these  panels  fatigue  cracks  occurred  in  the  stiffeners  at  stiffener  attach¬ 
ment  fastener  holes  and  propagated  along  the  loading  direction  as  illustrated 
in  Figure  2.2.  The  fatigue  failure  mode  for  flat  stiffened  panels  loaded 
in  compression,  however,  is  unique  to  this  design.  Crack  initiation  in  the 
skins  at  these  fastener  holes  is  also  possible  depending  on  the  local 
stresses  in  the  skin  and  in  the  stiffener.  Thus,  to  interrogate  all  possible 
modes  of  fatigue  failure  in  postbuckled  metal  panels,  additional  tests  on 
different  designs,  including  curved  panels,  need  to  be  conducted. 

In  order  to  develop  a  generally  applicable  fatigue  life  predic¬ 
tion  methodology  for  postbuckled  metal  shear  and  compression  panels,  analyt¬ 
ical  techniques  that  can  predict  the  local  stress  intensity  factors  are  re¬ 
quired.  This  in  turn  requires  a  knowledge  of  the  detailed  stress  field  in 
the  skins  and  the  stiffeners.  As  mentioned  before,  the  local  stress  field 
can  only  be  obtained  from  nonempirical  analyses.  Thus,  a  Rayleigh-Ritz 
type  analysis  in  conjunction  with  a  fatigue  crack  growth  law  such  as  that 
given  by  the  Forman  equation  can  be  readily  used  to  predict  the  fatigue  life 
of  postbuckled  metal  panels. 

In  contrast  to  the  state-of-the-art  in  fatigue  analysis  of  metals, 
fatigue  analysis  of  composites  is  still  in  its  infancy.  However,  a  sizeable 
fatigue  test  deta  base  for  postbuckled  composite  panels  is  available  and  in¬ 
dicates  that  composite  panels,  in  general,  are  extremely  durable. 
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Figure  2.1.  Metal  Shear  Panel  Fatigue  Data  (Reference  3) 
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The  fatigue  test  data  generated  in  some  of  the  investigations 
cited  in  Table  2.1  provide  a  good  Insight  into  the  durability  characteris¬ 
tics  of  postbuckled  composite  designs. 

Test  data  for  shear  panels  subjected  to  constant  amplitude  and 
spectrum  fatigue  loading  are  shown  in  Figure  2.3.  These  data  were  obtained 
from  several  test  specimens  in  various  government-funded  programs.  From 
Figure  2.3  it  can  be  seen  that  the  spectrum  fatigue  lives  are  considerably 
longer  than  constant  amplitude  fatigue  lives;  this  illustrates  the  rela¬ 
tively  high  severity  of  constant  amplitude  loading.  The  test  data  from 
Reference  16  appear  to  be  the  lower  bound  for  the  fatigue  data.  The  rela¬ 
tively  steep  S-N  curve  for  these  data  is  due  to  a  design  flaw  at  the  stiff¬ 
ener  skin  junction  where  no  ply  drop-offs  were  included  for  a  smooth  trans¬ 
ition.  This  design  drawback  when  corrected  for  (Reference  83)  yields  a 
fatigue  life  comparable  with  the  test  data  from  References  14  and  23.  The 
lack  of  tapered  stiffener  flanges  in  Reference  16  is  also  responsible  for 
making  the  R  =  0.1  fatigue  data  appear  more  severe  than  the  R  «  -1  data  of 
Reference  14.  In  all  cases,  however,  it  should  be  noted  that  the  fatigue 
endurance  limit  is  at  least  the  design  limit  load.  Thus,  the  data  indicate 
extremely  long  fatigue  lives  at  panel  design  limit  load.  In  these  designs, 
the  panels  were  prevented  from  buckling  during  the  level  flight  condition 
of  a  typical  V/STOL  aircraft.  The  minimum  gage*  requirements  resulted  in 
panel  failure  load  being  much  greater  than  the  required  ultimate  load,  a 
condition  which  is  typical  in  most  aircraft  applications. 

The  fatigue  response  of  composite  compression  panels  is  summa¬ 
rized  in  Figure  2.4.  The  data  indicate  that  extremely  long  fatigue  life 
can  be  expected  for  design  limit  strain  levels  of  2,500  yinches/inch.  Most 
postbuckled  panels  are  buckling-critical  and  not  strength-critical.  The 
current  design  practice  does  not  allow  the  average  compressive  limit  strain 
to  be  higher  than  3,000  yinches/inch.  Thus,  fatigue  for  composite  panels 
may  not  be  a  design  driver. 


*  Minimum  gage  defines  the  minimum  laminate  thickness.  In  current  industry 
practice  minimum  laminate  thickness  ranges  between  0.02  inch  and  0.04  inch. 
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Figure  2.3.  Composite  Shear  Panel  Fatigue  Response 
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Figure  2.4.  Composite  Compression  Panel  Fatigue  Response 
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The  available  test  data  for  composite  panels  loaded  under  com¬ 
bined  compression  and  shear  also  demonstrate  the  same  trend.  The  test  data 
(Reference  25)  presented  in  Figure  2.5  indicate  almost  no  loss  in.  panel 
strength  after  two  lifetimes  of  spectrum  fatigue  with  the  maximum  load  set 
at  71.6  percent  of  the  failure  load.  Preliminary  test  data  from  ongoing 
Navy  programs  (References  28,  32)  also  indicate  similar  trends.  The  domi¬ 
nant  failure  mode  in  all  these  fatigue  tests  was  stiffener/skin  separation 
which  is  a  direct  consequence  of  initiation  and  propagation  of  delaminations 
in  the  skin/stiffener  interface.  Analytical  prediction  of  fatigue  life  of 
postbuckled  composite  panels,  therefore,  requires  a  knowledge  of  the  inter¬ 
facial  stresses  and  a  fatigue  analysis  methodology  that  can  predict  damage 
propagation  in  composites. 

Several  analysis  methods  to  predict  skin/ stiffener  interfacial 
stresses  and  subsequent  interface  failure  have  been  proposed  (References 
20,  73,  83  and  84).  However,  experimental  validation  of  these  methods  has 
not  been  very  successful.  In  Reference  20,  Agarwal  has  proposed  a  stiffener/ 
web  interface  stress  analysis  using  a  two-dimensional  nonlinear  model  of  a 
diagonal  strip  from  a  shear  panel.  The  model  utilizes  a  Rayleigh-Ritz  type 
procedure  to  obtain  the  shear  and  normal  stresses  at  the  skin-stiffener 
interface.  Experimental  validation  of  the  model  was  attempted  by  testing 
metal  coupon  specimens  and  comparing  the  measured  failure  loads  with  pre¬ 
dictions.  The  coupon  tests  showed  good  agreement  with  the  predictions  which 
were  based  on  a  quadratic  failure  criterion.  However,  correlations  with 
data  from  tests  on  stiffened  composite  panels  have  not  been  successful  due 
to  uncertainties  in  the  interface  properties  and  the  validity  of  the 
quadratic  failure  criterion.  The  model  proposed  in  Reference  84  by  Tsai  is 
similar  to  the  beam  model  in  Reference  20,  except  that  Tsai  uses  experi¬ 
mentally  measured  out-of-plane  displacements  to  obtain  the  stiffener  web 
interfacial  peel  stresses. 

In  Reference  83  a  detailed  3-D  NASTRAN  stress  analysis  of  the 
stiffener/web  interface  has  been  performed.  These  results  although  useful 
for  comparison  with  other  simplified  analyses,  have  not  been  experimentally 
validated  and  the  method  itself  cannot  be  used  as  a  cost-effective  design 
tool.  A  more  rigorous  approach  of  first  predicting  the  nonlinear  post- 
buckled  response  of  a  stiffened  composite  panel  and  then  using  the  local 
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Figure  2.5.  Fatigue  Response  of  Flat  Panels  Under  Combined  Loading 
(Reference  25) 


stress  field  in  the  skin  to  predict  the  stiffener /web  interface  stresses  by 
a  linear  analysis  has  been  developed  in  Reference  73.  The  global  nonlinear 
analysis  (Reference  72)  and  the  local  (stiff cner/web  interface)  linear 
analysis  (Reference  73)  are  both  carried  out  using  the  Rayleigh-Ritz  method. 
The  analysis  takes  into  consideration  several  skin/stiffener  interface 
variables  that  are  not  accounted  _or  in  the  simple  beam  models,  and  at  the 
same  time  is  less  cumbersome  to  use  than  a  finite  element  analysis.  However, 
the  results  have  not  been  experimentally  validated.  Thus,  at  present  a 
fully  developed  stiffener/skin  interface  stress  analysis  methodology  is  not 
available  for  use  in  fatigue  analysis  of  composite  stiffened  panels.  The 
approach  proposed  in  Reference  73,  however,  seems  to  be  the  most  promising. 

An  extensive  survey  (Reference  85)  of  the  available  methods  for 
fatigue  life  prediction  of  composites  showed  that  these  can  be  broadly 
classified  as  empirical  techniques,  degradation  models,  and  damage  propa¬ 
gation  models.  A  summary  comparison  of  the  advantages  and  disadvantages 
of  these  methods  is  shown  in  Table  2.2  taken  from  Reference  85.  Among  these 
methpds,  the  damage  propagation  models  appear  to  be  well  suited  to  the 
fatigue  analysis  of  postbuckled  composite  panels.  In  particular,  the  strain- 
energy-release-rate  based  delamination  propagation  model  (Reference  86) 
appears  most  promising.  It  is  necessary,  however,  to  extend  this  model 
for  application  to  postbuckled  composite  panels. 

In  summary,  therefore,  the  observation  of  different  failure 
modes  in  metal  and  composite  postbuckled  structures  and  the  significantly 
different  response  of  composite  structures  makes  it  essential  that  separate 
methodologies  be  developed  to  predict  the  fatigue  life  of  metal  and  com¬ 
posite  panels.  A  key  prerequisite  in  both  cases  is  that  the  local  displace¬ 
ment  or  stress  field  in  the  postbuckled  regime  be  known.  In  addition,  for 
composite  panels  a  validated  methodology  to  predict  the  skin/stiffener 
interface  stresses  ia  required.  Finally,  development  of  a  fatigue  analysis 
methodology  by  coupling  the  local  analyses  with  a  crack  growth  or  delami¬ 
nation  growth  law  needs  to  be  carried  out. 

2.3  SELECTION  OF  ANALYTICAL  METHODS 

A  complete  static  analysis  of  postbuckled  structures  consists 
of  predicting  the  initial  buckling  load  of  the  skin,  failure  load  of  the 
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TABLE  2.2.  SUMMARY  OF  LIFE  PREDICTION  MOOELS  FOR  COMPOSITES 


Technique  or  model 

Method 

Advantages 

Disadvantages 

1 

2 

Empirical 

•  Experimental 

•  Simple  analysis 

•  Extensive  testing 

•  No  general  conclusions 

•  Conservative  design 

Miner's  rule 

•  Linear  cumulative  damage 

•  Simple  analysis 

•  Loading  sequence  not  accounted  for 

•  Poor  corralation 

3 

Linear  ftrength 
degradation 

•  Nonlinear  damage 
accumulation  based  on 
linear  strength  degradation 

•  Simple  analysis 

•  Load  sequence  accounted  for 

•  Assumed  strength  degradation  model 
does  not  agree  with  actual  degradation 
generally  observed  in  composites 

4 

Wear-out  model 

•  Strength  degradation 
faasad  on  fracture 
mechanics  of  metals 

•  Statistical 

•  Fatigue  life  and  residual 
strength  directly  related 
to  static  strength 

•  Extensive  testing  parameters  depend  on 
lammete  and  load  spectrum 

•  Cannot  be  used  for  life  prediction 
from  S-N  curves 

5 

Strength  degradation 
model  and  other 
statistical  modeis 

•  Assumed  strength 
degradation  laws 

•  Statistical 

•  Similar  to  wear -out  mode* 

•  Similar  to  wear-out  model 

6 

Dblami  nation 
propagation  model 

•  Oelammation  propagates 
under  interlaminar 

stresses 

a  Growth-rate  equation 
similar  to  crack  growth 
equation  in  metals 

•  Actual  damage  propagation 
modeled 

•  Constants  depend  on  resin 
system  only 

•  Correlates  data  well 

•  Not  applicable  to  nondelamination- 
prone  laminates 

•  Interlaminar  stress 
computation  time-consuming 

7 

Fracture  mechanics 
de  lamination  model 

•  Relating  delamination 
growth  rate  to  strain 
energy  r- lease  rate 

•  Both  modes  1  and  II 
considered 

•  Similar  to  delamination 
propagation  model 

•  Strain  energy  release  rate 
easier  to  obtain 

•  Oelammation  size  determined 
more  accuratly 

•  Not  applicable  tit  nondelamination- 
prone  laminates 

•  Actual  application  to  life  prediction 
not  investigated 

8 

Intralaminar  cracking 
model 

•  Strain  energy  density 
matrix  cracking 

•  Process  of  matrix  cracking 
can  be  modeled  by  a  single¬ 
strain  energy  density 
parameter 

a  Not  applicable  to  delammation- 
prqne  laminates 

•  Prediction  not  verif-ed  by  data 
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structure  after  skin  buckling.  The  scope  of  this  program  encompassed  cyl- 
indrically  curved  stiffened  panels  loaded  in  compression  or  shear  only. 

As  discussed  in  subsection  2.2,  several  analysis  methods,  rang¬ 
ing  from  closed  form  semiempirical  methods  to  extremely  sophisticated  large 
computer  codes,  are  available  to  predict  the  initial  buckling  loads  of 
curved  panels.  The  main  difficulty  in  accurately  predicting  the  initial 
buckling  load  arises  from  a  lack  of  definition  of  the  exact  boundary  con¬ 
ditions  and  due  to  the  presence  of  structural  imperfections.  In  view  of 
these  uncertainties,  semiempirical  methods  based  on  test  data  are  best 
suited  for  preliminary  analyses.  Furthermore,  due  to  the  lack  of  a  well- 
established,  rigorous  failure  analysis  methodology  for  postbuckled  panels, 
semiempirical  analysis  methods  have  to  be  utilized  for  predicting  the 
strength  of  curved  metal  or  composite  panels.  Thus,  to  meet  the  objectives 
of  the  present  program  semiempirical  analytical  methods  that  are  well  docu¬ 
mented  for  metal  panels  (e.g..  References  2  and  34)  were  selected.  Based 
on  available  test  data  for  composite  panels  loaded  in  shear  or  compression, 
the  semiempirical  analyses  were  modified  for  generic  application  to  com¬ 
posites. 

The  semiempirical  analysis  techniques  selected  are  detailed  in 
the  following  paragraphs.  The  modifications  that  extend  the  applicability 
of  the  analyses  to  composites  along  with  supporting  data  are  also  dis¬ 
cussed  . 

2.3.1  COMPRESSION  PANELS 

Analysis  of  postbuckled  curved  compression  panels  is  performed 
in  the  steps  outlined  below.  The  critical  parameters  in  the  analysis  are 
evaluated  in  terms  of  strain,  since  strain  is  more  convenient  to  use  for 
composite  panels  whereas  for  metal  panels  it  can  be  used  interchangeably 
with  stress.  The  analysis  proceeds  as  follows: 

(a)  Determine  the  buckling  strains  for  all  possible 
modes  of  instability.  These  include:  skin  buck-, 
ling  between  stiffeners,  Euler  buckling  of  the 
stiffened  panel,  and  stiffener  crippxing. 

(b)  Determine  the  failure  load  due  to  Euler  buckling, 
stiffenei  crippling  and  other  modes  of  failure 


► 


It 


i 

1 


peculiar  to  composite  or  metal  panels.  For  com¬ 
posite  panels  the  load  foi  stiffener /web  separa¬ 
tion  mode  of  failure  needs  to  be  calculated  where¬ 
as  for  metal  panels  loads  causing  permanent  set 
in  the  skin  or  the  stiffener  must  bs  calculated. 


(c)  The  load  carrying  capacity  of  the  panel  is  then 
determined  as  the  lowest  of  the  loads  calculated 
in  (b)  above. 


Calculation  of  Skin  Buckling  Strain/Load  -  The  buckling  stress 
for  curved  metal  sheet  panels  can  be  calr^atad  from: 


CR 


K  K2E 
c _ 

12(1-V2) 


(1) 


where , 

Fcr  buckling  stress,  psi 

t  thickness  of  the  skin,  in 

b  stiffener  spacing  measured  between  the 

fastener  lines,  in 

E,v  modulus  and  Poisson's  ratio  for  the 
sheet  material 


K  buckling  coefficient  determined  from 

C  Figure  2.6  (Reference  34  and  35) 


The  theoretical  value  of  K  is  obtained  from  the  buckling  equations  for  thin 

c 

cylindrical  shells  and  is  a  function  of  the  nondimensional  curvature  z  of 
the  panel  expressed  as 


b2  (l-v2)^ 
rt 


where  r  is  the  radius  of  the  cylindrical  panel.  Experimental  data  (Reference 
35)  have  shown  that  K£  is  also  a  function  of  the  r/c  ratio  for  the  panel. 

The  design  curves  cf  Figure  2.6,  obtained  from  test  data,  show  this  depend¬ 
ence  of  K  on  r/t. 
c 

Compression  buckling  strains  for  curved  composite  panels  can  be 
accurately  determined  through  the  uue  of  computer  codes  SS8  (Reference  45) 
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Compressive  Buckling  Coefficients  for  Long  Curved  Flates  (Reference  34) 
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and  BUCLASP-2  (Reference  41),  for  example.  However,  for  an  approximate  cal¬ 
culation  of  the  skin  buckling  strain  in  cases  where  the  stiffener  spacing  is 
realistic,  the  simplified  equation  given  below  can  be  used. 


Dll  +  2 


(D 


12 


+  2D66> 
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E  -  /  2v  E 
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E  E 
xw  yw 
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(2) 


where  D  are  the  terms  of  the  bending  stiffness  matrix  of  the  composite 

skin,  E  ,  E  ,  G  ,  V  and  t  are  the  web  elastic  constants  and  thick- 
xw  yw  xyw  xyw  w 

ness,  respectively,  L  is  the  panel  length,  b  is  the  width  of  the  skin,  R 

w 

is  the  radius  of  curvature  of  the  panel  and  n  and  m  are  integer  coeffi¬ 
cients  representing  the  number  of  half  buckle  waves  in  the  width  and  length 
direction,  respectively.  The  lowest  value  of  strain  for  various  values  of 
n  and  m  represents  the  buckling  strain  of  the  specimen. 


The  effective  width  of  the  skin,  bw>  was  assumed  to  be  equal  to 

the  distance  between  the  two  adjacent  stiffeners  measured  from  one  stiffener 

flange  edge  to  the  next  stiffener  flange  edge  as  shown  in  Figure  2.7.*  Note 

that  b  is  less  than  the  stringer  spacing  b  . 
w  s 

Equation  (2)  was  derived  in  Reference  87  from  the  equations  de¬ 
veloped  for  the  buckling  of  or th ■'tropic  complete  cylinders  by  making  sim¬ 
plifying  assumptions. 


Euler  Buckling  Strain  Calculations  -  The  Euler  buckling  strain 
for  a  stiffened  panel  is  calculated  by  treating  the  panel  as  a  wide  column 
with  the  width  set  equal  to  the  stiffener  spacing.  The  critical  strain  is 
calculated  using  the  standard  column  equation: 


*Note  that  this  definition  of  b^  was  used  initially.  The  test  data  in  Sec¬ 
tion  3  indicate  that  b  should  be  measured  between  stiffener  flange  center- 

w 

lines.  See  Section  5  for  details. 
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(3) 


_E  C  7T2EI 

e<*  =  ITT- 

where,  El  is  the  equivalent  bending  stiffness  of  the  panel,  EA  is  the  equiv¬ 
alent  axial  stiffness,  L  is  the  panel  length,  and  C  is  the  end  fixity  coef¬ 
ficient.  The  fixity  coefficient  depends  upon  the  support  conditions  at  the 
panel  ends.  Most  compression  panels  are  tested  by  flat  end  testing  and  the 
results  obtained  by  using  C  =  4  are  quite  unconservative;  therefore,  a  value 
of  C  =  3  is  recommended.  The  values  of  C  for  other  end  conditions  can  be 
obtained  from  Reference  34  (Section  A18.23). 

Stiffener  Crippling  Strain/Stress  Calculation  -  The  crippling 
strength  of  metal  stiffeners  is  calculated  using  the  well  established  Need¬ 
ham  and  Gerard  methods  documented  in  Reference  34.  In  the  present  program, 
the  Gerard  method  was  used  since  it  is  a  generalization  of  the  Needham  method 
and  was  derived  from  a  broader  data  base.  The  empirical  Gerard  equation  for 
calculating  the  crippling  6tress  for  2  corner  sections,  such  as  the  Z,  J  and 
channel  sections,  is: 


(4) 


where , 

F  =  crippling  stress  for  the  section,  psi 
cs 

F  =  compressive  yield  stress  of  the  material,  psi 

t  =  element  thickness,  in 

2 

A  =  section  area,  in 

A  design  curve  based  on  Equation  (3)  is  shown  in  Figure  2.8  taken  from  Refer¬ 
ence  34.  Additional  crippling  equations  that  apply  to  sections  other  than  2 
corner  sections  are  also  given  in  Reference  34. 
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In  order  to  calculate  the  crippling  strains  for  stiffeners  made 
of  composite  materials,  a  semiempirical  methodology  was  developed  in  the 
program.  The  methodology  consists  of  modelling  the  stiffener  in  terms  of 
interconnected  flat  plate  elements,  calculating  the  initial  buckling  and 
crippling  strains  for  each  element,  and  determining  the  crippling  strain 
for  the  stiffener  as  the  lowest  strain  that  causes  crippling  of  the  most 
critical  element  in  the  stiffener  section.  It  should  be  noted  here  that 
the  absolute  minimum  of  the  crippling  strains  for  the  various  plate  ele¬ 
ments  is  not  necessarily  the  stiffener  crippling  strain;  element  critical¬ 
ity  with  respect  to  stiffener  stability  has  to  be  considered  as  well.  The 
procedural  details  of  this  methodology  given  in  the  following  paragraphs 
provide  additional  clarifications  relating  to  the  determination  of  the  most 
critical  plate  element. 

The  first  step  in  calculating  the  stiffener  crippling  strain  is 
to  model  the  stiffener  as  an  interconnected  assembly  of  plate  elements.  As 
examples,  plate  element  models  of  a  hat-section  and  a  J-section  stiffener 
are  shown  in  Figure  2.9.  The  hat-section  stiffener  is  made  up  of  six  ele¬ 
ments,  whereas,  the  J-section  stiffener  consists  of  five  elements. 

The  crippling  strains  for  the  plate  elements  are  calculated  from 
empirical  equations  of  the  form 


where  > 

e  =  crippling  strain  of  the  plate  clement 
c  s 

Ecr  =  initial  buckling  strain  of  the  plate  element 

e  =  compression  ultimate  strain  for  the  plate  element 
laminate 

a, 8  =  material  dependent  coefficients  obtained  from 

test  data 


Equation  (5)  has  the  same  functional  form  as  that  used  by  Gerard  (Reference 
35)  for  metal  stiffeners.  The  coefficients  a  and  3  depend  on  the  plate  edge 
conditions  and  have  been  obtained  in  References  46  and  47  from  a  large  data 
base  for  plate  elements  that  are  connected  on  both  sides  (e.g.,  elements  2, 

3,  4  and  5  of  the  hat-section  stiffener  shown  in  Figures  2.9).  The  crippling 
strain  for  stiffener  plate  elements  connected  on  both  sides  is  given  by 
(Reference  47) : 


e 

cs 


0. 56867s 

cr 


0.47567 


(6) 


where  e  ,  the  buckling  strain  for  the  plate  element  is  given  by  (Reference 
89): 
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cr 


2tr 


2 

b^t  E 


(V» 


UD22  +  D12  +  “66 


(7) 


In  Equation  (7) 

b  =  plate  element  width 

t  =  plate  element  thickness 

E  =  compression  modulus  of  the  plate  laminate 
along  the  longitudinal  direction 


terms  from  the  laminate  bending  stiff  ness 
matrix,  (i,j  =  1,  2,  6) 


Equation  (7)  applies  to  plate  elements  for  which  the  length  to  width  ratio 
(L/b,  where  L  =  stiffener  length)  is  at  least  4. 


The  crippling  strain  for  plate  elements  that  are  connected  on 
one  side  only  is  calculated  using  the  following  equation: 
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cc 


0.4498S 

cr 


0.72715 


(8) 
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where, 


e 

cr 


12D66  +  2°11 
b2t  E  L2t  E 


(9) 


L  =  length  of  the  stiffener 

with  the  other  nomenclature  remaining  the  same  as  for  Equations  (6)  and  (7) . 

The  coefficients  in  Equation  (8)  were  obtained  by  fitting  Equa¬ 
tion  (5)  to  the  crippling  data  generated  from  tests  on  one-edge  free  plates 
in  References  46  and  47.  Data  for  two  material  systems,  T300/5208  and 
AS/3501  graphite/epoxy,  were  pooled  to  obtain  Equation  (8). 

In  Equations  (6)  through  (8) ,  the  thickness  of  plate  elements 
attached  to  the  skin  is  taken  as  the  sum  of  the  plate  element  and  the  co¬ 
cured  skin  thicknesses.  In  the  case  of  the  hat-section  stiffener,  cripp¬ 
ling  strains  for  plate  elements  representing  the  skin  only,  such  as  element 
5  in  Figure  2.9  are  also  calculated.  Another  consideration  in  calculating 
the  crippling  strain  for  stiffener  flange  elements  attached  to  the  skin  is 
the  choice  of  an  appropriate  element  width.  For  example,  in  most  practical 
designs  the  stiffener  flanges  attached  to  the  skin  are  tapered  by  dropping- 
off  plies  as  shown  in  Figure  2.10  for  a  hat-section  stiffener.  The  flange 
plate  element  width  in  this  case  is  defined  as  the  width  to  the  end  of  the 
taper  with  the  weighted  average  of  the  element  thickness  added  on  to  the 
attached  skin  thickness  to  obtain  the  total  thickness  for  use  in  Equations 
(6)  through  (8) . 

Equations  (6)  through  (8)  are  quite  general  in  nature  and  take 
into  account  ply  composition,  stacking  sequence,  and  material  characteris¬ 
tics.  The  oly  composition,  i.e.,  the  percentages  of  0°,  45°  and  90°  plies, 
is  reflected  in  the  compression  ultimate  strain  e  .  Stacking  sequence 
effects  are  accounted  for  in  the  expression  for  e  where  the  bending 

stiffnesses  D. .  are  used.  The  D. and  e  also  account  for  mechanical 
ij  i j  s  cu 

property  changes  from  one  material  system  to  another.  Use  of  strain 
rather  than  stress  for  crippling  calculations  provides  another  signifi¬ 
cant  advantage  in  that  laminate  non-linearity  (e.g,,  stress-strain  response 
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of  ±45°  laminates)  is  accounted  for  by  way  of  the  compression  ultimate 

strain  e  . 
cu 

Failure  Load  Calculation  -  The  failure  load  for  the  panel  is 
determined  as  the  lowest  of  the  loads  calculated  for  the  various  instabil¬ 
ity  modes  mentioned  above,  for  stiffener-web  separation  in  composite  panels, 
and  for  skin  or  stiffener  yielding  in  metal  panels**  The  methods  for  fail¬ 
ure  load  calculation  are  given  in  the  following  paragraphs. 

Failure  Load  Due  to  Euler  Buck? ing  -  The  failure  load  due  to 
Euler  buckling  is  calculated  using  the  following  equation: 


cr 


(E  A  +  E  b  t  ) 
xs  a  xw  w  w 


(10) 


where, 

E 

e  =  Euler  buckling  strain  determined  using 

cr  Equation  (3) 

E  =  Compression  modulus  of  the  stiflener  in  the 
loading  direction 

A  =  Cross-sectional  area  of  the  stiffener 

s 

E  =  Compression  modulus  of  the  web  (skin)  in  the 
loading  direction 

b  =  Stiffener  spacing 

w 

t  =  Skin  thickness 


Failure  Load  Due  to  Stiffener  Crippling  -  In  order  to  determine 

the  failure  load  due  to  stiffener  crippling,  it  is  necessary  to  determine 

the  load  carried  by  the  stiffener  and  the  panel  web  individually.  The 

load  carried  by  the  stiffener  (P  )  is  determined  as  follows: 

s 

1.  Determine  the  two  lowest  crippling  strains  (e  ,) 

ccl 

and  (ecc2^  a II  the  elements  making  up  the 
cross-section  using  Equations  (6)  through  (8) . 

2.  If  the  element  with  the  lowest  crippling  strain 
(£ccj)  is  normal  to  the  axis  of  least  bending 

stiffness  of  the  cioss- section,  the  stiffener 


50 


and  the 


will  fail  at  a  strain  equal  to  £c„1* 

corresponding  failure  of  the  stiffener  is  given 
by: 


E  A  e  , 

XS  X  ccl 


3.  If  the  element  with  the  lowest  crippling  strain 
is  parallel  to  the  axis  of  least  bending  stiff¬ 
ness  of  the  cross-section,  the  stiffener  will 
carry  additional  load  until  the  second  member 
in  the  cross-section  becomes  critical  due  to 
crippling.  In  this  case  the  load  carried  by 
the  stiffener  is  given  by: 


Pc  =  (EA),  (e 


ccl 


"  ecc2) 


+  e 


cc2 


E  A 
XS  s 


where  (EA)^  is  the  extensional  stiffness  of  the 
member  becoming  critical  first,  and  the  stiff- 

g 

ener  failure  strain  €  =  e  _. 

cc  cc2 


(11) 


(12) 


The  total  load  carried  by  the  panel  is  the  sum  of  the  load 
carried  by  the  stiffener  up  to  crippling  and  the  load  carried  by  the  buckled 
skin.  In  order  to  calculate  the  load  carried  by  the  skin,  the  effective 
width  concept  is  utilized.  The  effective  width  for  metal  panels  is  calcu¬ 
lated  using  the  semiempirical  equation  given  below  (Reference  34) : 


w  *  1.9t 

w 


(13) 


where , 

w  =  effective  width  of  the  skin  after  initial  buckling 

t  =  skin  thickness 
w 

F  „  =  stress  in  the  stringer 

st 

For  composite  panels,  in  the  absence  of  any  other  guidelines. 
Equation  (13)  expressed  in  terms  of  strain  is  used  to  compute  the  effective 
skin  width.  Thus, 

a  '°-5 

w  *  1.9t  (e8)  (13A) 

w 
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strain  in  the  stiffener. 


for  composite  skins  where,  e  *» 

Thus,  the  total  load  carried  by  the  panel  for  a  stiffener  crippling  mode 
of  failure  is  given  by: 


where , 


P  +  P 
s  w 


P  =  load  carried  by  the  nanel  at  stiffener  crippling 
cc 

P  »  stiffener  load  given  by  Equation  (12) 
s 


P  =  load  carried  by  the  skin 
w 


The  load  P  is  calculated  as: 
w 


F  w  t 
cs  w 


1.9t  /EF 

W  C3 


for  metal  panels,  and  for  composite  panels  as: 


=  1.9t2  E  (eS  ) 
w  xw  cc 


Failure  Load  bug  to  Stiffener /Web  Separation  -  Failure  of  com¬ 
posite  stiffened  panels  due  to  stiffener/web  separation  is  a  common  mode  of 
failure  in  the  postbuckling  range.  It  is  extremely  difficult  to  predict 
this  failure,  even  by  using  rather  sophisticated  analysis  methods.  The 
attempts  to  date  on  making  such  predictions  have  been  inconclusive.  A 
simple  empirical  equation  to  predict  such  failure  was  developed  in  this 
program.  The  correlation  of  experimental  data  with  the  predicted  failure 
loads  based  upon  this  equation  is  surprisingly  gocd.  The  empirical  equa¬ 
tion  was  derived  by  analogy  with  tha  crippling  data  for  plates  with  one 
edge  simply  supported  and  one  edge  free.  It  is  hypothesized  that  when  the 
panel  web  strain  reaches  the  crippling  strain  the  interfacial  stresses  be¬ 
come  high  enough  to  cause  failure.  The  equation  should  represent  the  level' 
bound  on  predicted  failure  loads.  Any  attempts  to  improve  the  interface 
(for  example,  by  stitching,  riveting,  etc.)  can  result  in  higher  failure 
loads . 
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(17) 


SS 


£__  (E  A_  +  E  b  t  ) 
SS  xs  S  xw  w  w' 


where. 


'SS 


0.4498  e 


cr 


(H 


.72715 


(18) 


£  =  Failure  strain  for  stiffener/web  separation 

8S 


Pgg  *  Failure  load  for  the  stiffener/web  separation  mode 


the  metal  compression  panel  analysis  methodology  outlined  in  the 
preceding  paragraphs  has  been  experimentally  validated  (e.g..  Reference  35) 
and  is  representative  of  current  usage.  In  the  case  of  composite  panels, 
experimental  validation  was  necessary  before  the  methodology  could  be  used 
in  designing  the  program  test  panels.  Composite  compression  panel  test  data 
available  from  some  of  the  studies  cited  in  Table  2.1  were  utilized  to  vali¬ 
date  the  semiempirical  analysis.  Results  of  the  correlation  between  the 
predictions  and  the*  test  data  are  given  in  the  following  subsection. 

Experimental  Validation  of  Composite  Compression  Panel  Analysis  - 
Experimental  verification  of  the  semiempirical  equatiens  was  accomplished  in 
two  parts:  (i)  test  data  on  stiffeners  of  various  shapes  (e.g.,  hat,  channel, 
2,  cruciform)  were  compared  with  predictions  nu.de  using  Equations  (6)  through 
(8);  and  (ii)  test  data  for  flat  and  curved  stiffened  composite  compression 
panels  were  compared  with  the  initial  buckling  and  failure  strain  predictions. 


Stiffener  local  buckling  and  crippling  test  data  for  channel,  Z, 
hat  cruciform  and  I  grahpite/epoxy  sections  were  obtained  from  References 
46,  47  and  48.  A  summary  comparison  of  the  predictions  with  the  test  re¬ 
sults  is  shown  in  Table  2.3.  A  comparison  of  the  predicted  and  measured 
failure  loads  for  the  stiffeners  as  a  function  of  the  strain  ration  Cp/e  , 

where,  £  is  the  strain  at  failure,  is  illustrated  in  Figure  2.11.  As  seen 
r 

in  the  figure,  a  majority  of  the  test  data  fall  on  or  above  the  peXp/^ani  = 
line,  indicating  conservatism  in  the  analysis  which  is  at  most  25  percent. 

A  few  data  points  <n  Figure  2.11  fall  below  the  peXp/pan}  =  1  line.  How¬ 
ever,  these  data  correspond  to  stiffeners  for  which  the  failure  strain  was 
very  close  to  the  compression  ultimate  strain  of  the  laminates  and  the  fail¬ 
ure  mode  was  column  buckling  rather  than  crippling.  Thus,  the  semiempirical 
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TABLE  2.3.  CORRELATION  OF  SEMI-EMPIRICAL  STIFFENER  CRIPPLING  PREDICTIONS 
WITH  TEST  DATA  (Continued) 
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ANALYTICAL  BUCKLING  STRAIN,  ec  -  ANALYTICAL  FAILURE  STRAIN,  -  LAMINATE  ULTIMATE  STRAIN 


TABLE  2.3.  CORRELATION  OF  SEMI-EMPIRICAL  STIFFENER  CRIPPLING  PREDICTIONS 
WITH  TEST  DATA  (Concluded) 
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-  LAMINATE  ULTIMATE  ALLOWABLE  STRAIN,  e  -ANA-:  ''TICAL  BUCKLING  STRAIN,  eF  -  ANALYTICAL  FAILURE  STRAIN 


Figure  2.11.  Comparison  of  Analytical  and  Experimental  Data 
for  Stiffener  Crippling 
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j.tiffener  crippling  predictive  methodology  is  substantiated  by  test  data  over 
a  wide  range  of  stiffener  geometries  and  laminate  lay-ups.  The  analysis  can 
be  used  as  a  design  tool  with  a  high  degree  of  confidence  due  to  the  built-in 
conservatism. 

A  similar  comparison  of  analysis  and  test  data  was  performed  for 
stiffened  compression  panels.  A  majority  of  the  data  available  in  the  lite¬ 
rature  pertained  to  flat  panels.  Table  2.4  summarizes  the  results  of  analyt¬ 
ical  failure  predictions  and  measured  data  for  a  series  of  flat  stiffened 
panels.  To  evaluate  the  accuracy  of  the  predictions  the  failure  load  data 
were  plotted  as  shown  in  Figure  2.12,  where  the  ratio  of  the  measured  fail¬ 
ure  load  to  the  analytically  determined  failure  load  (Ppxp/-Panp  is  Plotted 
against  the  failure  strain  to  ultimate  allowable  strain  ratio  (eF/ecu) •  T.ie 
data  trend  is  similar  to  that  observed  for  stiffener  crippling  predictions 

in  that  a  majority  of  the  data  fall  on  or  above  the  P  /P  ,  =  1.0  line, 

exp  anl 

indicating  conservatism  in  the  analysis  of  approximately  25  percent.  A 
data  point  corresponding  to  NASA  flat  panels  (Reference  26)  falls  approxi¬ 
mately  20  percent  below  the  analytically  predicted  failure  value.  However, 
this  panel  was  designed  such  that  failure  occurred  simultaneously  with  ini¬ 
tial  buckling  of  the  skin  and  as  such  the  panel  was  not  loaded  into  the 
postbuckling  range.  Secondly,  the  initial  buckling  and  failure  strain  of 
6700  yin/in  is  substantially  greater  than  the  current  design  allowable 
strain  levels  for  strength  critical  parts  which  in  turn  are  higher  than 
the  operating  strain  levels  for  postbuckled  designs.  Thus,  the  semiempiri- 
cal  analysis  methodology  for  composite  compression  panels  is  well  suited  to 
the  design  of  postbuckled  panels  where  the  operating  strain  levels  are  of 
the  order  of  2500-3500  yin/in.  As  a  design  tool,  the  semiempirical  method¬ 
ology  is  somewhat  on  the  conservative  side  and  can  be  used  with  a  high  de¬ 
gree  of  confidence. 

Automation  of  Stiffened  Composite  Panel  Design  Methodology  -  The 
semiempirical  compression  panel  analysis  documented  in  the  preceding  para¬ 
graphs  has  been  used  to  develop  a  computer  program  named  CRIP  to  provide  an 
effective  design  tool.  This  program  is  fully  documented  in  Reference  91 
where  its  use  is  also  demonstrated  by  a  design  example. 
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LAMINATE  ULTIMATE  ALLOWABLE  STRAIN,  €cr  -  ANALYTICAL  BUCKLING  STRAIN,  €f  -  ANALYTICAL  FAILURE  STRAIN 
TEST  DATA  SUPPLIED  BY  CAPT  BECKER 
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Figure  2.12.  Comparison  of  Analytical  and  Experimental  Data 
for  Composite  Compression  Panels 
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The  program  has  been  written  for  interactive  use  and  has  several 
built-in  stiffener  shapes  for  application  to  a  wide  variety  of  designs. 

2.3.2  Shear  Panels 

Flat  or  curved  shear  panel  analysis  is  accomplished  by  means  of 
the  semiempirical  tension  field  theory  developed  by  Kuhn  (Reference  2)  for 
metal  panels.  In  this  program  the  tension  field  theory  was  modified  for 
application  to  composite  shear  panels  by  taking  into  account  material  ani¬ 
sotropy. 

The  essential  elements  of  the  generalized  (for  application  to 
metals  as  well  as  composites)  tension  field  theory  and  its  application  are 
summarized  in  Figure  2.13.  Details  of  the  semiempirical  analyses  required 
to  perform  the  various  steps  in  Figure  2.13  are  given  in  the  following  para¬ 
graphs.  The  equations  as  presented  below  pertain  to  cylindrically  curved 
composite  panels  and  to  flat  composite  panels  if  terms  incorporating  the 
radius  of  curvature  R  are  set  equal  to  zero.  Use  of  the  appropriate  values 
for  elastic  constants  in  the  equations  permits  their  direct  application  to 
metal  panels.  The  analysis  procedure  is  based  entirely  on  the  theory  pre¬ 
sented  in  Reference  2  unless  specifically  noted. 

Computation  of  the  Diagonal  Tension  Factor  -  The  diagonal  ten¬ 
sion  factor  k  characterizes  the  degree  to  which  diagonal  tension  is  developed 
in  the  skin  of  stiffened  panels  loaded  in  shear.  A  value  of  k  =  0  charac¬ 
terizes  an  unbuckled  skin  with  no  diagonal  tension;  a  value  of  k  =  1.0  char¬ 
acterizes  a  web  in  pure  diagonal  tension.  The  diagonal  tension  factor  is 
computed  using  the  following  expression: 


k 


Tanh 


(0.5  +  300 


where, 

t  =  web  thickness 

w 

hr  =  ring  spacing 

h  =  stringer  spacing 

s 


(19) 
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Figure  2.13  Application  of  Tension  Field  Theory  to  Shear 
Panels . 


R  =  panel  radius 

T  =  applied  shear  stress 

Tcr  =  buckling  shear  stress  of  web 

The  shear  buckling  stress  or  strain  for  composite  webs  can  be  calculated 
using  program  SS8  (Reference  45) .  The  buckling  stress  for  curved  metal 
webs  can  be  calculated  using: 

X 

cr, elastic 


where , 

Ksl’Ks2  = 

R 

E 


Z 


V 

Computation  of  Diagonal  Tension  Angle  *q'  -  An  initial  value  is 
assigned  to  the  diagonal  tension  angle  'a'  that  defines  the  angle  of  the 
’folds’  in  the  buckled  skin.  For  curved  web  systems  Ct  =  30°  was  found  to 
be  a  convenient  starting  point.  The  actual  value  of  a  is  determined  by  the 
iterative  procedure  outlined  below. 

Using  the  assumed  initial  value  of  a,  a  'new'  value  for  a  is 
calculated  by  the  equation: 
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critical  shear  stress  coefficients  for 
simply  supported  curved  plates,  given 
in  Reference  2 
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Young's  modulus  for  the  material,  psi 
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Poisson's  ratio  for  the  material 
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For  eccentric  stringers  and  rings 
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(21e) 


In  Equations  (21),  e  is  the  skin  strain  in  the  diagonal  tension  direction, 

and  e  and  e  are  the  strains  in  the  stringer  and  the  ring  leg  attached  to 

the  web  averaged  over  their  lengths,  respectively.  E^,  E^  and  E^  are 

the  web  moduli  in  the  direction  of  the  tension  field,  stringers  and  rings, 

respectively.  G  is  the  web  shear  modulus.  EA  and  EA  are  the  effective 
rs  s  r 
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axial  stiffnesses  of  the  stringers  and  the  rings,  respectively.  El  is  the 
bending  stiffness  about  the  stiffener  neutral  axis  and  El  the  bending  stiff¬ 
ness  about  the  web  midsurfaca. 


Ip  general,  o^,  the  new  diagonal  tension  angle  will  not  equal 
the  initially  assumed  value  of  30°.  Therefore,  is  used  as  the  next  guess 
and  the  computations  of  Equations  (21)  are  repeated  until  the  process  con¬ 


verges  ,■  i .  e . ,  a 


new 


a 


old' 


Once  the  diagonal  tension  angle  has  been  determined  with  suffi¬ 
cient  accuracy,  the  next  step  is  to  compute  the  margins  of  safety. 

Computation  of  Stringer  and  Frame  Margins  of  Safety  -  The  diago¬ 
nal  tension  angle  value  computed  above  is  now  substituted  in  Equations  (21) 
to  obtain  the  diagonal  tension  strain  in  the  skin,  the  stringer  strain,  and 
the  ring  strain.  Next,  the  stringer  and  ring  strains  averaged  over  the 
cross  section  and  the  length  (C  .  )  and  the  maximum  strains  in  the  legs 

c i\  6 

attached  to  the  web  (e  )  are  computed  using  the  following  equations: 

max 
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(22) 


(23) 


(24) 

(25) 
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The  stringer  and  ring  crippling  mode  of  failure  is  then  analyzed  for  by  com¬ 
puting  the  stringer  and  ring  forced  crippling  strains  (e  and  e  ,  respec- 

OS  OXT 

tively)  using  the  following  equations: 


(26) 


(27) 


where  E  is  the  laminate  allowable  strain,  E  and  E  ,  are  the  modulus 
all  cs  cr 

of  the  stringer  and  ring  leg  attached  to  the  web,  respectively,  and  tys  and 
tur  are  the  thickness  of  the  stringer  and  the  ring  leg  attached  to  the  web. 

The  critical  stiffener  strains  corresponding  to  the  bending 
stiffness  required  for  stiffener  stability  are  calculated  using  Equations 
(28)  and  (29). 


(28) 


(29) 


where,  Egg  and  are  the  Euler  buckling  strains  for  the  stiffener  and  the 
ring,  respectively. 

The  margins  of  safety  can  now  be  computed  for  each  of  the  possi¬ 
ble  failure  modes  by  comparing  the  calculated  strain  values  with  the  allow¬ 
ables.  Thus,  to  ensure  positive  margins,  the  following  failure  modes  are 
examined  and  the  corresponding  inequalities  verified. 
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(i) 


For  stringer  and  ring  stability, 
i . e . ,  no  column  failure 


e  „  >  s 
sB  '  s 


(ii)  For  stability  of  the  entire 

panel,  i.e.,  to  prevent  buck¬ 
ling  of  the  web  as  a  whole, 
before  formation  of  the 
tension  field 


ave 


e  _  -s  e 
rB  ^  s 

ave 
-  3h 

El  >  E  tJ  (-^  -2)h 
s  sw  h  s 

r 


,  3h 

El  >  E  t  (~  -2)h 
r  ^  r  w  v  h  '  r 
s 


(30) 


(iii)  For  prevention  of  forced 
crippling  of  stiffeners 


e  •>  e 
os  ^  smax 

£  >  £ 

max 


An  additional  check  needs  to  be  performed  for  metal  panels  where 
yielding  or  permanent  set  in  the  web  is  likely  due  to  excessive  skin  deforma¬ 
tion.  The  only  available  criterion  for  permanent  set  check  has  been  empiri¬ 
cally  obtained  from  tests  on  flat  aluminum  metal  panels.  Its  applicability 
to  other  materials  or  curved  panels  has  not  been  verified.  Thus,  in  the  ab¬ 
sence  of  any  other  guidelines,  the  flat  panel  requirement  that  the  maximum 

allowable  value  of  the  diagonal  tension  factor  k  , ,  be  limited  to 

all 

kall  “  0,78  ~  (t-0.012)0,50  (31) 


at  design  ultimate  load  to  prevent  permanent  buckling  of  the  web  at  limit 
load,  is  used  in  the  present  analysis. 

Experimental  Validation  of  Shear  Panel  Analysis  -  The  semiempir- 
ical  analysis  outlined  above  has  been  experimentally  verified  for  metal 
panels  in  References  2  and  5.  In  order  to  validate  the  modifications  intro¬ 
duced  in  the  methodology  for  anisotropic  materials,  the  analysis  results 
weie  compared  with  available  test  data  for  composite  shear  panels.  The 
analysis  methodology  was  exercised  on  composite  shear  panels  designed,  fab¬ 
ricated  and  tested  in  References  13  through  17.  In  these  studies  a  total 
of  7  panel  configurations  with  different  stiffener  shapes,  web  thicknesses, 
web  laminate  orientations,  and  stiffener  spacings  were  tested  through  the 


use  of  various  test  setups.  Table  2.5  summarizes  some  of  the  key  parameters 
of  the  panels,  and  the  ratio  of  the  analytically  predicted  failure  loads  to 
the  experimentally  observed  failure  loads.  The  analytical  predictions  were 
based  on  stringer  forced  crippling  mode  of  failure.  The  failure  mode  for 
all  panels  tested  was  separation  of  the  stringers  from  the  panel  skin.  The 
good  agreement  between  the  measured  and  predicted  failure  strains  inspite 
of  the  difference  in  failure  modes  indicates  that  the  two  failure  modes  are 
closely  related  and  it  is  hypothesized  that  forced  crippling  of  the  stringers 
in  shear  panels  precipitates  skin/stiffener  separation.  Thus,  the  stiffener 
forced  crippling  criteria  can  be  used  to  predict  stiffener/skin  separation 
in  composite  shear  panels. 


TABLE  2.5.  CORRELATION  OF  SEMHEMPIRICAL  ANALYSIS  WITH 
TEST  DATA  FOR  8HEAR  PANELS 


INVESTIGATOR 

(REFERENCE) 

hs(INCH) 

h  ,h 
r/  s 

STIFFENER 

SHAPE 

EA 

s 

(KSI) 

El 

s 

LB/IN2  X  10 

tult 

Tcr 

PANL/PEXP 

NORTHROP/NAVY 

(14) 

10 

1.5 

HAT 

2.5 

0.40 

5 

0.93 

MCDONNELL /NAVY 
(16,94) 

6 

2.67 

HAT 

1.8 

0.37 

9.2 

0.91 

(1.03)* 

LOCKHEED  NAVY 
(15) 

6 

3.75 

I 

1.5 

0.90 

6 

1.08 

GRUMMAN /NAVY 
(17) 

7 

3.43 
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2.9 

1.0 

6 

1.025 

NORTHROP /IR&D 
(13) 

13 

1.15 
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2.9 

0.73 

10 

1.05 
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(13) 

13 

1.15 
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0.313 

10 

0.91 
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0.80 
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Failure  due  to  ring  crippling 
Stringer  spacing 

Ring  spacing 

Stringer  axial  stiffness 
Stringer  bending 


ULT 


cr 


'  ANL 

'exp 


Ultimate  failure  stress 
Buckling  stress 
Analytical  failure  load 
Experimental  failure  load 


It  3hould  be  noted  that  the  forced  crippling  equations  used  in 
the  analysis.  Equations  (26)  and  (27),  do  not  specifically  include  the  inter¬ 
face  material  properties.  Additional  verification  using  data  generated  for 
a  variety  of  material  systems  is  essential  before  the  application  of  theae 
equations  to  st Lf f ener/web  separation  prediction  can  be  generalized.  Fail¬ 
ure  of  the  stif fener/web  interface  does  not  necessarily  have  to  occur  if 
the  interfacial  strength  can  sustain  the  applied  stresses  due  to  forced 
crippling.  Thus,  the  forced  crippling  criterion  sterns  to  present  a  lower 
bound  for  the  failure  load  of  cocured  composite  stringers  by  stiffener/web 
separation.  This  information  is  of  significance  to  designers. 

All  panels  examined  above  consisted  of  composite  stringers  and 
metal  frames  (rings)  with  the  exception  of  the  panels  in  Reference  16.  The 
ring  crippling  failure  load  and  the  stringer  crippling  failure  load  for  the 
panels  in  Reference  16  were  nearly  equal  in  magnitude  with  analytically 
predicted  failure  due  to  stringer  separation.  Although  the  experimental 
data  showed  failures  due  to  frame  (ring)  separation,  subsequent  efforts  to 
improve  the  strength  of  these  panels  by  improving  the  ring/web  interface 
(Reference  83)  resulted  in  failure  due  to  stringer/web  separation  without 
much  increase  in  the  panel  failure  load,  indicating  both  modes  ot  failure 
to  be  quic;e  close  to  each  other  as  predicted. 

The  analytical  and  experimental  L-nr  relations  presented  above, 
thus,  mark  a  milestone  in  the  analysis  of  pnstbuckled  composite  shear  panels 
since  they  validate  the  modifications  to  the  tension  field  theory. 

Automation  of  Shear  Panel  Analysis  Methodology  -  The  modified 
tension  field  theory  outlined  above  has  been  incorporated  in  a  computer 
program  called  TENWEB  that  can  be  used  as  an  efficient  design  tool.  De¬ 
tailed  documentation  for  this  shear  panel  analysis  program  is  given  in 
Reference  91.  The  program  is  interactive  and  has  several  built-in  stiff ener 
profiles  for  design  flexibility. 

Program  TENWEB  was  used  to  design  the  shear  panels  tested  in 
this  study.  Details  of  program  operation  are  given  in  Reference  91. 
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SECTION  3 


EXPERIMENTAL  PROGRAM 


3.1  INTRODUCTION 

The  technology  assessment  documented  in  Section  2  showed  that  the 
current  design  methodology  for  postbuckled  panels  is  predominantly  empirical 
and  was  originally  developed  on  the  basis  of  flat  metal  panel  test  data.  Al¬ 
though  the  methodology  was  later  extended  to  curved  metal  panels,  the  data 
base  to  verify  this  extension  was  extremely  limited.  Since  the  near  term 
application  of  postbuckled  metal  or  composite  panel  designs  is  expected  to  be 
in  curved  fuselage  structures,  it  is  essential  that  the  differences  in  the 
static  rnd  fatigue  response  of  flat  and  curved  panels  be  understood  and  the 
available  data  base  on  curved  metal  and  composite  panels  be  expanded  for  veri 
ficetion  of  the  analysis  methodology. 

There  are  two  main  differences  in  the  postbuckling  behavior  of 
curved  and  flat  panels.  First,  evidence  exists  to  show  that  the  initial  buck 
ling  load  of  curved  metal  panels,  while  higher  than  flat  panels  of  the  same 
size,  is  reduced  significantly  after  repected  loading.  Test  data  do  not  show 
the  same  phenomenon  for  flat  metal  panels.  Test  data  show  thac  the  buckling 
load  of  composite  flat  panels  is  reduced  due  to  fatigue  loading. 

Second,  the  buckling  of  curved  shear  panels  produces  significant 
inward  normal  forces  on  the  stringers  and  frames.  These  forces  may  prove 
beneficial  for  cocured  composite  panels,  since  they  tend  to  delay  the  sepa¬ 
ration  of  stringers  and  frames  from  the  skin,  which  was  found  to  be  a  pri¬ 
mary  mode  of  failure  in  flat  composite  tension  field  panels.  These  differ¬ 
ences  in  the.  behavior  of  curved  and  fiat  panels  may  have  a  significant 
effect  on  their  ultimate  strength  and  fatigue  life.  Therefore,  in  this 
experimental  program,  curved  metal  and  composite  panels  were  tested  to 
establish  a  reliable  data  base  on  their  static  and  fatigue  response.  As  a 
first  step  in  developing  a  design  methodology  for  curved  postbuckled  panels, 
it  was  also  decided  that  the  panels  would  be  loaded  either  in  compress 
sion  or  shear  only,  so  that  once  panel  behavior  under  these  simpler 
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loading  conditions  was  understood,  the  more  complex  case  of  combined  load¬ 
ing  could  be  addressed  next. 


Thus,  the  specific  objective  of  the  test  program  was  to  conduct 
static  and  fatigue  tests  on  curved  metal  and  composite  panels,  loaded  well 
into  the  postbuckling  range  under  compression  or  shear  so  that  data  could  be 
obtained  to  fill  the  gaps  in  the  current  technology  related  to  postbuckled 
aircraft  structures.  In  order  to  define  a  cosl -effective  test  matrix  for  the 
program,  the  most  crucial  data  requirements  were  identified  from  the  technology 
assessment  of  Section  2.  These  data  gaps  are  summarized  in  Table  3.1  and  were 
used  in  selecting  the  test  matrix.  Selection  of  the  test  specimen  configura¬ 
tion,  and  the  design  criteria  was  based  on  the  geometric  and  loading  conditions 
encountered  in  actual  aircraft  fuselage  construction.  Design  of  the  test 
specimens  is  detailed  in  Section  3.2. 

A  detailed  rationale  for  the  selected  test  matrix  and  the  scope 
of  the  tests  is  given  in  Section  3.3.  The  other  significant  aspects  of  the 
experimental  program  such  as  fabrication  of  the  test  specimens,  the  test  fix¬ 
ture  and  instrumentation  used,  and  the  test  procedure  are  described  in  Section 
3.4  through  3.6.  The  test  data  obtained  are  summarized  in  Section  3.7. 

3.2  DESIGN  OF  CURVED  TEST  PANELS 

Aircraft  fuselage  structural  panels  are  rarely,  if  ever,  of  con¬ 
stant  curvature.  Typical  military  aircraft  fuselage  structures  range  approxi¬ 
mately  between  6  and  20  feet  in  diameter.  Stiffened  panels  used  in  con¬ 
structing  ;  large  diameter  fuselage  have  relatively  mild  curvatures,  and 
flat  panels  can  generally  be  used  to  simulate  f^ir  behavior.  Fuselage 
panels  in  fighter  aircraft  have  considerably  smaller  radii  of  curvature. 

In  order  to  duplicate  the  behavior  of  such  panels  and  to  evaluate  the  effect 
of  curvature  on  postbuckling  behavior,  panels  with  relatively  small  radii  of 
curvature  must  be  tested.  The  panels  selected  for  the  present  test  program 
fall  in  this  latter  category  and  have  a  radius  of  45  inches.  This  radius  of 
curvature  was  selected  since  it  is  representative  of  small  diameter  fuselage 
panels  and  to  enable  demonstration  of  the  most  significant  differences  be¬ 
tween  the  behavior  of  flat  and  curved  postbuckled  panels.  The  results  obtained 
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TABLE  3.1.  BATA  GAPS  IN  THE  CURRENT  TECHNOLOGY  FOR  DESIGN  AND 
LIFE  ANALYSIS  OF  CURVED  POSTBUCKLED  PANELS 
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will  be  directly  applicable  to  a  large  number  of  future  aircraft.  In  addi¬ 
tion,  the  metal  and  composite  panels  are  designed  to  satisfy  the  same  design 
criteria  so  that  their  relative  efficiencies  can  be  compared. 


The  design  of  metal  and  composite  shear  and  compression  panels, 
and  the  design  criteria  are  described  in  this  section.  The  analysis  methodol¬ 
ogy  used  for  this  purpose  has  been  detailed  in  Section  2.  The  resulting  panel 
configurations  were  used  in  the  test  program. 

3.2.1  Design  Criteria 

The  typical  compression  and  shear  loads  acting  on  an  aircraft 
fuselage  panel  can  have  a  relatively  wide  ranfi-  if  values  depending  upon  the 
panel  location  and  the  type  of  aircraft.  Hownver,  panels  allowed  to  buckle 
are  generally  lightly  loaded  and  thus  the  loading  range  is  significantly 
narrowed.  A  limit  load  intensity  range  of  300  to  800  pounds  per  inch  for 
shear  and  compression  panels  can  accommodate  a  large  number  of  fighter  as  well 
as  larger  aircraft.  Recent  studies  conducted  under  Navy  sponsorship  have 
concentrated  on  a  limit  load  intensity  of  400  pounds  per  inch.  In  order  to 
extend  the  range  of  currently  available  experimental  data,  the  panel  configura¬ 
tions  selected  for  this  program  were  designed  for  a  limit  load  intensity  of 
600  pounds  per  inch.  The  panels  are  designed  to  buckle  at  approximately  30 
percent  of  the  limit  load  and  to  withstand  ultimate  load  (1.5  times  design 
limit  load)  without  rupture  or  collapse.  The  design  loads  for  the  metal  and 
composite  panels  are  summarized  in  Table  3.2. 

Material  Selection 

The  composite  shear  and  compression  panels  were  fabricated  using 
a  combination  of  woven  and  unidirectional  graphite/epoxy  materials.  The  woven 
graphite/epoxy  material  selected  was  Hercules  A370-5H/3501-6,  whereas  Hercules 
AS/3501-6  graphite/epoxy  tape  was  used  for  the  unidirectional  material.  These 
material  systems  are  representative  of  the  composite  materials  currently  being 
used  in  fighter  aircraft  structures.  The  mechanical  properties  of  these 
materials  are  summarized  in  Table  3,3. 

The  metal  panels  were  fabricated  using  rolled  aluminum  sheet  and 
extruded  stringers.  The  elloy  used  was  7075  with  a  T6  heat  treatment.  The 
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TABLE  3.2.  DESIGN  LOADS  FOR  METAL  AND  COMPOSITE  TEST  PANELS 


DESIGN  LOADS 

COMPRESSION 

PANELS 

Nx  ,  Ibs/in 

SHEAR  PANELS 
Nxy.  Ibs/ln 

SKIN  BUCKLING  LOAD,  NCf 

200 

200 

DESIGN  LIMIT  LOAD,  DLL 

600 

600 

DESIGN  ULTIMATE  LOAD,  DUL 

900 

900 

TABLE  3.3.  GRAPHITE/EPOXY  MATERIAL  PROPERTIES 


PROPERTY 

AS/3501-6 

A370-5H/3501-6 

(FABRIC) 

-C 

E }  ,  psi 

18.7  x  106 

10.0  x  106 

_c 

E2,  psi 

1.87  x  106 

9.2  x  1C6 

G|2  ,  psi 

0.85  x  106 

0.9  x  106 

".2 

0.3 

0.055 

3 


properties  for  this  material  were  obtained  from  MIL-HDBK-5. 

3.2.2  Curved  Shear  Panel  Design 

A  flow  chart  summarizing  the  design  procedure  for  curved  composite 
and  metal  shear  panels  is  shown  in  Figure  3.1.  The  design  loads  for  the  shear 
panels  are  given  in  Table  3.2.  A  frame  spacing  (h  )  of  24  inches  was  selected 
for  the  shear  panels  since  it  is  representative  of  actual  fuselage  structures. 
The  stiffener  configuration  selection  was  based  on  consideration  of  structural 
efficiency,  manufacturing  feasibility  and  cost  (for  composites),  and  current 
design  practice. 

In  several  recent  studies  hat  section  stringers  have  been  chosen 
for  composite  shear  panels  because  of  their  superior  efficiency.  Design  appli¬ 
cation  and  fabrication  studies  have  shown  that,  due  to  their  higher  torsional 
stiffness  as  compared  to  open  sections,  hat  stiffened  panels  can  be  efficiently 
accommodated  in  fuselage  construction.  Thus,  a  hat  section  stringer  configura¬ 
tion  was  selected  for  the  composite  shear  panels.  The  frame  configuration 
selected  was  a  J-section  since  it  is  relatively  easy  to  fabricate,  while  at 
the  same  time  providing  ease  of  attachment  to  other  substructure. 


For  the  metal  shear  panels,  Z-section  stringers  and  frames  were 
selected  since  they  offer  the  best  cost  and  efficiency  advantages  as  demonstra¬ 
ted  by  their  widespread  use  in  many  existing  aircraft. 

The  overall  shear  panel  configuration  selected  consists  of  three 
stringers  and  two  frames  (rings) .  Sizing  of  the  composite  and  metal  panels  to 
meet  the  design  criteria  was  carried  out  as  follows: 

a.  Determine  the  optimum  stringer  spacing,  and  web 
configuration  to  satisfy  design  buckling  loads. 

b.  Size  stringers  and  frames  to  accommodate  ultimate 
panel  load. 

Details  of  the  procedure  used  are  given  in  the  following  paragraphs. 


74 


75 


Composite  Shear  Panel 

Two  possible  web  configurations  (  45  2>  90,  45  2 )  and  (  45  ,  0, 

90,  0,  45  ),  which  are  efficient  in  the  range  of  design  loads  being  considered, 
were  studied  to  determine  their  initial  buckling  loads.  The  first  configura¬ 
tion  consists  of  four  fabric  plies  and  one  unidirectional  tape  ply  and  has  a 
nominal  thickness  of  0.0572  inch.  The  second  configuration  consists  of  two 
fabric  plies  and  three  unidirectional  tape  plies  and  is  0.0416  inch  thick.  The 
second  ply  skin  configuration  barely  exceeds  the  minimum  gage  that  is  permitted 
in  sound  design  practice. 

The  next  step  in  the  design  procedure  was  to  determine  the  skin 

buckling  load  as  a  function  of  skin  thickness  (tw)  and  the  stiffener  spacing 

(h  )  in  order  to  permit  a  judicious  selection  of  values  for  these  two  parameters, 
s 

For  this  purpose  a  buckling  parameter  A,  equal  to  the  ratio  of  the  calculated 
buckling  load  and  the  design  buckling  load,  was  defined.  The  buckling  load 
was  calculated  using  computer  code  SS8  (Reference  45)  and  the  previously 
selected  frame  spacing  of  24  inches. 

Plots  of  buckling  parameter  "A"  for  the  two  web  configurations  for 
different  widths  are  presented  in  Figures  3.2  and  3.3.  These  plots  were  ob¬ 
tained  clamped  and  simply  supported  for  boundary  conditions  as  illustrated  in 
the  two  figures.  The  cylindrically  curved  edges  for  both  cases  were  clamped. 

In  a  stiffened  panel  the  exact  boundary  conditions  aie  not  known  and  it  is 
common  practice  to  assume  that  the  conditions  are  intermediate  between  the  two 
above  boundary  conditions  to  determine  the  buckling  load  of  the  panel  web  between 
the  stiffeners.  Thus,  the  stiffener  spacing  for  the  two  panel  configurations 
to  satisfy  the  design  buckling  requirements  should  be  10  and  5.25  inches, 
respectively. 

The  panel  web  configuration  with  a  (  45  2>  90,  4£  2)  lay-up  and 
with  the  larger  stiffener  spacing  is  much  more  desirable  than  the  Figure  3.3 
configuration,  since  it  will  result  in  substantial  manufacturing  cost  and  weight 
savings  and,  therefore,  was  selected  for  use  in  this  program. 

In  order  to  size  the  stringers  and  the  frames,  tension  field  theory 
as  applicable  to  composite  panels  was  used.  Details  of  the  tension  field  theory 
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Figure  3,2  Buckling  Load  of  a  Curved  Graphite-Epoxy  Plate 
(  ^45,  2 »  90 »  his  Onder  Shear  Loading. 
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Figure  3.3  Buckling  Load  of  a  Curved  Graphite/ Epoxy 

Plate  Under  Shear  Loading.  (.45./90/0/90445J 
Layup. 


7S 


I 

i  i 

with  an  accompanying  summary  flow  chart  are  given  in  Section  2  (Refer  to 

Figure  2.13).  The  iterative  design  procedure  was  implemented  via  computer  . 

program  TENWEB  (Section  2  and  Reference  91) .  The  program  run  used  to  size  these 

test  panels  is  given  in  Reference  91,  Section  3,  as  Figure  5.  The  resulting 

panel  design  shown  in  Figure  3.4.  The  calculated  design  values  and  failure 

modes  are  summarized  in  Table  3.4.  The  panel  edges  have  been  increased  in 

thickness  to  preclude  any  failures  due  to  load  introduction.  As  mentioned 

before,  the  hat  section  stringers  and  J-section  frames  have  been  used  pri-  ^ 

marily  due  to  their  efficiency  and  lower  fabrication  costs. 

Metal  Shear  Panel 

The  design  of  the  metal  shear  panel  proceeded  exactly  along  the 

same  lines  as  that  of  the  composite  shear  panel  described  above”.  The  frame 

and  stringer  spacing  wer2  selected  to  be  the  same  as  for  the  composite  panel 

(h  =  24",  h  =  10") ,  the  web  thickness  t  =*  0.063"  was  calculated  using 
r  s  w 

Figure  3.5  (Reference  34;  as  necessary  to  prevent  buckling  below  design 
buckling  load. 

The  stringers  and  the  frames  were  sized  using  computer  code  TEN- 
WE3  (Reference  91).  The  resulting  panel  configuration  is  shown  in  Figure  3.6. 

The  fastener  spacing  was  calculated  so  as  to  preclude  inter-fastener  buckling 
(Reference  34)  and  tc  prevent  bearing  failure  near  the  fastener  holes.  HYLOK 
fasteners  were  used  instead  of  rivets  to  reduce  fabrication  costs.  In  order 
to  use  flush  rivets  on  thin  akin  (0.063")  the  skin  has  to  be  dimpled,  whereas 
the  use  of  HYLOK  fasteners  does  not  necessitate  skin  dimpling.  It  should  be 
noted  that  the  panel  edges  were  not  initially  reinforced  since  the  web  thick¬ 
ness  is  sufficiently  large  to  prevent  any  static  failure  due  to  load  intro¬ 
duction.  The  computer  output  for  the  metal  shear  panel  design  is  included 
in  Appendix  A  for  reference  purposes.  The  additional  check  required  for  metal 
panels  where  yielding  or  permanent  set  in  the  web  is  likely  due  to  excessive 
skin  deformation,  was  performed  using  Equation  31.  The  predicted  failure  mode 
for  the  metal  panel  configuration  shown  in  Figure  3.6  was  permanent  set  in  the 
web.  The  calculated  design  values  and  the  failure  load  predictions  for  the 
metal  shear  panel  are  summarized  in  Table  3.4, 
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Curved  Compression  Panel  Design 

Design  of  the  curved  composite  and  metal  compression  panels  was 
carried  out  in  accordance  with  the  flowchart  shown  in  Figure  3.7.  It  is 
worthwhile  to  mention  here  that  a  curved  compression  panel  can  be  used  to 
simulate  the  behavior  of  a  cylindrical  built-up  fuselage  structure  provided 
the  cylindrically  curved  panel  is  of  sufficient  width  and  the  appropriate 
boundary  conditions  are  used.  Guidelines  for  determining  the  panel  width 
and  appropriate  boundary  conditions,  so  that  the  panel  buckling  load  will 
equal  the  buckling  load  of  a  cylinder  loaded  in  axial  compression,  have  been 
presented  by  Sobel  and  Agarwal  (Reference  92) .  It  was  shown  that  a  panel 
enclosed  by  an  angle  which  is  greater  than  100  degrees  results  in  a  buckling 
load  equal  to  the  complete  cylinder  load  for  any  arbitrary  boundary  conditions 
along  the  straight  edges.  At  the  same  time,  a  panel  which  is  enclosed  by  less 
than  20  degrees  results  in  a  much  higher  buckling  load  than  a  complete  cylin¬ 
der.  A  panel  enclosed  by  30  degrees  is  able  to  model  a  complete  cylinder  if 
the  appropriate  boundary  conditions  are  used  along  the  straight  edges, 
namely  SSI  (w  =  M  =  Ny  =  =  0) ,  SS3  (w  -  M  -  u  -  Ny  -  0)  or  CC1  (w  =  0  = 

Ny  =  N  =  0),  A  combination  of  boundary  conditions  SSI,  SS3,  and  CC1  can 
be  obtained  if  one  stringer  is  located  at  each  side  of  the  panel.  Thus,  for 
the  test  program,  a  cylindrically  curved  panel  enclosed  by  at  least  a  30- 
degree  angle  with  one  stiffener  at  each  side  was  used  to  simulate  the  complete 
cylinder  behavior. 

Selection  of  the  stringer  configuration  for  the  composite  panels 
was  based  on  an  experimental  and  analytical  evaluation  of  several  flat  com¬ 
pression  panels  with  different  stiffener  configurations  conducted  in  Reference 
44.  Figure  3.8  taken  from  Reference  44  shows  that  hat  stiffeners  are  the 
most  efficient  stiffeners  for  axially  loaded  panels  (all  panels  assumed  to 
be  buckling  resistant) .  Although  the  panels  in  the  present  program  were  sub¬ 
jected  to  loads  beyond  buckling,  the  stiffeners  carry  a  major  portion  of  the 
load  in  the  postbuckling  range  and,  therefore,  the  efficiency  comparison 
shown  in  Figure  3.8  should  be  equally  applicable.  Thus,  hat  section  stiffener 
configuration  was  selected  for  the  composite  compression  panels.  Z-section 
stiffeners  were  selected  for  the  metal  compression  panels  due  to  their 


Figure  3.7.  Compression  Panel  Design  Procedure 
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Figure  3.8  Compression  Load  Structural  Efficiency 
Comparison  for  Hat-,  J-,  and  Blade 
Configurations . 


widespread  use  in  current  design  practice.  The  frame  spacing  for  the  com¬ 
pression  panels  was  selected  as  20  Inches  based  on  the  typical  spacing  of 
15-20  inches  used  in  stringer  stiffened  fuselage  shells,  and  to  allow  for 
skin  buckling  and  stiffener  crippling  at  loads  reasonably  close  to  the  design 
loads. 

A  four  stiffener,  3-bay  configuration  was  selected  for  the  com¬ 
pression  panels  based  on  test  data  developed  in  Reference  93,  where  panels 
tested  with  two  or  three  stiffeners  resulted  in  poor  agreement  with  analytical 
solutions  due  to  distortion  of  edge  stiffeners.  However,  panels  tested  with 
four  stiffeners  resulted  in  good  agreement  with  the  analytical  solutions. 

The  reason  for  the  good  correlation  in  the  latter  case  was  that  the  distortion 
of  the  edge  stiffeners  did  not  affect  the  panel  center  bay. 

Design  calculations  for  the  composite  and  metal  compression  panels 
are  summarized  in  the  following  paragraphs. 

Composite  Compression  Panel 

In  selecting  the  web  configuration,  design  and  test  studies  con¬ 
ducted  in  References  30  and  31  were  used  for  guidance.  In  the  referenced 
studies  the  stiffeners  were  spaced  relatively  close  together  and  the  panel 
web  was  quite  thin  with  a  (±45,  +45)  lay-up.  In  the  present  program  the 
panel  web  was  made  slightly  thicker  and  the  stiffener  spacing  was  increased 
to  lower  manufacturing  costs  and  improve  panel  efficiency.  A  skin  lay-up  of 
[45  /0/90/0/  45]  with  a  nominal  thickness  equal  to  .0416  inch  was  selected 
in  conjunction  with  a  stiffener  spacing  of  12.2  inches  to  meet  the  design 
load  requirement  for  skin  buckling.  The  end  bays  of  the  four  stiffener, 

3-bay  compression  panel  were  made  narrower  to  preclude  early  failure  in  the 
end  bays  while  at  the  same  time  the  end  bay  width  was  sufficient  to  ensure 
skin  buckling  at  loads  much  lower  than  the  failure  load. 

The  initial  buckling  load  and  the  Euler  buckling  load  of  the 
composite  curved  panels  were  obtained  through  the  use  of  computer  code 
BUCLASP-2  (Reference  41) .  Since  it  is  easier  to  work  with  strain  for  com¬ 
posite  panels,  the  following  discussion  makes  extensive  use  of  strain  rather 
than  stress. 
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The  buckling  strains  obtained  from  the  computer  code  BUCLASP-2  were 
rec’uced  by  25  percent  to  accommodate  lower  buckling  strains  due  to  imperfection 
effects,  a  common  practice  for  metal  panels.  However,  the  percent  knockdown  is 
dependent  on  several  design  parameters  which  are  not  defined  for  composite 
panels;  therefore,  for  an  initial  design  estimate,  the  guidelines  for  metal 
panels  given  in  Reference  34  (see  also  Figure  2.7)  were  used. 

The  crippling  strain  (e  )  is  obtained  using  Equations  5  through  9 
and  the  failure  load  calculated  using  Equations  11  through  14.  The  total  load 
on  the  panel  at  stiffener  crippling  is  the  sum  of  the  loads  carried  by  the 
stiffeners  and  the  web.  In  order  to  obtain  the  load  in  the  web  the  effective 
width  method  was  used.  Equation  13A  given  in  Section  2  was  used  to  obtain  the 
effective  width. 

Detailed  design  calculations  for  the  composite  compression  panel  were 
conducted  using  computer  code  CRIP,  a  sample  run  for  which  is  included  in 
Reference  91,  Section  1,  as  Figure  3.  The  panel  configuration  obtained  using 
the  above  approach  is  shown  in  Figure  3.9.  A  summary  of  the  initial  buckling 
and  final  failure  strain  predictions  for  the  test  specimen  is  shown  in  Table 
3.5. 

Metal  Compression  Panel 

The  metal  panel  configuration  selected  for  preliminary  design  was  of 
the  following  geometry: 

Stiffener  spacing  bg  =  10  inches 

Web  thickness  t  -  0.05  inch 

w 

Panel  length  L  =  20  inches 

The  local  buckling  stress  for  the  center  bay  web  between  stiffeners 
was  calculated  using  the  equations  given  in  Reference  34.  In  accordance  with 
Reference  34,  the  buckling  stress  F  for  a  curved  sheet  in  compression  with 
simply  supported  boundary  conditions  is  obtained  as 


where 


cr 


K 


K  _2  E 

12(1 

13 


,2  E  /  t  \2 

7T _  /  W  \ 

Cl-V2)  \  bs  / 
(from  Figure 


2.6) 


F  -  3143  psi 
cr  r 
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Figure  3.9.  Composite  Compression  Panel  Configuration 


The  design  buckling  load  (N  ),  however,  *»  200  lb/in.  Hence  the  stiffener 

x }  cr 

area  (Ag)  required  is  obtained  as  follows: 


N 


x.cr  s 


rcr 


-  b  t 
s  w 


A  ■  .1363  inch 
s 


Assume  the  stiffener  configuration  to  be  AND  10138-1004  which  is  shown  in 

Figure  3.10.  Thus,  the  predicted  panel  buckling  load  N  ■  206  lb/in  which 

xy,cr 

results  in  a  3  percent  margin  of  safety. 

Failure  analysis  of  the  panel  was  carried  out  in  accordance  with 

the  procedure  outlined  in  Section  2.  Figure  2.6  was  used  to  calculate  the 

crippling  stress  F  for  the  stiffener  and  yielded: 
cs 


F  »  48.9  ksi 
cs 


The  effective  web  width  at  the  time  of  stiffener  crippling,  w, 
was  calculated  from  Equation  13A  as: 


w  =  1.2  inch 

The  total  load  at  panel  failure  P  is  calculated  using  Equation  (14)  which 
yields : 


P 


ult 


F 

cs 


(AS  +  wtw} 


-  48900  (.155  +  1.2  x  .05) 

P  ,  «  10500  lb. 

ult 

Hence,  the  ultimate  failure  load  per  unit  width  (N  x  . 

W  18 


\u  ■  Hr-  ■  1050  lb/ln 

s 

Thus,  the  panel  failure  load  allows  approximately  a  15  percent  margin  of 
safety. 
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0.75"  — | 


stiffener  area  ■  0.155  sq.  Inch 

4 

stiffener  M.O.I.  I  -  0.0236  Inch 

“i 

Figure  3.10  Z-Section  7075-T6  Aluminum  Stringer. 

AND  10138-1004  Configuration. 


For  overall  panel  instability  the  Euler  buckling  stress  was  cal¬ 
culated  using  Equation  3  in  the  following  form: 

2 

7T  El 
_ _ e 

Euler  T  2  A 
L 

e  t 

where,  Lg  is  the  effective  length  of  the  panel,  Afc  is  the  total 
area  of  the  panel  and  I  is  the  panel  moment  of  inertia  about  the  neutral 
axis.  Since  the  frame  spacing  for  design  purposes  was  assumed  to  be  20 
inches,  the  effective  length  "L"  for  Euler  buckling  is  10  inches  (C  =  4  in 
Equation  3)  assuming  fully  fixed  ends.  Thus,  the  calculated  Euler  buckling 
stress  for  the  panel  was: 

cl,  ..  *  90.63  ksi 

Euler 

The  actual  Euler  buckling  stress  will  be  lower  than  the  value  above  due  to 
eccentricity  effects  and  yielding,  but  it  is  still  well  in  excess  of  the  panel 
crippling  load.  Thus,  the  panel  design  meets  all  the  required  design  criteria. 

A  sketch  of  the  resulting  metal  compression  panel  configuration 
is  shown  in  Figure  3.11.  The  fastener  pitch  and  other  related  details  shown 
were  obtained  to  prevent  inter-fastener  buckling  and  failures  near  fastener 
holes  (Reference  34) . 

3.3  TEST  PLAN 

The  selection  of  a  cost-effective  and  suitable  test  matrix  for 
the  program  was  made  after  a  review  of  the  gaps  in  the  current  technology  that 
are  summarized  in  Table  3.1. 

The  principle  objectives  of  the  curved  panel  tests  were  to  gener¬ 
ate  static  test  data  for  analysis  methodology  verification  and  to  conduct 
fatigue  tests  for  failure  mode  identification  and  development  of  fatigue 
analysis  procedures.  One  key  concern  addressed  in  the  program  is  the  fatigue 
response  of  curved  metal  panels.  In  addition,  previous  studies  have  indicated 
that  metal  shear  panels  cannot  survive  constant  amplitude  fatigue  test  at 
maximum  loads  much  higher  than  approximately  50  percent  of  their  static 
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Figure  3.11  Metal  Compression  Panel  Configuration. 
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ultimate  strength.  Test  data  show  that  composite  compression  panels,  on  the 
other  hand,  can  endure  maximum  fatigue  load  amplitudes  as  high  as  7u  percent 
of  their  static  ultimate  strength  without  any  strength  degradation.  However, 
since  fatigue  data  on  composite  compression  and  shear  panels  are  sparse, 
additional  fatigue  tests  on  these  panels  were  required.  A  greater  emphasis 
was  placed  on  curved  shear  panels  in  this  program  due  to  their  anticipated 
sensitivity  to  fatigue  loading.  Static  tests  on  the  curved  panels  were  con¬ 
ducted  to  obtain  skin  buckling  strains,  sriffener  crippling  strains,  skin 
and  stiffener  strains  at  stiffener/web  separation,  and  the  strain  dxsLribu- 
tion  in  the  panels  as  a  function  of  the  applied  load.  These  data  were  re¬ 
quired  for  comparison  with  predictions  made  using  the  semiempirical  analysis 
and  the  strain  distribution  in  particular  for  comparison  with  the  non-empiri- 
cal  Rayleigh-Ritz  analysis  procedure  developed  in  the  program.  Additional 
details  of  the  test  program  are  given  in  the  following  paragraphs. 

3.3.1  Test  Matrix 

The  program  test  matrix  is  shown  in  Table  3.6.  As  indicated  in 
the  table,  a  total  of  26  panels  were  tested  in  the  program.  Four  sets  of 
panels  were  tested  to  obtain  the  initial  buckling  load,  postbuckling  behavior, 
ultimate  failure  load,  and  mode  of  failure.  The  four  sets  consisted  of: 
aluminum  compression,  aluminum  shear,  graphite/epoxy  compression,  and  graphite/ 
epoxy  shear  panels.  All  tests  were  conducted  in  a  room  temperature  dry  (RTD) 
environment.  A  greater  emphasis  is  placed  on  shear  panel  fatigue  tests  since 
fatigue  data  for  curved  composite  shear  panels  are  not  available  and  those 
for  metal  panels  are  limited  in  quantity.  Each  test  condition  was  replicated 
twice  to  demonstrate  the  repeatability  of  the  test  and  to  obtain  more  reliable 
test  data  for  analysis  verification. 

The  constant  amplitude  fatigue  tests  on  compression  panels  were 

conducted  at  an  R-ratio  (a  /a  )  of  10  with  the  maximum  fatigue  load  set 

min  max 

at  66  percent  of  the  static  strength  for  the  first  metal  panel  test  and  at 
70  percent  of  the  static  strength  for  the  first  composite  paneL.  Selection 
of  the  maximum  fatigue  loads  for  the  subsequent  tests  was  made  on  the  basis 
of  the  measured  panel  response  for  the  first  set  of  tests  and  the  need  to 
obtain  a  definition  of  the  S-N  curve  for  the  panels. 
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The  aluminum  shear  panel  fatigue  tests  were  also  conducted  under 
constant  amplitude  loading  with  R  =  0.  The  fatigue  load  amplitudes  selected 
for  these  tests  were  55  percent  and  45  percent  of  the  average  static  ultimate 
strength  for  similar  panels.  In  the  case  of  the  composite  shear  panels  the 
fatigue  tests  were  conducted  at  R  =  0.25  and  R  *  -0.25,  with  the  latter  R-ratio 
allowing  for  partial  reversal  of  the  shear  loading.  The  fatigue  load  ampli¬ 
tudes  for  the  composite  shear  panels  were  higher  than  those  used  for  the 
aluminum  panels  (70  percent  and  60  percent  of  the  average  static  strength  for 
the  panels)  due  to  their  much  superior  fatigue  response. 

The  test  panel  instrumentation  consisted  mainly  of  strain  gages 
and  in  some  select  cases  of  displacement  transducers .  A  more  complete  de¬ 
scription  of  the  instrumentation  is  given  in  paragraph  3.3.2.  It  is  noted 
here,  however,  that  prior  to  and  periodically  during  the  course  of  the  fatigue 
tests,  strain  surveys  up  to  the  maximum  fatigue  load  were  conducted  on  all 
fatigue  test  panels.  As  indicated  in  Table  3.1,  the  intermediate  strain  sur¬ 
veys  during  the  fatigue  tests  were  conducted  at  50,000  and  100,000  cycles. 

3.3.2  Instrumentation 

All  panels  in  the  test  program  were  instrumented  with  strain 
gages.  The  static  compression  test  panels  were  instrumented  with  LVDT's  in 
addition  to  the  strain  gages  so  that  out  of  plane  displacements  could  be 
monitored  during  the  course  of  the  tests.  Figure  3.12  shows  the  layout  of 
the  strain  gages  and  the  locations  of  the  out-of-plane  displacement  transducers 
for  the  compression  panels.  In  this  figure  the  gage  layout  for  the  less  ex¬ 
tensively  instrumented  compression  fatigue  panels  is  also  shown.  The  strain 
gage  layout  for  the  static  and  fatigue  tested  shear  panels  is  shown  in  Figure 
3.13.  As  noted  in  Figures  3.12  and  3.13,  all  gages  were  located  back-to-back 
on  the  convex  and  the  concave  surfaces  of  the  panels  in  order  to  determine 
bending  as  well  as  membrane  strains  due  to  postbuckling  deformations. 

In  all  static  tests,  a  visual  indication  of  the  out-of-plane 
displacements  in  the  postbuckling  regime  was  obtained  by  means  of  the  Moire' 
grid  technique.  For  this  purpose  the  composite  specimens  were  painted  white 
and  a  Moire'  grid  placed  within  0.25  inch  of  the  specimen  surface  in  the  case 
of  both  metal  and  composite  panels  to  obtain  the  fringes  associated  with  the 
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with  gages  ®  through  ®  only. 

Fatigue  test  panels  MC4,  CC5,  and  CC6  were  instrumented 
with  gages  ©  through  (g)  only . 


-r 


buckle  pattern. 

During  the  course  of  the  fatigue  tests  some  of  the  strain  gages 
were  damaged  due  to  fatigue,  thus  causing  loss  of  some  data.  However,  a 
majority  of  the  strain  gages  were  not  affected  and  valid  data  were  obtained 
at  the  50,000  cycle  and  100,000  cycle  strain  surveys.  The  gages  that  survived 
and  the  detailed  nomenclature  are  noted  in  Appendix  A  along  with  the  test  data. 

3.3.3  Test  Fixture 

The  static  and  fatigue  tests  on  the  compression  panels  were  con¬ 
ducted  in  a  100,000  pound  capacity  Tinius  Olsen  test  machine.  The  panel  ends 
were  potted  in  an  epoxy  compound  for  load  introduction.  A  full  view  of  the 
test  setup,  including  the  Kaye  data  acquisition  system,  is  shown  in  the  photo¬ 
graph  of  Figure  3.14.  A  close-up  view  of  a  composite  compression  panel  in  the 
test  machine  with  full  instrumentation  is  shown  in  Figure  3.15. 

The  test  fixture  used  for  curved  shear  panel  tests  is  a 
fixture  designed  and  developed  at  Northrop  under  Independent  Research  and 
Development  funds.  This  fixture  results  in  the  application  of  extremely  uni¬ 
form  shear  stress  in  the  panel  with  no  adverse  stresses.  The  loading  mechanism 
and  the  test  fixture  are  schematically  illustrated  in  Figure  3.16.  A  view 
of  the  test  fixture  with  a  metal  panel  installed  for  testing  is  shown  in  the 
photograph  of  Figure  3.17.  The  glossy  appearance  of  the  specimen  is  due  to 
the  Moire'  grid  which  has  been  positioned  close  to  the  specimen  surface. 

The  curved  panel  is  enclosed  by  two  flat  dummy  panels  making  up 
a  triangular  tube.  The  two  flat  panels  are  considered  part  of  the  test  fix¬ 
ture  and  are  connected  at  a  point  midway  between  the  test  panel  center  of 
curvature  and  the  test  panel.  One  end  of  the  tube  is  clamped  against  all 
degrees  of  freedom.  The  other  end  is  connected  to  the  loading  frame  plate. 
Loading  frame  support  plates  are  slotted  to  allow  free  rotation  of  the  load¬ 
ing  frame  shaft  about  the  tube  centroid.  The  shear  load  is  applied  by  a 
torque  introduced  by  two  load  cylinders  moving  in  opposite  directions.  The 
two  torque  application  cylinders  are  each  of  50,000  pound  capacity  and  the 
torque  arm  is  74.0  inches.  Operation  of  the  test  fixture  was  verified  using 
an  aluminum  panel  which  was  tested  to  failure.  The  fixture  design  permits 


100 


101 


Figure  3.16.  Shear  Panel  Test  Fixture  Schematic 
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3.17.  View  of  Test  Fixture  Showing  Loading  Arrangement 
and  Test  Panel. 
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testing  of  curved  panels  in  a  wide  range  of  sizes  and  curvatures  with  mini¬ 
mal  alterations. 

3.4  FABRICATION  OF  COMPOSITE  TEST  PANELS 

Special  testing  and  curing  fixtures  were  designed  and  built  to 
fabricate  the  composite  shear  and  compression  panels.  The  tooling  for  the 
shear  panels  was  more  complex  due  to  the  curvature  of  the  J-section  frames. 

For  the  shear  panels  a  steel  template  formed  to  a  45-inch  radius  was  used  to 
lay  up  the  skin.  The  template  was  marked  to  indicate  the  peripheral  net  trim 
area;  stringer,  frame  and  doubler  (edge  as  well  as  under  the  stiffener)  loca¬ 
tions.  An  orientation  rosette  was  also  marked  on  the  template  to  ensure  angu¬ 
lar  accuracy  of  the  various  plies. 

The  stiffeners  were  shaped  and  fabricated  using  Dow  Corning 
Silastic-J  RTU  rubber  mandrels.  The  mandrels  themselves  were  cast  in  sheet 
metal  molds  with  the  frame  mandrel  molds  designed  to  allow  for  opening  and 
shifting  of  hat  cavities  when  the  mandrel  is  bent  into  a  45-inch  radius  after 
being  cast  straight.  After  cure  of  the  frame  mandrel,  slots  were  cut  into  it 
to  allow  for  expansion  of  the  rubber.  The  cauls  for  the  stringers  were  fab¬ 
ricated  using  two  plies  of  graphite/ epoxy  sandwiched  between  two  layers  of 
Air-Tech's  Airpad  black  rubber.  These  cauls  were  in  three  pieces  of  which 
two  were  used  on  the  short  stringer  ends  outside  of  the  frames  (Refer  to 
Figure  3.4)  whereas  the  third  one  was  used  for  the  portion  of  a  stringer 
between  the  frames.  The  cauls  on  the  stiffeners  extended  to  the  surface  of 
the  frames  and  were  tapered  to  prevent  excessive  mark-off.  This  kept  the  hat 
stiffeners  straight,  and  prevented  them  from  rolling,  bowing  and  distorting 
during  cure.  Fiberglass  cauls  were  fabricated  for  the  top  of  the  J-section 
frames.  These  kept  the  frames  circumferentially  straight,  and  eliminated 
wrinkles  in  the  cap  of  the  frame.  Figure  3.18  shows  the  stringer  and  curved 
frame  mandrels  in  place  on  the  skin  template.  The  graphite-epoxy  cauls  used 
to  compact  the  stringers  are  also  shown  in  the  figure.  The  cuts  in  the  frame 
mandrels  to  accommodate  the  hat  section  stringers  are  illustrated  in  Figure 
3.19  which  shows  a  photograph  of  the  partially  laid-up  frame  mandrel. 
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Figure  3.18.  Intermediate  Step-in  Composite  Shear  Panel  Fabrication  Illustrating 
Frame  and  Stringer  Mandrels  Located  on  Skin  Template 


Figure  3.19.  Detail  View  of  Frr.ue  Mandrels 


The  panel  -  /^rication  procedure  consisted  of  laying  up  the  skin 
and  the  stringers  separately  on  a  flat  template  and  then  locating  the  preformed 
stringers  onto  the  skin.  The  subassembly  was  debulked  under  vacuum  for  30 
minutes,  then  placed  on  the  curved  steel  template  used  as  the  curing  fixture 
and  taped  in  place  to  avoid  movement  in  subsequent  operations.  The  frames  were 
laid  up  on  their  mandrels  and  located  over  the  subassembly  in  the  curing  fix¬ 
ture.  The  stringer  cauls  were  then  installed,  followed  by  installation  of 
the  frame  cauls.  The  panel  assembled  up  to  this  stage  is  illustrated  in 
Figure  3.20.  '.his  assembly  was  then  covered  with  bleeder  and  breather  plies 
as  required  and  bagged  for  cure.  The  panel  was  cured  and  postcured  in  accord¬ 
ance  with  Northrop  specification  MA-133.  The  cured  composite  shear  panel  is 
shown  in  Figure  3.21. 

Fabrication  of  the  composite  compression  panel  was  considerably 
simpler  due  to  the  absence  of  the  curved  frames.  The  procedure  followed  in 
fabricating  these  four  stringer  panels  (Refer  to  Figure  3.9)  was  identical  to 
that  used  for  the  shear  panels  up  to  the  stringer/skin  subassembly  stage.  The 
finished  composite  compression  panel  is  shown  in  Figure  3.22.  The  panels  then 
are  potted  for  compression  load  application. 

The  composite  compression  and  shear  panels  were  nondestructively 
inspected  by  means  of  ultrasonic  C-scan  to  ensure  defect-free  panels. 

3.5  TEST  RESULTS 

3.5.1  Compression  Panel  Static  Tests 

The  composite  and  metal  compression  panel  static  test  results 
are  summarized  in  Table  3.7.  The  metal  compression  panels  MCI  and  MC2  failed 
due  to  stiffener  crippling.  Web  buckling  and  inter-rivet  buckling  was 
observed  in  some  areas  prior  to  failure.  These  static  test  results  are 
compared  against  predictions  in  Section  5.  Development  and  progression  of 
the  buckle  pattern  for  the  metal  panels  is  illustrated  in  Figure  3.23a 
through  e  where  Moire  grid  pattern  photographs  for  panel  MCI  are  shown.  Just 
beyond  the  initial  buckling  load  the  web  is  seen  to  have  buckled  in  two 
half  waves  along  the  load  axis  as  well  as  across  the  width.  This  pattern 
becomes  more  easily  visible  as  the  load  is  increased  further.  However, 
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Figure  3.21.  Cured  Composite  Shear  Panel 


Figure  3.22. 


(a)  Load  -  0  lb 


(b)  Load  =  16K  lb 

Figure  3.23.  Progression  of  Buckle  Pattern  with  Load  for 
Metal  Compression  Panel  MCI 
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(c)  Load  =  25K  lb 


ft 


Po^tbnck ! *d  A'.'j»'.fr>w  Far?: 


(d)  Load  =  3 OK  lb 

Figure  3.23.  Progression  of  Buckle  pattern  with  Load  for 
Metal  Compression  Panel  MCI  (Continued) 


f  Buckle  Pattern  with  Load  for 
ission  Panel  MCI  (Concluded) 


between  25,000  and  30,000  lb*  the  pattern  changes  to,  one-half  wave  across 
the  width  of  the  panel  which  remains  unchanged  to  failure.  A  photograph 
illustrating  the  final  failure  mode  of  the  panel  by  stiffener  crippling  is 
shown  in  Figure  3.24.  Strain  data  for  the  two  metal  panels  are  given  in 
Appendix  A  and  were  used  to  determine  the  buckling  loads  given  in  Table  3.8. 

The  composite  compression  panels  CC1  and  CC2  failed  by  stiffener/ 
web  separation.  Final  failure  was  abrupt  and  resulted  in  the  separation  of 
all  stiffeners  from  the  skin  and  in  secondary  rupture  of  all  stiffeners. 
Development  and  progression  of  the  buckle  pattern  in  panel  CC1  is  shown 
in  Figure  3.25a  through  f.  Panel  CC2  also  displayed  the  same  buckle 
pattern.  The  strain  data  given  in  Appendix  A  were  used  to  determine  the 
panel  buckling  loads.  The  failure  mode  for  these  panels  is  illustrated  In 
Figure  3.26. 

3.5.2  Compression  Panel  Fatigue  Tests 

The  metal  and  composite  compression  panel  fatigue  tests  were 
conducted  at  an  R-ratio  of  10  .  The  test  data  are  summarized  in  Table  3.8. 
Metal  panel  MC3  tested  under  constant  amplitude  loading  with  the  maximum 
fatigue  load  set  at  61  percent  of  the  average  static  strength,  developed 
sizeable  cracks  (~2.5  in)  after  16,000  cycles.  The  cracks  in  the  skin  were 
parallel  to  the  stiffener  and  appeared  to  be  caused  by  the  web  bending 
against  the  stiffener.  Unlike  the  failure  mode  shown  in  Figure  2.2  for 
flat  compression  panels,  the  skin  cracks  in  the  curved  panel  were  located 
along  the  stiffener  edge  away  from  the  fasteners. 

Panel  MC4  tested  at  a  load  amplitude  equal  to  51  percent  of  the 
average  static  strength  also  developed  cracks  in  the  skin  along  the  stiffe¬ 
ners  at  43,000  cycles.  The  panel,  however,  retained  its  fatigue  load 
carrying  capacity  to  100,000  cycles.  The  initial  2.5-inch  crack  grew  to 
4.75  inches  in  the  last  57,000  cycles.  The  panel  was  statically  tested  for 
residual  strength  but  did  not  show  an  appreciable  loss  in  strength.  The 
statically  failed  panel  with  the  initial  fatigue  cracks  marked  is  shown  in 
Figure  3.27. 
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(a)  Skin  Rupture 


(o)  Stiffener  Crippling 


Figure  3.24.  Failed  Metal  Compression  Panel  MCI 
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TABLE  3.8.  COMPOSITE  AND  METAL  COMPRESSION  PANEL  FATIGUE  TEST  RESULTS 
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Additional  strain  surveys  at  150K  and  200K  cycles. 

R  denotes  run-out.  ++  Panel  tested  under  Northrop  IRAD 

Prefix  M  in  panel  number  denotes  metal  panel  and  prefix  C  denotes  graphite/ epoxy  level 


(b)  Load  =  18K  lb 

Figure  3.25.  Progression  of  Buckle  Pattern  with  Load  for 
Composite  Compression  Panel  CC1 


v 


w 


(c)  Load  =  25K  lb 


(d)  Load  =  35K  lb 


Figure  3.25.  Progression  of  Buckle  Pattern  with  Load  for 
Composite  Compression  Panel  CCl  (Continued) 
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(e)  Load  =  60K  lb 


(f)  Load  =  8 OK  lb 


Figure  3.25.  Progression  of  Buckle  Pattern  with  Load  for 
Composite  Compression  Panel  CC1  (Concluded) 


(a)  Stiffener /Web  Separation  and  Failure 


(b)  Skin  Failure 

Figure  3.26.  Composite  Compression  Panel  CC1  Failure  Mode 
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A  curved  metal  compression  panel  with  rounded  stiffener  edges 
( I  Cl)  was  fatigue- tested  under  constant  amplitude  loading  with  R  =  10  and  at 
a  load  amplitude  equal  to  51  percent  of  the  ultimate  static  strength.  Panel 
I Cl  was  tested  under  Northrop' s  IRAD  program  to  determine  if  rounding  of 
stiffener  edges  is  effective  in  improving  the  panel  fatigue  life.  The 
stiffener  edges  were  rounded  to  eliminate  crack  initiation  adjacent  to  the 
stiffeners  due  to  web  bending  against  sharp  stiflener  edges  which  was 
observed  during  fatigue  tests  on  panels  MC3  and  MC4.  The  rounded  stiffener 
edges  were  effective  in  eliminating  crack  initiation  at  the  two  center 
stiffeners.  Creek  initiation  at  one  of  the  edge  stiffeners,  however, 
occurred  after  approximately  the  same  number  of  cycles  (40,000)  as  in  the 
case  of  panels  MC3  and  MC4  although  the  initial  crack  length  was  smaller 
(.15  inch  versus  2.5  inc'.us).  Panel  IC1  completed  100,000  cycles  of  fatigue 
loading  without  catast  ophic  failure.  During  these  100,000  cycles,  the 
crack  at  the  edge  stiffener  gt  w  to  a  length  of  4.75  inches.  The  panel  was 
residual  strength-tested  and  failed  by  stiffener  crippling  at  41,000  lbs. 
showing  no  reduction  from  the  ultimate  static  strength.  An  additional  ob¬ 
servation  from  these  tests  was  that  just  prior  to  failure,  a  sizeable  crack 
did  appear  at  one  of  the  center  stiffeners  indicating  that  rounded  stiffener 
edges  were  effective  in  reducing  the  web  bending  stresses  and  thus,  delay¬ 
ing  crack  initiation. 

Curved  composite  panels  CC3  and  CC4  were  tested  in  fatigue  with 
the  maximum  load  set  at  70%  of  the  average  static  strength  determined  from 
panels  CCl  and  CC2.  The  panels  failed  due  to  separation  of  the  stiffeners 
from  the  skin  after  61,640  and  12,758  cycles,  respectively.  In  the  case  of 
panel  CC3  a  single  stiffener  separated  whereas  for  panel  CC4  the  damage  was 
more  extensive  and  resulted  in  the  failure  of  three  stiffeners.  Panels  CC5 
and  CC6  were  fatigue-tested  at  60  percent  of  the  average  static  strength 
aud  both  survived  the  first  100,000  cycles  of  fatigue.  The  load  was  subse¬ 
quently  increased  to  65  percent  of  the  static  strength  and  the  fatigue 
test  continued.  Panel  CC5  survived  an  additional  100,000  cycles  at  this 
load  without  any  significant  loss  in  residual  strength.  Panel  CC6,  however, 
failed  after  approximately  6,000  cycles  at  the  increased  load.  The  results 
for  panels  CC3  through  CC6  are  discussed  in  Section  4. 
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3.5.3  Shear  Panel  Static  Tests 

The  test  results  for  curved  metal  and  composite  shear  panels 
are  summarized  in  Table  3.9.  The  lead  values  shown  in  the  table  are  the 
loads  applied  by  the  torque/cylinder.  Metal  panel  static  failure  was  due  to 
permanent  buckling  of  the  web.  The  failure,  load  was  well  above  the  design 
ultimate  and  the  panel  was  able  to  sustain  loads  higher  than  the  loads  at 
which  permanent  set  in  the  web  occurred.  The  diagonal  buckle  pattern 
representative  of  out-of -plane  skin  displacements  for  the  metal  shear 
panels  is  shown  in  Figure  3.28. 

The  composite  shear  panels  CS1  and  CS2  failed  at  nearly  identical 
loads  by  stiffener  web  separation.  Panel  CS1  failed  at  a  torque  cylinder 
load  of  16,000  lbs.  A  photograph  of  the  failed  panel  is  shown  in  Figure  3-29 
with  a  close-up  view  of  the  failure  area  shown  in  Figure  3.30.  The  buckle 
pattern  for  composite  shear  panels  was  similar  to  that  shown  in  Figure  3.28. 


Due  to  the  nature  of  the  test  set  up  for  the  shear  panels 
direct  measurement  of  the  shear  flow  Nxy  in  the  test  area  is  not  possible. 
The  cross-sectional  area  of  the  torque  box  cannot  be  directly  measured  due 
to  the  presence  of  attachment  hardware  at  the  corners.  In  addition, 
friction  inherent  in  the  test  arrangement  means  that  all  the  applied  torque 
is  not  converted  to  shear  flow  in  the  test  panel.  Hence,  the  first  metal 
shear  panel  static  test  was  used  to  calibrate  the  applied  torque  to  panel 
shear  flow  relationship.  Since  the  properties  of  7075-T6  aluminum  are 
well  established,  this  calibration  can  be  carried  out  by  plotting  the 
measured  shear  strain  in  the  panel  web  prior  to  web  buckling  versus  the 
applied  torque  as  shown  in  Figure  3.31.  Using  this  plot,  the  calibration 
proceeds  as  follows: 


T  74  Pm 
M  -  _ =  _ i. 

*y  2A  2A 


t 


where,  T  is  the  applied  torque,  A  the  cross-sectional  area  of  the  torque 

tube,  P,j,  the  applied  cylinder  load,  the  torque  arm  is  74  inches,  t  is  skin 

thickness  and  y  is  the  measured  shear  strain  in  the  web.  Frou  the 
xy 

above  equation  _ 


A 


f* 
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Figure  3.27.  Fatigue  Cracks  in  Metal  Compression  Panel  MC4 


Figure  3.28.  Diagonal  Buckle  Pattern  in  Metal  Shear  Panels 


Figure  3.30.  Close-Up  of  Composite  Shear  Panel  Failure  by  Stiffener/Web  Separation 
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or,  using  the  data  shown  in  Figure  3.31, 

1  AT 

A  =  2Gt  Ay 

Thus  the  cross-sectional  area  was  determined  from  Figure  3.31  to 

2 

be  462  in  .  This  yields  the  following  conversion  from  cylinder  load  to 
panel  shear  flow 

Nxy  =  0.08  PT 

where,  Nxy  is  in  lb/in  and  P^  in  lb.  The  strain  data  for  all  statically 
tested  shear  panels  is  given  in  Appendix  A. 


3.5.4  Shear  Panel  Fatigue  Tests 

The  shear  panel  fatigue  test  data  are  summarized  in  Table  3.10. 
Metal  panels  MS 3  and  MS4  were  tested  at  a  maximum  fatigue  load  equal  to 
80  percent  of  the  design  ultimate  strength  since  the  wide  scatter  in  the 
static  strength  of  MSI  and  MS2  made  it  difficult  to  define  a  meaningful 
average  static  strength. 

In  panel  MS3  test  the  bolt  in  the  corner  of  the  test  bay  failed 
after  7,200  cycles.  Fatigue  cracks  grew  soon  after  the  fastener  failure  in 
the  panel  web  adjacent  to  the  comer  hole  as  well  as  at  the  corner  hole.  The 
panel  failed  due  to  web  rupture  after  8,700  cycles.  The  rupture  was  caused 
by  cracks  which  grew  normal  to  the  direction  of  diagonal  tension. 

In  panel  MS4  the  corner  fastener  which  failed  during  the 
specimen  MS3  tests  was  replaced  by  a  higher  strength  fastener.  Panel  MS< 
sustained  12,500  cycles  before  any  cracks  were  visible.  The  first  crack 
appeared  in  the  stringer  near  the  intersection  of  the  frames,  where 
significant  stringer  bending  was  observed.  The  crack  was  located  in  the 
stiffener  heel  and  ran  along  the  length.  A  second  crack  developed  in  the 
panel  web  due  to  bending  of  the  web  after  14,800  cycles.  The  location  of 
this  crack  was  at  the  point  where  the  buckles  stop  in  the  diagonal  tension 
corner.  This  crack  resulted  in  panel  failure  after  additional  1200  cycles. 
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Another  crack  was  observed  after  15,000  cycles  along  the  stiffener  edge 
similar  to  the  cracks  in  the  compression  panels.  Complete  web  rupture 
occurred  at  about  16,600  cycles. 

Metal  shear  panels  MS5  and  MS 6  were  subjected  to  constant 
amplitude  fatigue  with  R  =  0  and  at  a  maximum  fatigue  load  equal  to  45  percent 
of  the  average  static  strength  obtained  from  static  tests  on  panels  MSI  and 
MS2.  The  absolute  fatigue  shear  load  amplitude  was  8,000  lbs.  as  compared 
with  10,000  lbs.  for  panels  MS3  and  MS4.  The  edges  of  panels  MS5  and  MS6 
were  thickened  by  means  of  an  externally  (convex  side)  bonded  aluminum 
doubler  to  facilitate  load  introduction  and  avoid  premature  cracking  of  the 
web  at  the  panel  corners.  Addition  of  the  doubler  did  not  affect  the  panel 
response  as  manifested  by  the  initial  buckling  loads  of  2,894  lbs.  and  3,811 
lbs.  for  panels  MSS  and  MS6,  respectively.  As  shown  in  Table  3.10  these 
buckling  loads  are  in  the  same  range  as  those  for  panels  MSI  through  MS4 
where  bonded  edge  doublers  were  not  used. 

Panels  MS5  and  MS6  showed  very  similar  modes  of  failure  although 
there  was  some  scatter  in  their  fatigue  life.  In  panel  MS5,  web  cracks  were 
first  observed  at  37,000  cycles  parallel  to  the  mid-bay  longitudinal  stiffe¬ 
ner  and  transverse  to  the  tension  field  at  the  vertical  frame  fastener  holes 
as  illustrated  in  Figure  3.32.  Final  failure  of  panel  MS5  occurred  by  web 
rupture  due  to  propagation  of  a  dominant  crack  in  a  direction  transverse 
to  the  tension  field  after  52,000  cycles.  The  dominant  crack  initiated  in 
the  lower  bay  and  propagated  across  the  second  bay  as  shown  in  Figure  3.32. 
Fatigue  cracks  in  panel  MS6  were  first  observed  at  22,000  cycles  and  final 
rupture  of  the  web  occurred  at  29,000  cycles. 

Panels  CS3  and  CS4  were  tested  under  constant  amplitude  fatigue 
load-ng  (R*-0.25)  .  The  maximum  fatigue  load  was  set  equal  to  60%  of  the 
static  ultimate  failure  load.  The  panels  completed  100,000  cycles  of 
fatigue  loading  without  any  detectable  damage.  The  panels  failed  at  loads 
of  15,200  and  16,400  lbs.,  respectively,  during  the  residual  strength  test. 

The  failure  load  and  mode  of  failure  were  almost  identical  to  the  static 
tests. 
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Figure  3.32.  Fatigue  Failure  Mode  of  Aluminum  Shear  Panel  MS5 
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The  maximum  load  amplitude  for  panels  CS5  and  CS6  was  increased 
to  70%  of  the  static  ultimate  load.  Panels  were  loaded  under  constant 
amplitude  fatigue  loading  (R  *  -0.25).  Both  panels  failed  during  fatigue 
after  72,000  cycles  and  19,500  cycles,  respectively.  The  failures  were 
sudden  and  not  detected  during  the  inspection  periods  which  were  set  to 
10,000  cycles.  Thus  the  time  of  damage  initiation  to  final  failure  was 
quite  rapid.  The  fatigue  failures  at  such  high  load  magnitudes  are  to  be 
expected  and  the  data  continue  to  demonstrate  excellent  durability  of 
composite  shear  panels. 

Panel  CS7  was  subjected  to  constant  amplitude  fatigue  loading 
(R  =  +0.25).  No  load  reversal  for  this  panel  was  chosen.  The  maximum  load 
amplitude  was  set  at  70  percent  of  the  static  ultimate  load.  The  panel 
developed  stiffener /web  separation* over  a  0.5-inch  length  after  7,000 
cycles.  This  disband  grew  rapidly  resulting  in  total  panel  failure  after 
8,700  cycles.  The  slightly  lower  life  with  less  severe  fatigue  loading 
points  towards  a  weakly-bonded  region  at  the  stiffener /skin  interface. 

Panel  CL*  was  tested  at  a  fatigue  load  amplitude  equal  to  70% 
of  the  average  static  strength  measured  for  similar  panels*  The  panel 
survived  100,000  cycles  of  fatigue  loading  without  any  detectable  damage. 

The  residual  strength  of  panel  CS8  was  14,500  lbs.  which  is  10%  lower  than 
the  average  static  strength.  The  failure  mode  of  the  panel  was  identical 
to  that  seen  in  the  static  tests  and  was  by  separation  at  the  stringer/web 
interface.  This  test  was  identical  in  all  respects  to  that  conducted  on 
panel  CS7. 

The  fatigue  test  results  for  the  two  panels  show  significant 
scatter  although  the  failure  modes  are  identical.  However,  the  reason  for 
this  scatter  can  be  deduced  from  the  observation  that  panel  CS7  developed 
stiffener/web  separation  very  early  in  the  fatigue  test  at  approximately 
7,000  cycles  and  thereafter  this  disbond  propagated  very  rapidly.  In  the 
case  of  panel  CS8,  no  such  stiffener /web  separation  was  observed  even  after 
100,000  cycles.  This  suggests  that  in  panel  CS7  a  flaw  in  the  form  of  a 
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void  or  incomplete  resin  cure  at  the  stiffener/web  interface  was  inadvertently 
introduced  during  fabrication  and  was  the  direct  cause  of  the  reduced  fatigue 
life. 

The  last  set  of  curved  composite  shear  panels  (CS9  and  CSJO) 
were  also  tested  under  constant  amplitude  fatigue  loading  at  an  R-ratio  of 
0.25.  The  fatigue  load  amplitude  for  these  two  panels,  however,  was  60%  of 
the  average  static  strength  for  similar  panels.  Panel  CS9  completed  100,000 
cycles  of  fatigue  and  survived.  The  residual  strength  of  this  panel  was 
not  significantly  reduced  from  the  average  static  strength.  Panel  CS10 
failed  by  stiffener/web  separation  after  69,100  cycles  of  fatigue.  The 
scatter  in  the  fatigue  test  results  for  these  two  panels  is  suspected  to  be 
due  to  variability  in  fabrication  of  the  panels. 

The  significance  of  the  test  results  presented  in  this  section 
is  discussed  in  Section  4. 
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SECTION  4 

DEVELOPMENT  OF  ANALYSIS  METHODS 

4.1  INTRODUCTION 

The  semi  empirical  analysis  of  postbuckled  panels  presented  in 
Section  2  was  done  in  steps  due  to  the  complexity  of  the  solution  process. 
These  analyses,  although  useful  for  design  purposes,  do  not  yield  the  stress 
and  displacement  fields  required  for  a  fatigue  analysis  of  the  panels. 
Furthermore,  with  the  semiempirical  analyses,  the  complex  mechanisms  of  load 
transfer  before  and  after  buckling  such  as  the  interaction  forces  between 
the  stiffeners  and  the  web,  and  the  effect  of  stiffness  changes  in  the  post- 
buckling  range,  cannot  be  modelled.  Most  of  these  analysis  methods  were 
developed  before  the  advent  of  high-speed  computers  and  result  In  desigi.s 
that  are  conservative  by  as  much  as  50  percent.  The  current  technology  and 
recent  advances  in  computation  methods  make  it  feasible  to  model  all  sig¬ 
nificant  aspects  of  postbuckled  panels  in  a  single  analysis  with  increased 
accuracy,  thereby  reducing  the  conservatism  in  the  final  design.  Thus,  de¬ 
velopment  of  such  an  analysis  to  accurately  predict  the  postbuckling  be¬ 
havior  of  the  panel  as  a  whole,  including  the  web  and  the  stiffeners,  was 
undertaken. 

It  was  envisioned  that  the  analysis  methodology  would  be  used 
as  a  design  tool.  Therefore,  ease  of  application  and  low  computational 
costs  were  prime  considerations  in  its  development.  Based  on  the  survey 
of  nonempi- .cal  static  analysis  methods  presented  in  Section  2,  the  total 
postbuckling  behavior  of  stiffened  panels  loaded  in  compression  or  in  shear 
was  modelled  using  the  principle  of  minimum  potential  energy.  The  analysis 
methodology  for  each  loading  case  is  generic  in  that  it  applies  to  curved, 
flat,  metal  and  composite  panels.  Details  of  the  analysis  for  compression 
and  shear  panels  are  discussed  in  the  following  paragraphs. 

4.2  COMPRESSION  PANEL  ANALYSIS 

The  governing  equations  in  this  formulation  were  derived  for 
orthotropic  laminates  that  are  balanced  and  symmetric.  It  is  assumed  that 
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the  webs  between  adjacent  stringei's  deform  in  an  identical  fashion.  The 
analysis  is  also  applicable  to  isotropic  materials  provided  the  appropriate 
constitutive  relations  are  used. 

A . 2 . 1  Geometry  and  Boundary  Conditions 

The  relationship  between  the  compression  panel  configuration  tested 
in  the  program  and  the  panel  geometry  used  in  the  analysis  is  shown  in  Figure 
4.1.  Since  adjacent  bays  are  assumed  to  deform  in  an  identical  fashion,  a 
single  bay  was  analyzed.  As  shown  in  Figure  4.1*  one  curved  edge  of  the  panel 
is  assumed  to  be  fixed  with  compression  load  applied  to  the  other  curved  edge 
At  the  loaded  edge,  due  to  the  presence  of  the  stiff  frames,  the  displacement 
in  the  Z-direction  (w)  is  restricted  to  zero  and  no  displacement  is  permitted 
in  the  Y-direction  (v  =  0  at  this  edge) .  The  lattar  boundary  condition  is  a 
realistic  representation  of  the  edge  conditions  in  panels  under  pure  compres¬ 
sion.  The  stringer  deformations  determine  the  boundary  conditions  at  the 
straight  edges  of  the  panel.  The  stringers  are  assumed  to  be  initially  straight 
and  deformation  relative  to  the  mid-surface  of  the  panel  in  t?)e  Z-direction  (w) 
occurs  only  due  to  stiffener  crippling  or  Euler  buckling  of  the  panel  which 
results  in  catastrophic  failure.  The  u,  and  v  displacements  at  the  straight 
edges  are  not  restricted. 

4.2.2  Strain  Energy  Expressions 

The  analysis  employs  the  principle  of  minimum  potential  energy. 

The  method  is  ideally  suited  for  this  analysis  since  it  simplifies  the  handling 
of  discrete  stiffeners  and  imposition  of  the  boundary  conditions. 

According  to  the  principle  of  minimum  potential  energy,  a  solution 
of  the  present  problem  renders  the  total  potential  a  relative  minimum.  The 
total  potential  energy,  II,  is  the  sum  of  the  strain  energy  stored  in  the  web, 

Uw,  in  the  stringers,  Ug,  in  the  frame,  U^,  and  the  potential  of  the  external 
load,  ft,  i.e. 

II  =  U  +  U  +  Up  +  ft  (32) 

w  s  f 
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Curved  Compression  Panel  Geometry  and  Coordinate  System 


In  applying  the  principle  of  minimum  potential  energy,  the  first 
step  is  to  assume  kinematically  admissible  displacement  functions  with  unknown 
coefficients.  These  displacement  functions  are  selected  to  satisfy  the  geo¬ 
metric  boundary  conditions.  The  total  potential  energy  is  computed  next  as  a 
function  of  the  unknown  coefficients.  The  governing  equations  are  then  obtained 
in  terms  of  the  unknown  coefficients  by  minimizing  the  total  potential  energy 
with  respect  to  these  coefficients.  Finally,  the  set  of  non-linear  equations 
obtained  by  minimization  of  the  total  potential  energy  are  numerically  solved 
to  determine  the  displacement  coefficients. 

The  strain  energy  stored  in  the  skin  and  the  stiffeners  can  be 
expressed  in  terms  of  the  assumed  displacements  u,  v  and  w  using  nonlinear 
strain-displacement  relations  for  curved  laminates  as: 

Strain  Energy  of  the  Skin: 
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where,  £  and  n  are  normalized  variables  as  shown  in  Figure  4.1,  and  u,  v  and  w 

* 

are  displacements  in  the  x,  y  and  z  directions,  respectively.  A  are  the 

ft 

elements  of  the  laminate  axial  stiffness  matrix  and  are  the  elements  of  the 
flexural  stiffness  matrix.  Commas  denote  differentiation  with  respect  to  the 
subscripted  variables.  The  strain  energy  due  to  shear  and  stretching  coupling 
is  ignored  since  the  laminate  is  balanced  and  symmetric  about  the  midplane, 
i.e.,  Alfi  =  A26  =  0. 


Strain  energy  in  the  stringer: 
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where,  Ag  is  the  cross-sectional  area  of  the  stringer,  Es  is  the  modulus  and 
I  is  the  equivalent  moment  of  inertia  of  the  stringer  about  the  z-axis. 


Strain  energy  in  the  frames: 
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where,  A^  is  the  cross-sectional  area  of  the  frame,  E^  is  the  average  modulus 
and  1^  is  the  equivalent  moment  of  inertia  of  the  frame  about  the  z-axis. 
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The  potential  of  the  external  load  is  expressed  as: 
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(36) 


where,  P  is  a  reference  load  and  X  is  an  unknown  load  control  parameter. 

4.2.3  Governing  Equations 

The  governing  equations  are  obtained  by  substituting  a  set  of  as¬ 
sumed  displacement  functions  with  unknown  coefficients  into  the  expressions 
for  the  total  potential  II  and  then  minimizing  it  with  respect  to  the  unknown 
coefficients.  As  discussed  in  Section  4.2.1,  the  boundary  conditions  that  the 
assumed  displacement  functions  must  satisfy  are: 
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where,  u(x,y),  v(x,y)  and  w(x,y)  are  the  displacements  in  the  x,  y  and  z  direc¬ 
tions,  respectively.  The  stiffeners  along  the  straight  edges  are  restricted 
to  in-plane  deformations. 

The  admissible  displacement  functions  for  the  compression  panel 
were  assumed  in  the  following  form: 
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and 


f 


x  =  a£  y  =  br| 


The  total  potential  energy  for  the  compression  panels  takes  the 


form: 


H  -  II(A^  ,  ,  0^^  ,  a^ ,  b^  •  X) 


(39) 


Using  the  indicial  notation  where  summation  over  repeated  indices  i9  implied, 
the  total  potential  energy  can  be  written  as  follows: 
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functions  obtained  from  the  evaluation  of  integrals 
which  appear  in  the  total  potential  energy  expression. 
Detailed  expressions  are  given  in  Appendix  B. 


The  conditions  for  obtaining  the  values  of  the  coefficients  that 
minimize  the  total  potential  energy,  H,  are 
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After  performing  the  differentiation  operation,  the  following  five 
nonlinear  simultaneous  equations  are  obtained: 
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The  governing  equations  are  obtained  by  linearizing  Equations  42 
through  46  by  first  expressing  the  total  potential  energy  as: 
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where  q  is  the  vector  of  the  unknown  coefficients.  The  equilibrium  equations 
related  to  the  stationary  condition  of  the  total  potential  energy  are: 
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The  resulting  equations  are  a  system  of  nonlinear  coupled  algebraic  equations. 

Linearization  of  these  equations  is  required  in  order  to  ensure  a  systematic 

solution  scheme.  Expanding  these  equations  about  a  known  position  vector  q 

•»o 

and  retaining  the  linear  terms  in  the  expansion,  the  following  set  of  equations 
is  obtained. 
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Defining 
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<*q  ~  ax 


~  3q 


a<i 


AA+£(<}o;  Ao>‘0 

the  equation  above  becomes, 


KA^PAX-RC^.Xo) 


(49) 


(50) 


The  matrix  K  is  called  the  tangent  stiffness  matrix  and  it  provides  information 
about  the  stability  of  the  structure.  At  the  bifurcation  point  K  becomes  sing- 
ular.  R,  F  and  are  called  the  residual  internal  and  external  load  vectors, 
respectively.  The  full  set  of  equations  to  be  solved  numerically  then  appear 
as  follows; 
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The  submatrices  or  the  global  matrix  in  Equation  51  can  be  expressed  ao: 
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Numerical  solution  of  Equations  51  is  accomplished  using  the  method 
described  in  Section  A. 2. 4. 


4.2.4  Numerical  Solution  Procedure 

Among  the  variety  of  numerical  solution  schemes,  the  following 
method  described  in  Reference  94  is  used  in  the  calculation  of  the  fundamental 
and  the  post-bifurcation  paths. 


The  full  set  of  equations  to  be  solved  numerically  can  be  written 


in  the  form: 

K  Aq  -  PAA  =  -  R 

-**■  "*»  * 
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(76) 


154 


t 


I 

l 


I 


where  the  unit  vector  t  is  colinear  with  the  vector  r  and  is  defined  as: 


k-1 


!t=  r/|rl  j  and  r(k'1)=¥  A»<J> 

~  -  j  ~ 


J=l 


(77) 


(k-1) 

where,  Ax  a  {Aq,AX}  and  r  is  tha  vector  of  the  total  of  instantaneous 

displacements  at  the  (k-l)1"^  iteration. 

Equation  76  can  be  rewritten  in  a  condensed  fora  as: 

H  A  x  =  S  (78) 


The  matrix  H  is  constructed  by  addition  of  the  column  -  P  and  the  row  t  to 


the  K  matrix  as  shown  below: 
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If  one  of  the  components  q^  of  the  vector  g  is  selected  as  a  control  parameter. 


then  t  initially  will  be: 


t  —  {0,....  0,1, 


0) 


(80) 


where,  the  1  appears  at  the  qQ  position.  Assuming  qQ  =  0  as  the  initial  equi¬ 
librium  position  for  zero  applied  load  (X  =  0) ,  the  extended  residual  vector  is: 

S(o)  =  {0,  As) 

•w 

=  {q^k\o>  for  k>l  and 
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Equation  81  governs  the  equilibrium  on  both  the  fundamental  (pre¬ 
buckling)  and  the  post-buckling  (post-bifurcation)  paths.  The  solution  method 
outlined  above  can  be  used  for  determining  equilibrium  points  on  both  parts; 
however,  to  ensure  convergence  of  the  solution  procedure  to  equilibrium  points 
on  the  post-buckling  path,  an  orthogonality  condition  has  to  be  imposed.  This 
is  achieved,  as  illustrated  in  Figure  4.2,  by  requiring  that  the  vector,  tp, 
tangent  at  point  B  to  the  fundamental  path,  be  non-colinear  with  the  vector, 
t„,  tangent  at  point  B  to  the  post-buckling  path. 

I 


Figure  4.2.  Tangent  Vectors  tp  and  tg 

The  vectors  t_  and  t„  can  be  approximated  as: 
r  d 

tf  =  A  xm/|Axml  and  t  =  a{a+ut,}  (82) 

T 

where,  a  =  l/}tD  |,  p  =  -a  t*,  a  is  the  eigen-vector  of  the  matrix  K  at  point 

li  i 

B. 

The  imposition  of  the  above  condition  prevents  the  solution  from 
returning  to  the  fundamental  path  of  equilibrium  in  the  postbuckling  regime. 
Once  the  first  equilibrium  point  is  obtained  on  the  postbuckling  path,  the 
subsequent  increments  follow  the  procedure  outlined  for  the  fundamental  path. 
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The  above  procedure  is  used  to  first  determine  the  buckling  load 
and  the  corresponding  mode  shape*  In- the  postbuckling  solution  only  the  un¬ 
knowns  corresponding  to  the  buckle  mode  shape  are  retained  and  the  rest  are 
ignored . 

4. 2.5  Displacement  Predictions 

The  compression  panel  analysis  has  been  coded  in  program  COMPAN 
documented  in  Reference  91  along  with  the  user  instructions.  The  program  was 
used  to  analyze  the  metal  and  composite  compression  panels  described  in  Sec¬ 
tion  3.  Actual  program  t uns  for  these  two  panels  are  given  in  Reference  91. 

Solutions  for  the  metal  and  composite  panels  were  obtained  by  first 
calculating  the  buckling  load  and  mode  shape  and  then  retaining  in  the  assumed 
displacements  the  unknown  coefficients  corresponding  to  this  mode  shape  only. 
Thus,  the  number  of  unknowns  were  significantly  reduced  from  over  20  to  5. 

The  calculated  mode  shape  for  composite  panels  was  six  half  waves  along  the 
load  direction  -.ad  one  half  wave  transverse  to  the  load  direction.  The  metal 
panel  was  predicted  to  buckle  into  five  half  waves 'along  the  load  direction 
and  one  half  wave  transverse  to  the  load  direction. 

The  postbuckling  predictions  consisted  of  out-of-plane  displace¬ 
ments  as  a  function  of  the  applied  load,  the  end  shortening,  strains  in  the 
stringers,  and  membrane  strain  distribution  in  the  skin.  The  predicted  end¬ 
shortening  for  the  metal  and  composite  panels  is  shown  in  Figures  4.3  and  4.4, 
respectively.  These  predictions  illustrate  the  nonlinear  postbuckling  response 
of  the  compression  panels .  End-shortening  data  were  not  measured  in  the  test 
program,  therefore,  verification  was  not  possible.  Comparison  of  the  other 
predictions  with  test  data  is  carried  out  in  Section  5,  where  the  accuracy  of 
the  predictions  is  also  discussed. 

4.3  SHEAR  PANEL  ANALYSIS 


The  shear  panel  analysis  closely  followed  the  approach  of  Denke 
(Reference  58)  and  Kudva  (Reference  75).  In  the  latter  study,  Denke' s  analysis 
for  isotropic  flat  shear  panels  was  extended  to  anisotropic  flat  shear  panels. 
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APPLIED  COMPRESSION  LOAD,  lbs/: 


Figure  4.4.  End-Shortening  for  Composite  Compression  Panels 


The  present  analysis  is  a  further  improvement  on  the  Reference  75  analysis  in 
that  it  applies  to  curved  anisotropic  panels .  The  von-Karman  strain  displace¬ 
ment  relations  used  in  the  present  analysis  include  the  strain  term  that 
accounts  for  panel  curvature. 

4.3.1  Panel  Geometry  and  Assumptions 

The  panel  geometry  for  analysis  formulation  and  its  relationship 
to  the  three  stringer  panels  tested  in  the  experimental  program  is  shown  in 
Figure  4.5.  The  horizontal  and  vertical  stiffeners  correspond  to  the  stringers 
and  the  frames,  respectively. 

The  panel  thickness  is  small  compared  with  the  other  dimensions  of 
the  panel.  This  allows  the  use  of  von  Karman  plate  theory.  When  subjected  to 
shear  loads  above  the  critical  load,  a  stiffened  curved  panel  undergoes  large 
deflection  and,  to  account  for  this,  nonlinear  strain-displacement  relations 
are  considered. 


The  nonlinear  strain-displacement  relations  are: 


(83) 


and  z  is  the  distance  from  the  middle  surface.  Superscript  o  indicates  the 
mid-plane  and  comma  denotes  differentiation  with  respect  to  the  subscript. 

The  solution  method  employs  the  concept  of  principle  of  minimum  potential 
energy.  The  total  potential  energy,  II,  is  the  sum  of  the  strain  energy  stored 
in  the  web,  Uw>  in  the  stringers,  Ug,  in  the  frame,  Uf,  and  the  potential  of 
the  external  load,  Q.  i.e., 

n  =  Uw  +  ug  +  uf  +  R  (84) 
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4.5.  Curved  Shear  Panel  Geometry  and  Coordinate  System 


The  total  potential  is  evaluated  in  terms  of  the  unknown  in-plane 
and  out-of -plane  deformation  parameters.  A  two  term  trignometric  expression 
is  assumed  for  the  out-of-plane  displacement.  The  unknown  parameters  associated 
with  the  assumed  displacement  function  w(x,y)  given  in  Equation  85  below  are 
m,  n  and  f. 

w(x,y)  »  8f  cos  px  cos  qy  cos  (m  —  +  n  •£)  (85) 
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.  —  2ir  .  -  2-/r  . 

where  m  =  —  a  sina,  n  =  y-  b  cos  a 

'  o  o 


7T  TT 

—  and  q 


2b 


a  and  X^  are  the  diagonal  tension  angle  and  wave  length  of  the  buckles,  respec¬ 
tively.  f  is  the  amplitude  of  the  buckles.  A  giaphic  illustration  of  the  as¬ 
sumed  displacement  functions  is  shown  in  Figure  4.6.  The  additional  in-plane 
deformation  parameters  assumed  are  e^,  e^,  e^  and  Yq  which  are  strain  quanti¬ 
ties  with  physical  significance  as  illustrated  in  Figure  4.7.  These  strain 
quantities  are  a  result  of  the  following  integrations: 


(86) 


The  out-of-plane  displacement  w(x,y)  is  kinematically  admissible 
since  w(x,lb)  *  w  (+a,y)  =  0. 

-V 
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DISPLACEMENT  AT  y  *  b 


A. 6*  Assuaed  Out-of-Plane  Displacement  Parameters 
a,  X  and  f 


63 


In  addition  t-.o  the  above  assumed  displacements,  the  following  are 
assumed  regarding  the  in- plane  stress  resultants: 

(a)  is  a  constant  independent  of  x  and  y  and  is  equal  to  the 

externally  applied  load. 

(B)  is  a  function  of  y  only. 


(C)  Ny  is  a  function  of  x  only. 

The  latter  two  assumptions  ensure  that  the  in-plane  equilibrium 
equations  are  satisfied. 


4.3.2 


Strain  Energy  Expressions 


The  strain  energies  Uw,  U  %  and  and  the  potential  of  the  external 
forces  are  expressed  as  snown  below. 

Strain  energy  of  the  web: 
a  b 

,J„  ‘  I  /  /  Ay  “i  "j  +  DtJ  Ki  d*  “y  <87> 

-a  -b 

The  first  and  the  second  term  of  the  integrand  represent  the  stored 
strain  energy  associated  with  in-plant,  and  bending  deforaatioijs ,  respectively. 
The  constitutive  relations  for  the  panel  under  consideration  are: 

K  )■  fAUA12°1  \ 


*12  A22  ° 


v,y  +  w,y/2  +  w/R 


0  0  A6fi  +  +  w.  w, 


x  'y 


The  average  stress  resultants  N  and  N  are  defined  by  Equation 

x  _  y 

(89)  below.  The  average  shear  stress  resultant  N^y  equals  the  applied  shear 
stress  due  to  assumption  (A)  above. 
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Nx  =  A  Nx(y)  dy  and  Ny  =  fzf  Ny(x)  dx  (89) 

•'-b  J- a 

Using  the  average  stress  resultants,  the  constitutive  relations  can  be 
rewritten  as  follows: 


where , 
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The  expressions  for  N  (y)  and  N  (x)  can  be  obtained  by  performing  the  inte 

x  y 


gration  shown  below: 
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(92) 
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Knowing  the  explicit  form  of  N  (y)  and  N  (x) ,  the  strain  energy  due  to  in- 

x  y 

plane  deformation  can  be  obtained  from  the  following: 
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where  a^  =  .  Also,  since  the  laminate  is  symmetric  and  balanced  about 

the  midplane  A.,  and  A„,.  are  equal  to  zero, 
lb  zb 

Performing  the  integration  of  Equation  (93)  leads  to  the  follow¬ 
ing  expression: 


"i-2*"  [ei*uVife62+^u)] 


(94) 


where , 
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The  bending  strain  energy,  i.e.,  the  second  term  in  Equation  (87)  can  be 
rewritten  as: 
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(96) 


168 


Substituting  the  expression  for  w(x,y)  into  Equation  (96)  and  carrying  out 
the  integration  results  in: 
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(98) 


Strain  Energy  of  the  Stringers  and  Frames: 

The  axial  strain  energy  of  a  stringer  is  given  by: 

US  -  1  Vs  <2a>  4  (99) 

The  axial  strain  energy  of  a  frame  is  given  by: 

UF  =  \  Vf  (2b)  ev  (100) 

In  addition,  the  deflection  shape  of  the  stringer  is  assumed  in 
the  form  written  below: 

Sy  -  2ey  b  cos2  px  (101) 
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is: 


The  bending  strain  energy  for  this  deflection  shape  in  a  stringer 


> 

where , 


S 

E„ 


“S  EsIs  b2 

2a3  S  S 


e 


2 

y 


Cross-sectional  area  of  the  stringer 
Young's  modulus  of  the  stringer 
Cross-sectional  area  of  the  frame 
Young's  modulus  of  the  frame 
Moment  of  inertia  of  the  stringer. 


(102) 


Total  Potential  Energy: 

In  addition  to  the  strain  energy  contributions  from  the  various 
components,  the  potential  of  the  externally  applied  shear  load  needs  to  be 
taken  into  account.  This  is  expressed  as: 


Q  =  -N  Y  4  ab 
xy  o 


(103) 


where. 


N 


xy 

Then, 


=  applied  in-plane  shear  load  per  unit  length 
the  total  potential  II  is: 


n  =  uI  +  oB  +  us  +  usl3  +  uF  +  a 


004) 


4.3.3  Governing  Equations  and  Solution  Procedure 

The  governing  equations  are  obtained  by  minimizing  the  total  po¬ 
tential  with  respect  to  the  deformation  parameters.  Hence, 


M_  =  M  =  in_  =  iflaiiL  =  iIL  =  ilL_.o 

3y  rT;  9e  3e  3e 

o  r  dm  dn  h  v  y 


(105) 


17C 
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Equations  (106)  werei  numerically  solved  using  standard  mathemat¬ 
ical  library  routines.  The  specific  library  was  IMSL  and  the  solution  routine 
used  was  ZSPOW. 

4.3.4  Displacement  Predictions 

The  shear  panel  analysis  was  coded  in  a  computer  program  called 
SHRPAN1  which  is  documented  in  Reference  91  along  with  user  instructions. 

The  metal  and  composite  shear  panel  designs  of  Section  3  were  analyzed  using 
SHRPAN1.  Actual  program  runs  for  these  two  panels  are  given  in  Reference  91. 

The  shear  panel  solutions  provide  out-of-plane  displacements, 
skin  strains,  buckle  wavelength,  and  diagonal  tension  angle  as  a  function 
of  the  applied  loads.  The  solution  is  initiated  above  the  buckling  load 
at  some  pre-selected  value  of  the  displacement.  The  predicted  out-of-plane 
displacements  as  a  function  of  the  applied  load  are  shown  in  Figures  4.8  and 
4.9  for  the  metal  and  composite  panels,  respectively.  The  out-of -plane  dis¬ 
placement  contours  at  a  constant  load  for  the  metal  and  composite  panels  are 
shown  in  Figures  4.10  and  4.11,  respectively.  These  postbuckling  displace¬ 
ment  contours  illustrate  the  diagonal  buckling  pattern.  Verification  of  the 
out-of-plane  displacements  was  not  possible  since  these  displacements  were 
not  measured  during  the  tests.  However,  the  predicted  strains  and  the  diag¬ 
onal  tension  angle  are  compared  with  the  test  data  in  Section  5  where  the 
accuracy  of  the  solution  is  also  discussed. 

4.4  FATIGUE  ANALYSIS  APPROACH 

The  fatigue  tests  conducted  on  metal  and  composite  panels  in  this 
program  were  useful  in  identifying  the  panel  failure  modes.  Based  on  this 
evidence  and  on  test  data  from  other  sources  such  as  Reference  80,  approaches 
to  performing  fatigue  life  analyses  of  metal  and  composite  panels  were  de¬ 
veloped.  These  approaches  utilize  the  results  of  the  nonempirical  analysis 
developed  in  this  program  and  in  addition  require  some  fracture  property 
data.  These  fatigue  life  prediction  methodologies  are  outlined  in  the  fol¬ 
lowing  paragraphs. 
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Plane  Displacement  at  Panel  Center  as  a  Function  of  Applied  Shear 
>r  Composite  Panel.  SHRPANl  Predictions 


SHRPANl  Predictions 


4.4.1 


Metal  Panels 


The  fatigue  life  prediction  methodology  for  metal  panels  loaded 
in  compression  or  shear  is  summarized  in  Figure  4.12.  The  observed  failure 
modes  for  metal  compression  panels  are  cracks  in  the  skin  parallel  to  the 
stringers,  but  away  from  fastener  holes,  or  cracks  in  the  stringer  or  skin 
itself  at  fastener  holes  (see  Figure  2.2).  In  the  case  of  metal  shear 
panels,  the  dominant  cracks  are  those  that  initiate  in  the  skin  at  stiffener 
to  skin  attach  fastener  holes  and  propagate  transverse  to  the  diagonal  ten¬ 
sion  direction.  For  fatigue  failures  initiating  at  fastener  holes,  the 
analysis  approach  is  the  same  for  shear  and  compression  panels  with  the 
fatigue  life  being  governed  by  crack  growth  at  the  fastener  holes.  In  the 
case  of  skin  cracks  parallel  to  the  stringers  and  away  from  fastener  holes, 
the  analysis  approach,  is  somewhat  different  and  the  fatigue  life  is  governed 
by  crack  initiation  and  growth  in  the  skin. 

As  shown  in  Figure  4.12,  a  durability  rather  than  a  damage  toler¬ 
ance  approach  is  adopted  in  the  analysis  of  fatigue  failures  initiating  at 
fastener  holes.  An  initial  0.01-inch  corner  flaw  is  assumed  to  exist  at  the 
hole.  The  stress  intensity  factors  for  this  initial  flaw  are  computed  using 
available  analysis  methods.  Flaw  growth  in  shear  panels  occurs  due  to  diago¬ 
nal  tension  stresses  and,  therefore,  the  principal  tensile  stress  in  the 
skin  is  used  for  computing  the  stress  intensity  factor.  In  compression 
panels  transverse  tensile  stresses  are  caused  at  the  apex  of  the  fastener 
hole  by  the  remotely  applied  compression  stress  and  are  equal  to  it  in  mag¬ 
nitude.  This  stress  is  used  in  computing  the  stress  intensity  factor  for 
the  initial  flaw. 

Once  the  stress  intensity  factors  have  been  determined,  the 
Forman  crack  growth  equation  (Reference  95)  can  be  used  to  determine  the 
crack  growth  life  for  the  metal  panels. 

In  the  case  of  skin  cracks  parallel  to  the  stringers,  the  local 
skin  stresses  have  to  be  computed  using  the  nonempirical  analysis  described 
in  the  preceeding  paragraphs.  The  fatigue  life  is  estimated  as  the  sum  of 
crack  initiation  life  and  the  crack  growth  life.  The  crack  initiation  life 
(to  0.01")  can  be  predicted  using  the  cumulative  damage  analysis  of 
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COMPRESSION 


Figure  4 .12.  Life  Prediction  Methodology  for  Postbuckled  Metal  Panels 


Reference  96.  Crack  growth  life  can  be  predicted  using  the  Forman  equation 
(Reference  95) .  A  bending  correction  may  have  to  be  used  to  obtain  the 
effective  stress  intensity  factor. 

4.4.2  Composite  Panels 

The  predictive  methodology  for  composite  panels  is  more  complex 
due  to  the  lack  of  a  static  analysis  method  that  can  predict  the  stresses 
at  the  stiffener/skin  interface.  The  life  prediction  approach  for  composite 
panels  shown  in  Figure  4.13  addresses  the  stiffener/web  separation  mode  of 
failure  observed  in  panel  fatigue  failures  by  way  of  the  nonempirical  static 
analysis  developed  in  this  program  and  simple  beam  model  of  the  stiffener/ 
web  attach  area. 

The  approach  shown  in  Figure  4.13  applies  to  both  shear  and  com¬ 
pression  panels  and  is  based  on  computing  the  strain  energy  release  rate  at 
the  stiffener/web  interface.  In  this  approach  the  nonempirical  analysis  is 
used  to  determine  the  maximum  skin  deflection  in  the  postbuckling  regime  due 
to  fatigue.  This  maximum  deflection  is  then  applied  to  the  beam  model  shown 
which  represents  a  strip  between  the  stiffener  and  the  center  of  the  skin. 

An  initial  delamination  of  length  bo  is  assumed  to  exist  at  the  interface  of 
the  doubler  and  the  skin  laminate  in  the  beam  model.  The  applied  maximum 
displacement  is  then  used  to  compute  the  strain  energy  release  rate  at  the 
top  of  the  initial  delamination.  This  strain  energy  release  rate  G  is  the 
driving  force  for  delamination  growth  and  can  be  used  to  predict  delamination 
growth  under  static  as  well  as  fatigue  loading.  Fatigue  life  prediction  is 
accomplished  by  using  the  nonlinear  growth  law. 

i  '  CWG  -  ioth,n 

where,  C,  n  and  AG^  are  constants  determined  from  fatigue  tests  on  simple 
specimens  as  shown  in  the  box  on  the  extreme  right  hand  side  of  Figure  4.13. 
Static  failure  prediction  requires  a  knowledge  of  the  critical  strain  energy 
release  rate  which  can  be  determined  from  static  tests  on  specimens  iden¬ 
tical  to  the  fatigue  test  specimens. 

The  essential  requirements  for  this  life  prediction  approach  are 
the  development  of  an  analysis  for  the  beam  model  and  the  appropriate  fracture 
properties.  These  developments  are  needed  to  establish  a  fatigue  life  pre¬ 
diction  methodology  lor  postbuckled  composite  panels. 
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.13  Life  Prediction  Methodology  for  Postbuckled  Composite  Panels 


SECTION  5 


DISCUSSION  OF  RESULTS 


5.1  INTRODUCTION 

The  static  and  fatigue  test  data  presented  in  Section  3  and  Appendix 
B  were  analyzed  to  correlate  the  measured  initial  buckling  and  ultimate  strength, 
and  strain  values  with  predictions.  The  static  test  data  were  also  used  to 
determine  panel  stiffness  change  due  to  postbuckling.  The  fatigue  life  data 
were  utilized  to  establish  S-N  curves  for  metal  and  composite  panels.  These 
results  are  discussed  in  the  following  paragraphs. 

5.2  CURVED  PANELS  UNDER  STATIC  COMPRESSION  LOAD 


Metal  compression  panel  static  test  results  and  their  correlation 
with  the  semiempirical  analysis  given  in  Section  2  are  summarized  in  Table 
5.1.  In  this  table,  the  skin  buckling  data  obtained  from  the  fatigue  test 
specimens  are  pooled  together  with  the  static  test  data.  As  can  be  seen  from 
Table  5.1  the  local  skin  buckling  predictions  which  were  based  on  the  use  of 
the  empirical  correction  factor  Kc  in  the  expression 

K  7T2E 

F  =  - 

cr  12(1-V2) 

with  determined  from  Figure  2.6,  are  quite  conservative.  The  modified 
values  of  Kc  account  for  the  imperfection  sensitivity  of  curved  panels.  A 
teexamination  of  the  equivalent  Kc  for  the  data  in  Table  5.1  showed  that  the 
present  data  are  much  closer  to  predictions  based  on  the  theoretical  value 
of  30.  Thus,  the  theoretical  curve  shown  in  Figure  2.6  is  reasonably  accu¬ 
rate  for  use  in  predicting  the  local  skin  buckling  loads  and  strains.  The 
comparison  of  the  theoretical  and  the  semiempirical  predictions  is  shown  in 
Figure  5.1.  The  semiempirical  predictions  take  into  account  the  influence 
of  r/t  ratio  on  the  buckling  load.  Metal  panel  failure  under  static  load 
occurred  primarily  due  to  stiffener  crippling  and  as  a  consequence  the  fail¬ 
ure  data  shown  in  Figure  5.1  agree  very  well  with  predictions  which  were 
based  on  the  calculated  stiffener  crippling  loads. 
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TABLE  6.1.  CURVED  METAL  COMPRESSION  PANEL  STATIC  DATA 


IRAD  program 


Figure  5.1.  Comparison  of  Metal  Compression  Panel  Test  Data  with  Predictions 


The  skin  buckling  loads  were  determined  using  mid-bay  back-to-back 
strain  gage  data.  Figure  5.2  shows  the  typical  response  of  these  back-to-back 
strain  gages.  These  gages  were  located  at  the  panel  centroid  and  illustrate 
that  there  was  a  change  in  the  buckle  wave  length  at  approximately  30  kips. 

This  wave  length  change  in  the  postbuckled  regime  is  also  corroborated  by  the 
out-of-plane  displacement  data  shovm  in  Figure  5.3.  The  actual  variation  in 
the  number  of  skin  buckles  indicated  by  the  strain  and  displacement  data  depends 
on  the  location  of  these  gages  and,  therefore,  measurement  at  a  single  location 
is  not  sufficient  to  fully  describe  the  buckle  pattern  progression  with  load. 

For  instance,  the  displacement  plot  for  6^  in  Figure  5.3  shows  that  the  buckle 
pattern  changed  thrice  prior  to  attaining  the  final  configuration  shown  in 
Figure  3.23e.  In  Figure  5.3  the  out-of-plane  displacements  predicted  using 
program  COMP AN  are  also  shown  for  comparison.  The  displacement  corresponding 
to  in  the  figure  is  in  reasonable  agreement  with  the  test  data  in  the  low 
postbuckling  range.  However,  with  increasing  load  the  test  data  indicate 
changes  in  buckle  mode  shapes  which  are  not  accounted  for  in  the  analysis. 

Thus,  the  discrepancy  between  test  data  and  predictions  is  significant.  For 
displacements  at  an<*  63  which  were  symmetrically  located  with  respect  to 
panel  centerline,  the  disparity  in  the  test  data  and  predictions  is  signifi¬ 
cant.  In  addition  to  changes  in  the  buckle  mode  shapes,  another  reason  for 
this  discrepancy  could  be  the  extreme  sensitivity  of  the  out-of-plane  dis¬ 
placements  to  variations  in  the  location  of  measurement  on  the  panel.  From 
these  comparisons  it  is  apparent  that  the  compression  panel  analysis  for  metal 
panels  needs  additional  refinement. 

The  axial  strain  in  the  stiffeners  was  approximately  bilinear  up  to 
failure  with  a  distinct  change  in  the  slope  after  skin  buckling.  Axial  strain 
variation  with  applie  1  load  for  the  four  stringers  of  metal  panel  MCI  is  shown 
in  Figure  5.4.  The  average  of  the  four  stringer  strains  was  used  in  determin¬ 
ing  the  panel  gross  stiffness  change  due  to  skin  buckling.  Panel  MCI  showed  a 
postbuckled  stiffness  that  was  64  percent  of  the  prebuckling  stiffness,  whereas, 
for  panel  MC2  this  number  was  56  percent. 

Correlation  of  buckling  and  failure  strain  data  for  composite  com¬ 
pression  panels  with  theoretical  and  semiempirical  predictions  is  shown  in 
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Figure  5.2.  Metal  Compression  Panel  Web  Static  Strain  Response 


COMPRESSION  PANEL  MCI 


■Plane  Displacements  for  Metal  Compression  Panels 


Stiffener  Axial  Strains  for  Metal  Compression  Panels 


Table  5.2  and  Figure  5.5.  The  measured  web  buckling  strains  are  considerably 
lower  than  the  predicted  values.  This,  however,  was  not  surprising  since  the 
effective  width  of  the  web  was  not  known  a  priori  and  was  assumed  as  the  dis¬ 
tance  between  the  adjacent  stiffener  flanges  as  shown  by  "CURRENT"  in  Figure 
5.5.  The  test  data  indicate  that  this  assumption  is  unconservative.  An 
effective  width  equal  to  the  distance  between  the  centers  of  adjacent  stiff¬ 
ener  flanges,  indicated  by  "PROPOSED"  in  Figure  5.5,  yields  excellent  corre¬ 
lation  between  the  predictions  and  the  test  data. 

Composite  compression  panel  strength  measurements  closely  agree  with 
semiempirical  predictions  based  on  the  stiffener  crippling  mode  of  failure. 

Thus,  the  stiffener /web  separation  mode  of  failure  seen  in  the  static  tests 
is  induced  by  stiffener  crippling. 

Figure  5.6  shows  the  back-to-back  mid-bay  strain  gage  response  for 
composite  panel  CC1  and  is  typical  of  that  observed  for  the  other  five  panels. 

The  buckle  pattern  c'>  inges  shown  are  significant  from  the  point  of  view  of 
developing  a  nonempirical  analysis.  The  regions  where  the  buckle  pattern  changes 
are  apt  to  cause  numerical  difficulties  in  predicting  the  postbuckled  response 
and  the  nonlinear  analysis  may  have  to  be  performed  piecewise  with  the  regions 
selected  so  that  no  change  occurs  in  the  buckle  pattern.  The  out-of-plane 
displacement  data  for  composite  panel  CC1  are  shown  in  Figure  5.7.  These  data 
corroborate  the  buckle  pattern  progression  indicated  by  the  strain  gages.  For 
comparison,  the  predicted  values  of  out-of-plane  displacements  are  also  shown 
in  Figure  5.7.  The  trends  in  the  predicted  displacements  at  62  and  6^  match 
the  test  data.  However,  numerically  there  is  a  significant  amount  of  dis¬ 
crepancy. 


The  stiffener  axial  strain  response  for  composite  panel  CC1  is  shown 
in  Figure  5.8  and  is  representative  of  that  seen  for  the  remaining  composite 
compression  panels.  In  Figure  5,8,  COMP AN  predictions  of  stringer  strains  are 
also  shown.  As  opposed  to  the  skin  displacements  the  stringer  strain  variation 
with  applied  loads  is  reasonably  well  matched.  The  strain  values  predicted, 
however,  are  considerably  higher  than  the  test  data.  The  stiffener  strain 
data  were  used  to  determine  the  panel  stiffness  changes  due  to  buckling  shown 
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TABLE  5.2.  COMPOSITE  COMPRESSION  PANEL  STATIC  TEST  DATA  CORRELATION  WITH  PREDICTIONS 


**  Panels  CC3  through  CC6  were  fatigue  test 


i 


1 


* 


”a/d3  MO  ^3/a33 
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gure  5.5.  Comparison  of  Composite  Compression  Panel  Test  Data  with  Predictions 


6000 
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Figure  5.6.  Composite  Compression  Panel  CC1  Web  Strain  Response 


COMPOSITE  COMPRESSION  PANEL  CC1 


-Plane  Displacements  for  Composite  Compression  PaneL  CC1 


Figure  5.8.  Stiffener  Axial  Strains  for  Composite  Panel  CC1 


in  Table  5.2.  The  reduction  in  stiffness  after  buckling  for  composite  panels 
is  similar  to  that  for  metal  panels  and  is  approximately  40  percent. 

5.3  CURVED  PANELS  UNDER  COMPRESSION  FATIGUE  LOADING 


The  fatigue  data  for  metal  and  composite  compression  panels  are  shown 
in  the  S-N  diagram  of  Figure  5.9.  The  curves  were  faired  to  represent  the  data 
trend  and  due  to  the  limited  number  of  data  points  a  definitive  threshold  for 
100,000  cycles  of  constant  amplitude  fatigue  cannot  be  established.  The  data 
for  composite  compression  panels,  however,  are  consistant  with  those  obtained 
from  other  tests  and  summarized  in  Figure  2.4.  Therefore,  composite  compres¬ 
sion  panel  fatigue  does  not  appear  to  be  a  concern  in  the  2500-3500  yin/in 
operating  strain  level  typically  seen  in  postbuckled  structures. 

The  metal  panel  fatigue  tests  were  useful  in  identifying  the  failure 
mode  that  needs  to  be  accounted  for  in  developing  a  fatigue  life  prediction 
methodology  (see  Section  4.4).  The  fatigue  data  also  show  that  the  metal 
panels  are  quite  sensitive  to  fatigue  and,  as  illustrated  in  Figure  5.10,  are 
inferior  to  composite  panels  designed  to  the  same  loading  conditions.  Addi¬ 
tional  metal  panel  tests,  however,  should  be  conducted  to  accurately  define 
their  S-N  response. 

The  periodic  strain  surveys  conducted  during  the  fatigue  tests  were 
used  to  determine  if  repeated  buckling  of  the  panels  influenced  the  initial 
buckling  load  or  panel  stiffness.  The  data  showed  that  repeated  loading  did 
not  influence  the  initial  buckling  load  or  the  panel  stiffness. 

5.4  CURVED  PANELS  UNDER  SHEAR  LOADING 


The  metal  shear  panel  test  data  analysis  is  summarized  in  Table  5.3. 
Comparison  of  the  data  with  predictions  is  shown  in  Figure  5.11.  The  predic¬ 
tions  were  based  on  Figure  3.5  (Reference  34)  assuming  the  skin  width  bw  equal 
to  the  stiffener  pitch.  The  data  show  that  with  this  definition  of  the  skin 
width,  the  buckling  load  predictions  are  reasonable  estimates  considering 
the  scatter  in  the  initial  buckling  loads ,  and  the  semiempirical  approach 
of  Reference  34  can  be  readily  used  for  design  purposes.  The  failure  load 
predicted  for  metal  shear  panels  (Nu^  =  900  lbs/in)  was  based  on  the 
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Comparison  of  S-N  Data  for  Metal  and  Composite  Compression  Panels 


Figure  5.10.  Normalized  Fatigue  Response  of  Compression  Panels 


TABLE  5.3.  METAL  SHEAR  PANEL  STATIC  TEST  DATA  CORRELATION  WITH  PREDICTIONS 
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Predictions  based  on  stiffener  forced  crippling  mode  of  failure 
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stringer  forced  crippling  mode  of  failure.  The  ultimate  shear  flow  calculated 
for  shear  panels  MSI  and  MS2  using  web-maximum  shear  strain  data  at  failure, 
compares  well  with  these  predictions.  A  comparison  of  the  measured  stringer 
axial  strains  with  that  predicted  for  stringer  forced  crippling,  however,  shows 
that  the  measured  strains  are  considerably  less  than  the  predicted  crippling 
values.  The  actual  failure  of  the  metal  shear  panels  was  by  permanent  set  in 
the  skin.  This  was  confirmed  by  an  examination  of  the  measured  principal 
strains  (Gage  5  in  Figure  3.13)  in  the  skin  which  showed  that  the  6500  ]iin/in 
yield  strain  for  7075-T6  aluminum  had  been  exceeded.  Thus,  it  is  not  surpris¬ 
ing  that  the  measured  stringer  axial  strains  did  not  exceed  the  predictions  for 
forced  crippling. 

In  order  to  verify  if  permanent  set  can  be  predicted  using  the  flat 
metal  panel  criterion,  the  allowable  diagonal  tension  factor  for  the  present 
panels  was  calculated  with  the  aid  of  Equation  31 .  The  maximum  value  of  k  was 
calculated  to  be  0.554  which  translated  into  an  ultimate  load  to  initial  buck¬ 
ling  load  ratio  of  2.2.  However,  the  data  show  that  for  these  panels,  the 
ratio  is  of  the  order  of  4.  Thus,  the  flat  metal  panel  permanent  set  criterion 
is  very  conservative  for  curved  panels.  A  criterion  needs  to  be  developed 
for  curved  panels  by  additional  testing. 

The  web  shear  strain  variation  with  applied  load  is  shown  in  Figure 
5.12.  These  data  were  used  to  compute  the  change  in  panel  stiffness  after  buck¬ 
ling.  As  shown  in  Table  5.3,  the  metal  panels  retain  a  large  percentage  of 
their  initial  stiffness  in  the  postbuckling  range.  The  maximum  reduction  in 
stiffness  was  approximately  17  percent.  The  stringer  and  ring  axial  strain 
variations  with  applied  load  are  shown  in  Figures  5.13  and  5.14.  For  comparison, 
predictions  from  SHRPAN1  are  also  shown.  The  reasons  for  the  large  discrepancies 
are  explained  in  Section  5.6.  The  measured  diagonal  tension  angle  was  approxi¬ 
mately  27°  and  is  less  than  that  predicted  by  tension  field  theory  but  is 
Larger  than  the  20°  angle  predicted  by  SHRPANl.  A  plot  of  the  diagonal  ten¬ 
sion  angle,  calculated  using  the  mid-bay  strain  rosettes  (Gages  2  and  5), 
versus  the  applied  cylinder  load  is  shown  in  Figure  5.15,  along  with  the 
predictions  from  SHRPANl.  The  predicted  values  agree  reasonably  well  with 
the  test  data. 
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Figure  5.15.  Diagonal  Tension  Angle  as  a  Function  of  Applied  Load  for  Metal  Shear  Panels 


The  composite  shear  panel  test  data  analysis  is  summarized  in  Table 
5.4.  Comparison  of  the  initial  buckling  and  failure  strain  data  with  predic¬ 
tions  is  shown  in  'Figure  5.16.  The  predictions  were  based  on  a  skin  width 
equal  to  the  stringer  spacing  and  are  unconservative.  In  the  case  of  composite 
shear  panels,  a  more  realistic  definition  of  the  skin  width  is  the  distance 
between  adjacent  stiffener  flanges.  Use  of  this  definition  yields  better 
correlation  between  the  test  data  and  the  predictions.  Failure  of  the  composite 
panels  was  predicted  by  forced  crippling  of  the  rings.  In  Section  2  it  was  also 
noted  that  forced  crippling  strains  for  the  stiffeners  correspond  to  stiffener/ 
skin  separation  strains.  The  data  comparison  in  Figure  5.16  along  with  the  ring/ 
web  separation  mode  of  failure  observed  in  composite  shear  panels  CS1  and  CS2 
substantiate  the  above  hypothesis.  The  failure  predictions  made  using  the 
modified  tension  field  theory  are  on  the  conservative  side  by  approximately 
35  percent.  The  measured  maximum  shear  strains  in  the  composite  panel  web 
shown  in  Figure  5.17  were  used  to  determine  panel  stiffness  change  in  the  post- 
buckling  regime.  The  composite  shear  panels  show  a  dramatic  loss  in  stiffness 
after  initial  buckling  of  the  skin.  As  indicated  in  Table  5.4  the  postbuckled 
stiffness  for  these  panels  is  approximately  45  percent  of  the  initial  stiffness. 
These  data  are  of  significance  in  the  design  of  postbuckled  composite  panels 
since  the  stiffness  has  a  direct  influence  on  the  aeroelastic  response  of  the 
panels.  Therefore,  for  composite  shear  panels  verification  of  the  design  for 
aeroelastic  response  criteria  will  be  essential. 

The  hat  section  stringers  showed  significant  bending  during  the  static 
tests.  The  back-to-back  strain  gages  were  placed  on  the  stiffener  skin  flange 
and  the  crown  flange  and  due  to  the  local  bending  of  the  crown  flange,  separa¬ 
tion  of  axial  and  bending  strains  for  the  hat  section  stringers  was  not  possible. 
In  Figure  5.18  and  5.19,  the  stringer  and  ring  strains  measured  on  the  skin 
flange  of  these  stiffeners  are  shown.  These  data  do  not  indicate  the  true 
axial  strains  due  to  the  reasons  cited  above;  however,  they  do  show  the  buildup 
of  strain  in  the  stiffeners  with  increasing  load  after  diagonal  buckling.  The 
stiffener  strain  predictions  shown  in  Figure  5.18  were  obtained  from  SHRPANl. 

The  diagonal  tension  angle  variation  with  load  for  the  composite  panels 
is  shown  in  Figure  5.20.  The  predictions  are  in  reasonable  agreement  with  test 
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**  Lowest  shear  strain  as  measured  by  gages  2  and  5  in  Figure  3.13. 
t  Quantities  with  subscript  ’O'  denote  predicted  values, 

Y  =  770  yin/ in,  Yn  =  3500  yin/in,  (e  )  =  -2718  yin/in,  (e  )n  =  -3192  yin/in 

cr  vj  or  u  OS  u 


000 


Figure  5.17.  Composite  Panel  Web  Maximum  Shear  Strains  as  a  Function  of  Applied  Load 
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Strains  for  Composite  Shear  Panels 


Figure  4.19.  Ring  Strains  for 


Figure  5.20.  Diagonal  Tension  Angle  as  a  Function  of  Applied  Load  for  Composite  Shear  Panels 


data.  The  discontinuity  in  predictions  at  an  applied  cylinder  load  of  5,000 
lbs  is  due  to  the  instability  of  SHRPAN1  solution  in  a  close  neighborhood  of 
the  buckling  load.  As  the  load  increases,  the  diagonal  tension  angle 
decreases,  indicating  the  tendency  of  the  diagonal  buckles  to  merge  and  cross 
over  the  stringer  into  the  adjacent  bay.  After  buckling  the  diagonal  tension 
angle  is  approximately  18  degrees  and  compares  favorably  with  the  predictions. 

5.5  CURVED  SHEAR  PANELS  UNDER  FATIGUE  LOADING 


The  fatigue  test  data  for  metal  and  composite  shear  panels  are  shown 
plotted  as  a  function  of  the  maximum  fatigue  load  normalized  to  their  respec¬ 
tive  static  strengths  in  Figure  5.21.  The  data  demonstrate  the  sensitivity  of 
metal  shear  panels  to  fatigue  loading  and  that  their  response  is  much  inferior 
to  that  for  composite  panels  designed  for  the  same  loading  condition.  The  com 
posite  shear  panel  fatigue  response  is  not  affected  by  the  partial  reversal  of 
the  shear  loading  as  is  evidenced  by  a  comparison  of  R  =  +0.25  and  R  =  -0.25 
fatigue  data  in  Figure  5.21. 

The  fatigue  tests  were  also  useful  in  identifying  metal  and  composite 
panel  failure  modes.  Crack  initiation  and  propagation  in  the  metal  panels  as 
shown  in  Figure  3.32  and  stiffener /web  separation  in  composite  panels  are  the 
critical  modes  that  are  addressed  in  the  fatigue  analysis  methodology  proposed 
in  Section  4.4. 


Analysis  of  the  periodic  strain  survey  data  for  the  metal  panels  did 
not  show  any  influence  of  repeated  loading  on  initial  buckling  loads  or  panel 
stiffness  in  the  postbuckled  range.  A  majority  of  the  relevant  gages  on  the 
composite  shear  panels  were  lost  due  to  fatigue  damage  and  the  data  obtained 
could  not  be  used  for  a  meaningful  interpretation  of  repeated  buckling  effects 
The  buckling  loads  measured  in  these  strain  surveys,  however,  do  not  show  any 
influence  of  repeated  loading  on  panel  stiffnesses. 


5.6 


DISCUSSION  OF  ANALYSIS  AND  TEST  DATA  CORRELATION 


The  out-of-plane  displacement  and  strain  predictions  from  programs 
COMPAN  for  compression  panels  and  SHRPAN1  for  shear  panels  are  significantly 
different  from  the  measured  values. 


211 


l n 

CNi 

LH 

• 

CSJ 

CD 

1 

CD 

II 

II 

oc 

ac 

o 

•v 

II 

CO 

CO 

_ 1 

_ i 

QC 

UJ 

UJ 

<c 

N 

CO 

CL¬ 

CL¬ 

-J 

UJ 

UJ 

UJ 

z: 

h- 

H- 

>  « 

»■  < 

Q_ 

CO 

CO 

o 

o 

— t 

Q- 

Q_ 

<C 

sr 

SI 

h— 

Q 

o 

UJ 

CD 

UJ 

s: 

◄  < 


l n 

CD 


*2, 


CO 

UJ 

CD 

>- 

CD 

LL. 

O 

cn 


% 


CD 


cd 


H19N3U1S  3UV1S/QV01  3n9IlVd  XVU 


fn 


212 


igure  5,21.  Normalized  Fatigue  Response  Metal  and  Composite  of  Shear  Panels 


There  are  several  reasons  for  these  discrepancies.  First  of  all  the 
results  of  the  energy  based  analyses  are  very  sensitive  to  the  assumed  dis¬ 
placement  functions.  Secondly  to  account  for  mode  shape  changes  in  the  post- 
buckling  path  either  the  assumed  displacements  should  be  a  superposition  of 
several  modes  or  if  a  limited  number  of  terms  are  used,  the  mode  shape  para¬ 
meters  should  be  treated  as  unknowns.  The  latter  choice  is  more  desirable 
since  assumed  displacements  with  a  large  number  of  unknowns  inevitably  lead  to 
numerical  difficulties  in  the  solution  of  the  resulting  nonlinear  equations. 

Several  numerical  problems  were  also  encountered  in  solving  for  the 
metal  and  composite  compression  and  shear  panels.  In  the  case  of  compression 
panels,  the  solution  procedure  required  that  the  starting  point  be  zero  load 
and  zero  displacement.  The  solution  then  progressed  by  marching  up  the  pre¬ 
buckling  path  to  a  load  value  slightly  less  than  the  buckling  load  and  then 
switching  over  to  the  postbuckling  path  with  a  mathematical  artifice  of  ortho¬ 
gonal  vectors.  The  multivalued  nature  of  the  postbuckling  path  at  bifurcation 
of  equilibrium  near  the  buckling  load  leads  to  nonconvergent  solutions  due  to 
numerical  oscillations.  Thus,  the  orthogonal  vector  approach  works  only  for 
the  simplest  of  assumed  displacement  functions.  In  order  to  circumvent  these 
problems,  a  solution  scheme  that  starts  from  an  initially  guessed  pair  of  load 
and  displacement  values  in  the  vicinity  of  the  postbuckling  path,  but  suffic¬ 
iently  greater  than  the  load  and  displacement  values  at  bifurcation,  should  be 
utilized.  In  the  shear  panel  analysis  this  latter  scheme  was  adopted  with  the 
buckling  loads  being  computed  externally  by  programs  such  as  SS8 .  The  solution, 
however,  was  extremely  sensitive  to  the  initially  guessed  values  for  the  un¬ 
knowns  and  considerable  expertise  was  required  to  select  initial  values  that 
led  to  converged  solutions.  One  possible  strategy  that  may  be  used  in  se¬ 
lecting  the  initial  displacements  is  given  in  Reference  91, 

Based  on  the  present  experience,  it  is  recommended  that  in  future 
attempts  to  model  the  postbuckling  behavior  of  compression  and  shear  panels 
the  following  techniques  be  used: 

(A)  Introduce  Initial  imperfection  or  a  transverse  load  to  elim¬ 
inate  numerical  problems  in  the  vicinity  of  the  bifurcation 
point. 
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(B)  Treat  mode  shape  or  buckle  wave  lenom  a  =  , 

oucKie  wave  length  as  an  unknown  in  the 

assumed  displacements. 

(C)  Evaluate  functions  other  than  trignometric  functions  fo;  the 
assumed  displacements . 
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SECTION  6 


CONCLUSIONS 

The  significant  conclusions  from  this  program  are  summarized  in 
the  following  paragraphs. 

6.1  DESIGN  AND  ANALYSIS  METHODOLOGY  FOR  CURVED  COMPOSITE  POSTBUCKLED 

PANELS 

1 .  A  semiempirical  static  analysis  methodology  was  developed  for 
curved  composite  panels  loaded  in  compression  or  in  shear. 

2.  Experimental  verification  data  provided  guidelines  for  deter¬ 
mining  the  skin  dimensions  to  be  used  in  calculating  initial 
buckling  loads  for  both  compression  and  shear  panels. 

3.  Ultimate  load  predictions  based  on  the  semiempirical  analysis 
for  compression  panels  are  very  accurate  and  can  be  readily 
used  for  design  purposes. 

4.  The  modified  tension  field  theory  is  applicable  to  curved 
composite  shear  panels.  Ultimate  load  predictions  are  con¬ 
servative  by  approximately  35  percent. 

5.  The  postbuckled  stiffness  of  compression  panels  is  decreased 
by  approximately  40  percent  from  the  initial  stiffness. 

6.  In  shear  panels  the  loss  in  stiffness  after  buckling  is 
approximately  55  percent. 


7.  Stiffener/web  separation  was  the  observed  failure  mode  for 
static  and  fatigue  loading  of  compression  and  shear  panels. 


8.  An  empirical  equation  was  developed  to  predict  the  static 
stiffener/web  separation  strain  for  shear  panels. 

9.  Stringer  or  ring  forced  crippling  strains  correspond  to  the 
stiffener/web  separation  strain  for  shear  panels. 

10.  Fatigue  loading  is  not  a  concern  for  compression  or  shear 
panels  and  sufficient  data  exist  to  determine  their  safe 
operating  strain  levels. 

11.  Based  on  the  fatigue  failure  modes  observed  in  the  tests, 
an  approach  to  life  prediction  for  compression  and  shear 
panels  was  developed. 

12.  Repeated  buckling  had  no  influence  on  panel  initial  buckling. 

13.  Non-empirical  Rayleigh-Ritz  analyses  of  postbuckled  com¬ 
pression  and  shear  panels  have  been  developed.  The  analyses 
although  capable  of  predicting  the  detailed  displacement  and 
stress  field  in  postbuckled  panels  require  further  refinement 
to  ensure  numerical  accuracy. 

14.  The  program  results  were  used  to  develop  a  design  guide  for 
compression  and  shear  panels. 

6.2  DESIGN  AND  ANALYSIS  METHODOLOGY  FOR  CURVED  METAL  POSTBUCKLED 

PANELS  ~  ~  "  ~  . 

1.  Applicability  of  the  tension  field  theory  to  curved  shear 
panels  was  verified  by  test  data.  Skin  permanent  set  was 
seen  to  be  the  primary  failure  mode  in  these  panels.  A 
need  to  obtain  a  permanent  set  criterion  for  curved  panels 
was  identified. 

2.  Ultimate  load  predictions  based  on  the  available  analysis 
methods  were  found  to  be  quite  accurate  for  compression  panels. 
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3.  The  postbuckled  stiffness  reduction  for  compression  panels 
is  the  same  as  that  for  composite  panels .  In  the  case  of 
shear  panels  the  stiffness  change  after  buckling  was  seen  to 
be  only  about  15  percent. 

4.  Fatigue  sensitivity  of  compression  and  shear  panels  was  found 
to  be  the  greatest  concern. 

5.  Skin  cracking  parallel  to  the  stringers  and  away  from  fastener 
holes  was  identified  as  a  failure  mode  in  the  compression 
panels. 

6.  Skin  cracks  originating  at  stiffener  attach  fastener  holes 
and  propagating  transverse  to  the  tension  field  direction 
was  identified  as  the  failure  mode  in  shear  panels. 

7.  Based  on  these  failure  modes,  a  fatigue  life  prediction 
approach  was  formulated  for  compression  and  shear  panels. 

8.  Repeated  buckling  had  no  influence  on  the  initial  buckling 
load  or  stiffness  of  compression  or  shear  panels. 

9.  The  program  results  were  used  to  develop  a  design  guide. 

6.3  RECOMMENDATIONS  FOR  FUTURE  WORK 

1.  Develop  a  design  guide  for  panels  under  combined  load. 

2.  Complete  the  development  of  the  fatigue  life  prediction  meth¬ 
odologies  for  metal  and  composite  panels  and  extend  the  meth¬ 
odology  to  panels  operating  under  combined  loads. 
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APPENDIX  A 


COMPRESSION  AND  SHEAR  PANEL  STRAIN  DATA 

A.l  COMPRESSION  PANEL  STRAIN  DATA 

The  strain  data  obtained  from  all  compression  panel  static 
tests  and  fatigue  strain  surveys  are  tabulated  in  this  section  of  the  appen¬ 
dix.  Correspondence  Table  A-l  should  be  used  to  correlate  the  gage  numbers 
in  the  strain  data  tables  with  the  locations  shown  in  Figure  3.12.  All  gages 
were  axial  gages. 


» 

TABLE  A -1.  GAGE  NUMBER  CORRESPONDENCE  TABLE  FOR  STATIC 
COMPRESSION  TEST  PANELS 


GAGE  NO.  IN  FIGURE  3.12* 

GAGE  NO.  IN  DATA  TABLE* 

IE 

1 

11 

2 

2E 

3 

21 

4 

3E 

5 

31 

6 

4E 

7 

41 

8 

5E 

9 

51 

10 

6E 

11 

61 

12 

7E 

13 

71 

14 

8E 

15 

81 

16 

9E 

17 

91 

18 

10E 

19 

101 

20 

11E 

21 

111 

22 

12E** 

23 

121** 

24 

*Fatigue  panel  MC3  instrumented  with  gages  1  through  8. 

Fatigue  panels  MC4,  CC3  and  CC4  instrumented  with  gages  1  through  6. 
Fatigue  panels  CC5  and  CC6  instrumented  with  gages  1  through  4  and  6 
Gage  6  data  corresponds  to  column  numbers  9  and  10  in  the  data  table 


**Gage  on  midbay  stringer  of  panel  CC2. 
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SHEAR  PANEL  STRAIN  DATA 


uj  i  i  - 4iumrjthoJ<r 
<3  1  *  *  » 

n  •  •  . . .  .  •  .  - 

*  5»oot5(ywM^«n^NU}nnpH(fl^n<du)(i)tfN 

fu5c5ra^^i2><^cohj-*T^;'*r,im.^'cr>*ir>ajd»«3 
•H^wr)nt  ru  cu  g3  0o4i*j>c»c*p'co^  <c 

w  -  •  •  •  * . . * 

w  i  i  i^c«r-^*^r^r~^uT’M\ji7)wnIicouiC75»-'',*r«^occ 
<|  i  i  i  i  t  t  ^rotuprWnpcnpMWcu  i 


w  i  t  »  i  i  »  ;  i  r*  “ 


-«ir-pr^tr,r,,^^0oor*  Q^cncnn 
r^vttptoCoN^cooicQpjfy^ca 
i  rWmt*  o>v ©os wr-H r»  C5 


©®^  u>gw©*r 

p*  «^Tconfu©S©  «/> 

—4  •-*  «-*  *-4(\> 


T’r'tr-r^cs-'%'u^^«g>rrr“riV-rnr^toa)U)<s>c^corN-cn 
i  i  i  -rarr  ©r-cn<*>©foo»ift»-»»u)6ar-r'^» 

t  )  H»4««4  I  H  OJ  CO  r»  *4  «H  ■•*♦  l/^  I  CO 


■  »/>o>pruc$oppc^cQ  Tfncoc^r^  to-4coc* 
~  o  u .  ri  if>  rj  *r  B  *■  f*-  tft  **  ®  cit  w  «-« 

^cvjmrnpc>mr)r>t*ir>ri?u  i 


^  vwinfnm<»o»ou)covnt  vr^-f'ioor-**  *-» 
to  roco«x>iv  f^  tor- <?>  retain  njoovrtcftcoO’tf-ra 

*3  '^•^•ror.^uu/)C»,ircofvimcDrrin-Nu3'r'C^,,e 

<r  ^'Hjttwajnnnnnw.''  f 


u j  «  uJcr>r-rumn>oo»«.0’rtJsfoa>a>c^if;Po^r^r*-0« 

§  1,1,1 


s  ui  •« -J  J7)  ur>  v  **>  ru  0*  «qr> 
>  r>  *c  <v  t  -v  ^  ^ 


i  »  »  i  \  i  i 


! 


to 

UJ 


©nnnofw 
I  IA  COl/l 
I  I  ^fU 
I  « 


(UWV'f 
•  III 


U5h 

s" 


c 

■H 

<1" 


s 

PS 

< 

W 

P 

O 

PS 

o 

H 


C 

•H 

A 


25 
* — 1 


25 

M 

2 

H 

CO 


XS^^u)rum</jyjtorn?nryrur-«>c3P)(jDr-ioo*r-pvnnJ 
uj  i  i  ■^rQr\ja5roajr<tDfio^po56cy--‘b>NQ>aj£»'<3-^ 

g  iiit  -*^rj)lua  (jjajw^4 


rr  io  » fu  <•*  tfi  •«  -«  ©  to  •«  u 

»4  )  H  »4  »4  03  Q 

I  I  1  I  I  1 


r  <*  ft!  <u  ru  to  «p  «»  ^  ru  a  •*  o>  o  n  £■  (*}  to  -*  t  » 
^Prunja3cv3r>rWKi«-«^^v9P  u  fA  u>  •»*  to  r- 

-•inv*»»CD«&to,,-c>tt)£»r-^~« 
i  i  i  ~4-«(\jruru(u«*  i 


a 

<E 

£ 

*< 

S2 


j  i  i 


-«  V  w  oi  in  in  to  r>  fu 
♦  fl&r-  ^,rr-n‘vr^«coiA -rty  A 
i  fuii)  co  ru  i  hJtuwp^Ww®© 


i©npr-^evcor)f 
j  ~'~3mc0for 


stSftiN  {xvra  for  op^cif^N  nc*  m-ter  ss, ess  cvcics  o?  »‘<vn«u£ 


Srswt/>aja*meor*Mj^ 
aj*\j  t  10  «>r-  ~*tuu> 
*<ft)wtiftU)n®wnn 
I  I  S  I  I  I  *4w*w**+ 


O’ttacocar-tnrJtD^^ 
w  w  w>  uj  r- «  w  oj  to  ^  tn 
^CUD’TL’M^OI^IXl'Tin 

<  )  I  9  S  I  I 

9  t  I  • 


e>c*toruLT>^*4»-*^fv 
a  r*  u>oj  cn  ^  ir>  id 
i  **nirnTi/>«fie-wec«9 
i  i  i  i  »  .  *  i  i  <* 
• 


©^mu)i/>J3«dC»rntoco 
•cfucv)'e  i/>ws»<»  ^rurn 
i  i  i  i  •  •  i 

«  i  i  i 


t>  to  e»  *r  •«*  cn  cj  m  e  &  r*  n 

**«<a  •^pufu^i/i^ed’Vf^ 
-n..^vyHdc»-Mru^t® 

111*1*1  «-e***«1 
till 


r^c&orufjp^fuf^mom 

r-  <u  v  r-  a  r-  rr  e»i/>  -* 
-•w  n  <*•  t/>  r-  <ft -* ro  ^  tfi 

1  I  «  J  3  1  I  •*< 

I  »  »  I 


w  fur- t-  c*-  w  <v  ®  w  <* 

HWrJ^®lft»tC9P)0 

r*fyn^WV0«©^DlO 

*11* 


to  (\i  ^  oj  o  *jj  m  r- •«»  <°>  fn 
cu  <xi  r>  o  co  r- to  «*  to  «■*•  to 
ni  n  «r  u>  r- « <u  *r  Ui 

i  i  i  *  I  I  I 

•  ill 


r-caco©r>r-toto<oco® 
w  <u  ?n  <?  tr.  co  c*  ®  c\j 

i  i  i  i  i  t  i  f 

t  i  < 


oaj^toco©ru’«rtooo©nj 

C  r<rt^r<HW(U 


STRAIN  OftTfl  FOR  SPtCIfCH  RC4  ftFTEB  i*® ,9&S  CVfct.ES  OF  FfitlCUE 


u>  ru noi^  a  cv  in 
4»4  iw  r-y»^rnwr>ojM^>^  «r 

rJ  •♦futnn^iftu'^c»p«onH.r« 
«  t  f  I  I  I  I  -rH  W  fa  (\J  PJ  I  I  t 
*»  *  I  «  I  ?  I 


>)oj  ^vp-ciwoc.  i-no:nf.c» 
W  n^'r  wps>oc»^efj--«^G3na 
W  W<^ClVU)W'3SOV  OVlftMD 
«  OO  (  I  «  I  !  o4^^H 

<5  t  i 


uvpwu>vnu)>«)©LJ'innT 
iou)i/>r*fM-u>»tf>«x\n.«>ccr‘ 
^•f^n  ^  w>  u> » ru -r  ia  oo  ^  « 

i  i  i  t  i  •  i  *4+-'  <«**«w«n 
t  i  i  i  i  i  i  i 


i^^T'w©o^-^»ow>©r)c\im 
©<»_*cor-c>»v*~trjcuu-)*p'vco 
•^cyni^»^ru»ims-viPC3»€do‘®v^r 
i  i  i  i  i  i  i  i  i  i  »rurum 
i  i  i  :  i 


^mnninujwfTflujwnco^ 
t  wnTui«3cca^wf»)Lnvwn 
i  I  i  i  i  i  i  ****a»i\jfu 
*  i  i  i  i  i  i 


t^vino»cow>©nx>'ri/)»p*gn 

©,-,,,^^T?G>tjOLnCOm  *■  «r  iam 

w  v  *  *  utta  o  *-*  m  u>  co  oo  «  in 
I  t  t  i  I 

;  i  i  <  i  t  i  i  i 


<U'V»Ar-CTLrHS>CDVfVtpru^Lf>Wf» 

1111(1  hhhwh^NM 

I  I  I  •  I  *  i  »  I 


3Q*’Ct®(\i<"«r-~aMs'C»r-<0<0o<£* 
•  ®0>®»hOl^©tflC9^C6^n 

**•  nj  cu  ▼  *r  19  o>  <s»  a*  r>  v  r-  r-  no  tn 

I  I  I  |  *  I  |  Rrt  *1 

I  I  I  I  I  I  1  I 


»-4foa)coo*-«r-tra><D<3-«coLf) 
i  Air.7iA0}®«4(mA(oS^N(u 
i  i  s  ;  i  ^HHHMraajwn 
i  i  i  i  i  c  i  t  i 

W<^fUiO®  <*»;*»  r-  -wn^w 
^njajfuajfnr-  Qsrr^csc^otiniD 
-««f>vwrvW©^(\jn(6ovN 
i  i  i  t  i  i  i 

•  iiiiiii 


0(UY<oa®(uvu)ooo(\uoa>,r® 


,\.  2  SHEAR  PANEL  STRAIN  DATA 


The  strain  data  obtained  from  all  shear  panel  static  tests  and 
fatigue  strain  surveys  are  tabulated  in  this  appendix.  The  nomenclature 
for  the  strain  gages  in  the  following  tables  differs  from  that  shown  in 
Figure  3.12.  The  correspondence  Table  A-2  below  should  be  used  to  correlate 
the  gage  numbers  in  the  strain  data  tables  with  the  locations  shown  in 
Figure  3.13. 
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TAI'.l  h  A- 2.  SHEAR  i’ANI'T,  CACli  CORRIOSI'ONDIONCIO  TAlii.lO. 
SHOWN  IN  FIGURE  3.13. 


CACKS  OR  I  KNTIOI)  AS 


NO.  IN  FIGURE  3.13 
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TYPE  OF  GAGE 
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Axial 
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Axial 


Axial 


Axial 


Axia] 


Axial 


Axial 


Axial 


Axial 


Axial 


Axial 


Axial 


Shear 


GAGE  NO.  IN  DATA  TABLES 


1 

2 

3 


A 


Shear 


Shear 

Shear 


NOTES:  (1)  Gages  8  through  13  omitted  from  fatigue  test  panels  and 

gages  5  and  7  assigned  numbers  27  through  34  in  the  data 
tab  tes . 

(2)  After  50,000  cycles  of  fatigue  on  specimen  CS4,  gages  310, 
31  and  the  first  kg  of  2E  wore  lost.  All  remaining  strain 
channels  numbered  onserutively . 
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